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PREFACE 


The  theory  of  turbulence  has  made  so  much  progress  during  these  last 
years  that  it  is  of  interest  to  state  exactly  the  obtained  results ,  the 
hypotheses  on  which  these  results  are  based,  and  the  directions  in  which 
new  research  is  being  conducted. 

Messrs.  Bass  and  Agostini  have  undertaken  this  work  and  have  summar- 
ized our  actual  knowledge  of  turbulence  in  a  series  of  conferences  which 
took  place  at  the  Sorbonne,  within  the  Paris  Institute  of  Mechanics.   The 
reader  will  find  in  the  following  pages  the  text  of  these  conferences, 
perfected  and  revised  by  Mr.  Bass. 

In  view  of  the  magnitude  of  the  subject  and  its  simultaneously  phy- 
sical and  theoretical  aspects  it  had  seemed  advisable  to  entrust  this 
work  to  a  team  formed  by  a  mathematician  and  a  physicist.  Mr.  Bass  had 
taken  the  responsibility  for  the  theoretical  part,  Mr.  Agostini  for  the 
physical  part. 

Initially,  the  report  was  intended  to  contain  three  theoretical  and 
two  physical  chapters  followed  by  a  chapter  on  the  technique  of  the 
measurements  and  on  the  appratus  used:   anemometers  and  statistical- 
measurement  apparatus . 

The  unexpected  death  of  Mr.  Agostini  in  August  19^9  unfortunately 
made  modifications  of  the  original  project  necessary.   This  premature 
death  deprived  us  of  a  highly  valuable  physicist,  at  the  peak  of  intel- 
lectual maturity,  whose  current  work  on  these  problems  showed  particu- 
larly remarkable  promise . 

Mr.  Agostini  had  only  just  begun  drawing  up  the  two  last  chapters; 
Mr.  Bass  had  to  take  up  the  editorial  work  and  to  complete  it  according 
to  Mr.  Agostini 's  notes.   The  chapter  on  the  technique  of  the  measurements 
has  been  omitted  and  will  form  the  object  of  a  later  publication. 

In  order  to  make  up  for  the  gap  in  the  experimental  part,  the  text 
was  supplemented  by  some  curves  furnished  by  Mr.  Favre  which  will  allow 
utilization  for  numerical  calculations  and  will  enable  the  reader  to 
judge  the  agreement  between  theory  and  tests. 


A.  Fortier, 
Professor  at  the  Sorbonne 


NACA  TM  1377 


TABLE  OF  CONTENTS 


INTRODUCTION 


CHAPTER  I 
GENERAL  ASPECTS  OF  TURBULENCE.   THE  STATISTICAL  METHOD 

1 .  Definition  of  turbulence   6 

2.  Average  values  -  statistics  8 

3.  Random  variables  and  the  laws  of  probability 10 

k.   The  concept  of  random  point  -  velocity 

field  -  turbulent  diffusion 13 

5.  Equations  of  development  of  the  laws  of  probability 15 

6.  Random  velocity.   Hydrodynamic  equations   25 

7-  Systems  of  molecules 29 

CHAPTER  II 
CORRELATIONS  AND  SPECTRAL  FUNCTIONS 

8.  Introduction  -  correlations  in 

space  -  homogeneity,  isotropy 33 

9.  Properties  of  the  functions  f,  g,  a,  b,  c. 

Incompressibility 37 

10.  Spectral  decomposition  of  the  velocity J+3 

11 .  Spectral  tensor  and  correlation  tensor 4  5 

12.  Spectral  tensor  of  isotropic,  incompressible  turbulence  ....  50 

13.  Energy  interpretation  of  the  spectral  function  F(k)  .  55 

Ik.   Relations  between  spectral  function  F(k) 

and  correlation  functions  f(r)  and  g(r) 58 

15.  Lateral  and  longitudinal  spectrum  62 

CHAPTER  III 
DYNAMICS  OF  TURBULENCE 

16.  Introduction 6k 

17.  Fundamental  equation  of  turbulent  dynamics   65 

18.  Case  of  isotropic  turbulence 68 

19.  Local  form  of  the  fundamental  equation 73 

20.  Solution  of  the  fundamental  equation,  when  the 

triple  correlations  are  disregarded  76 

21.  Solutions  involving  a  similarity  hypothesis  80 

ii 


NACA  TM  1377 


22.  Transformation  of  the  fundamental 

equation  in  spectral  terms  86 

23.  First  theory  of  Heisenberg  89 

2k.   First  theory  of  Heisenberg  (continued). 

Space-time  correlations   93 

25-  Second  theory  of  Heisenberg   99 


CHAPTER  IV 
THEORY  OF  LOCAL  ISOTROPY  AND  STATISTICAL  EQUILIBRIUM 

26.  Introduction 106 

27.  Definition  of  local  homogeneity  and  local  isotropy  108 

28.  Similarity  hypotheses.   Statistical  equilibrium   Ill 

29.  Case  of  high  Reynolds  numbers Ilk 

30.  Validity  of  the  similarity  laws 118 

31.  Interpretation  of  the  laws  of  statistical 

equilibrium  in  spectral  terms  -  Weizsacker's 

and  Heisenberg' s  theories   121 

CHAPTER  V 
DECAY  OF  THE  TURBULENCE  BEHIND  A  GRID 

32.  History   . 124 

33«  Initial  and  final  phase  of  turbulence   126 

34.  Concepts  regarding  the  structure  of  the 

final  phase  of  turbulence 132 

35«  The  concept  of  "dynamic  statistical  equilibrium"  135 

36.  Synthesis  of  the  results  relating  to  the 

structure  of  the  spectrum  of  turbulence 141 


INDEX  OF  PRINCIPAL  NOTATIONS,  FUNDAMENTAL  FORMULAS, 

AND  DIMENSIONAL  EQUATIONS 1^3 

APPENDIX li+7 

Some  experimental  results   1^7 

REFERENCES 150 


iii 


Digitized  by  the  Internet  Archive 

in  2011  with  funding  from 

University  of  Florida,  George  A.  Smathers  Libraries  with  support  from  LYRASIS  and  the  Sloan  Foundation 


http://www.archive.org/details/theoriesofturbulOOunit 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  MEMORANDUM  1377 


THE  THEORIES  OF  TURBULENCE* 
By  L.  Agostini  and  J.  Bass 

INTRODUCTION 


The  theory  of  turbulence  reached  its  full  growth  at  the  end  of  the 
19th  century  as  a  result  of  the  work  by  Boussinesq  (1877)  anci  Reynolds 
(l893)-   It  then  underwent  a  long  period  of  stagnation  which  ended  under 
the  impulse  given  to  it  by  the  development  of  wind  tunnels  caused  by  the 
needs  of  aviation.  Numerous  researchers,  mathematicians,  aerodynamic ists, 
and  meteorologists  attempted  to  put  Reynolds'  elementary  statistical 
theory  in  a  more  precise  form,  to  define  the  fundamental  quantities,  to 
set  up  the  equations  which  connect  them,  and  to  explain  the  peculiarities 
of  turbulent  flows.   This  second  period  of  the  science  of  turbulence  ended 
before  the  war  and  had  its  apotheosis  at  the  1938  Congress  of  Applied 
Mechanics . 

During  the  war,  some  isolated  scientists  -  von  Weizsacker  and 
Heisenberg  in  Germany,  KoLmogoroff  in  Russia,  Onsager  in  the  U.S.A.  - 
started  a  program  of  research  which  forms  the  third  period.  By  a  system 
of  assumptions  which  make  it  possible  to  approach  the  structure  of  tur- 
bulence in  well-defined  limiting  conditions  quantitatively,  they  obtained 
a  certain  number  of  laws  on  the  correlations  and  the  spectrum.   These 
results,  once  they  became  known,  caused  a  spate  of  new  researches,  the 
most  outstanding  of  which  are  those  by  the  team  Batchelor-Townsend  at 
Cambridge . 

The  analysis  of  these  works  became  the  subject  of  a  series  of  lec- 
tures at  the  Sorbonne  in  February-March  19^9>  which  subsequently  were 
edited  and  completed.   The  mathematical  theory  of  turbulence  had  already 
been  published  in  19^6  (ref .  3)  but  practically  ignored  all  publications 
later  than  19^-0.   Since  the  late  reports  have  Improved  the  mathematical 
language  of  turbulence,  it  was  deemed  advisable  to  start  with  a  detailed 
account  of  the  mathematical  methods  applicable  to  turbulence,  inspired 
at  first  by  the  work  of  the  French  school,  above  all  for  the  basic  prin- 
ciples, then  the  work  of  foreigners,  above  all  for  the  theory  of  the 
spectrum. 

_  _ 

Les  Theories  de  la  Turbulence.   Publications  Scientif iques  et 
Techniques  du  Ministere  de  L'Air,  No.  237,  1950. 
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The  first  chapter  deals  with  the  precise  and  elementary  definition 
of  turbulence  (sections  1  and  2)  and  describes  the  tools  of  mathematical 
statistics  on  which  the  ultimate  developments  are  based.   Starting  from 
paragraph  3,  chapter  I,  the  reader  should  be  familiar  with  the  defini- 
tions of  the  calculus  of  probabilities,  the  theorem  of  total  probabili- 
ties, and  the  theorem  of  compound  probabilities.   This  chapter  is  entirely 
theoretical,  and  its  aim  is  to  review  the  methods  suggested  by  the  theory 
of  random  functions  and  which  seem  likely  to  be  applied  to  turbulence . 
Only  the  use  of  Navier's  equations  has,  so  far,  produced  positive  results, 
and  chapters  III,  IV,  and  V  are  largely  devoted  to  it.   However,  it  should 
be  pointed  out  that  there  are  theories  less  familiar  to  hydrodynamicists 
which  have  been  proved  in  other  branches  of  physics  (kinetic  theory 
of  gases,  quantum  mechanics).  The  purely  random  method  is  described  in 
paragraph  5>  and  its  adaptation  to  molecular  systems  (according  to  Born 
and  Green)  in  paragraph  7-   The  statistical  method  has  the  advantage  of 
furnishing  a  remarkable  demonstration  of  the  general  equations  of  hydro- 
dynamics (paragraph  6)  and  of  providing-  an  exact  classification  of  the 
statistical  parameters  of  turbulence  (paragraph  k) ,  which  is  interesting 
to  keep  in  mind  when  studying  the  foreign  reports,  too  exclusively  devoted 
to  spatial  correlations.   In  any  case,  the  reader  who  wants  to  read 
chapters  IV  and  V  can  pass  up  most  of  chapter  I,  except  perhaps  para- 
graphs 1  and  2,  without  major  trouble. 

Chapter  II  deals  with  the  kinematics  of  statistical  mediums  and, 
particularly,  isotropic  mediums.   It  seemed  practical  to  include  at  the 
same  time  the  velocity  correlations,  the  theory  of  which  has  been  given 
in  almost  final  form  by  Karman,  at  the  end  of  the  second  period,  and  of 
the  spectrum,  the  theory  of  which,  due  to  Taylor's  initiative,  has  only 
been  achieved  very  recently.   Only  paragraphs  10  and  11  refer  to  statis- 
tical functions,  and  their  detailed  knowledge  is  not  indispensable  for 
reading  the  rest  of  the  chapter.   The  results  and  the  formulas  of 
chapter  II  are  constantly  applied  in  the  subsequent  chapters,  but  it  is 
not  necessary  to  know  the  proofs  which  are,  in  most  cases,  a  simple  matter 
of  calculation.   The  most  important  of  these  formulas  are,  moreover,  com- 
piled in  a  special  section  following  chapter  V. 

Chapter  III  is  a  mathematical  study  of  the  application  of  Navier's 
equations  to  turbulent  motion.   Paragraphs  17,  18,  19,  and  22  are  funda- 
mental.  Their  main  purpose  is  to  recall  Karman 's  results  of  1938  with 
some  improvements  and  some  supplements  of  more  recent  date.   The  para- 
graphs 20  and  21  review  a  certain  number  of  physically  reasonable  solu- 
tions of  the  Karman-Howarth  fundamental  equation.   Some  of  these  assume 
particular  importance  in  chapter  V  but,  first,  it  seems  advisable  to  give 
an  impartial  view  of  the  whole  and  to  proceed  progressively  into  the 
domain  of  the  concrete.   The  paragraphs  23  and  2k   deal  with  an  equation 
by  Heisenberg  which  involves  time  correlations  and  from  which  probably 
not  all  possible  results  have  been  extracted.   It  is  not  indispensable 
to  have  knowledge  of  this  in  order  to  continue.   Paragraph  25  contains  a 


NACA  TM  1377 


mathematical  account  of  Heisehberg's  numerical  theory  of  the  spectrum, 
which  is  taken  up  again  in  chapter  V  from  a  more  physical  point  of  view 
and  whose  examination,  contrary  to  paragraphs  23  and  2k,   proved  useful 
before  attacking  chapter  V. 

Chapters  IV  and  V  deal  with  new  physical  theories  involving  simi- 
larity hypotheses  and  producing  numerical  laws.   Chapter  IV  reviews  the 
works  of  Kolmogoroff  and  Weizsacker,  chapter  V  those  of  Heisenberg, 
Batchelor,  and  Townsend  on  the  decay  of  turbulence  created  by  grids . 

Finally,  in  an  appendix,  the  theoretical  discussions  of  chapters  III, 
IV,  and  V  are  illustrated  by  some  correlation  curves  and  spectrum  curves 
measured  directly  in  the  wind  tunnel  by  A.  Favre,  in  the  laboratory  of 
the  mechanics  of  the  atmosphere  at  Marseille,  or  derived  from  experi- 
mental curves  by  elementary  transformations. 

An  exhaustive  study  of  modern  theories  of  turbulence  calls  for  some 
knowledge  of  the  calculation  of  tensors,  probabilities,  and  statistical 
analysis,  besides  the  classical  conceptions  of  differential  and  integral 
calculus . 

As  regards  the  tensors,  knowledge  of  the  definitions  and  fundamental 
operations  with  rectangular  cartesian  coordinates  is  sufficient.   There 
are  a  number  of  articles  on  this  subject,  but  they  generally  lean  toward 
the  tensor  analysis  with  curvilinear  coordinates  for  which  there  is  no 
need.   (it  should  be  noted  that  the  tensor  analysis  plays,  in  contrast, 
an  important  part  in  the  theory  of  the  boundary  layer  around  an  airfoil. ) 
Incidentally,  there  is  available  a  little  book  recently  published,  by 
Lichnerowicz,  entitled:   Elements  of  tensor  calculus  (collection 
Armand  Colin).   On  mathematical  statistics,  the  book  by  Darmois,  pub- 
lished by  Doin  (1928),  can  he  consulted.   For  the  elementary  theory  of 
random  functions,  consult  the  first  part  of  Bass'  report  (ref .  3)  and 
the  appendix  to  d'Angot's  "complements  of  mathematics"  (editions  of  the 
Revue  d'Optique),  edited  by  Blanc -Lapierre .  More  detailed  information 
on  statistical  functions  can  be  found  in  Levy's  book:   "Stochastic 
Processes  and  Brownian  Motion"  (19^-8). 

The  present  report  deals  only  with  general  theories  which  are  valid 
whatever  the  physical  or  geometric  causes  of  turbulence  may  be,  as  is 
shown  in  chapter  IV.   Only  in  chapter  V  the  assumptions  are  limited  and 
the  study  involves  a  problem  of  decay  of  turbulence  that  is  compatible 
with  the  turbulence  in  wind  tunnels.   These  theories  are  probably  appli- 
cable to  mediums  of  extremely  diverse  scales,  from  the  microturbulence 
to  the  terrestrial  atmosphere  (Dedebant  and  Wehrle),  to  stellar  atmos- 
pheres and  to  interstellar  matter  (Weizsacker).   However,  among  the 
hypotheses  there  is  always  that  of  the  incompressibility,  which  precludes 
the  application  to  sonic  or  supersonic  flows.   This  is  not  an  indispensi- 
ble  hypothesis,  but  it  simplifies  the  calculations  considerably,  and  the 
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consequences  are  easy  to  check.   So  the  incompressible  turbulent  motions 
are  practically  the  only  ones  studied  up  to  now. 

Among  the  other  hypotheses  worthy  of  discussion,  those  referring 
to  processes  of  energy  transport  (paragraph  2k   and  chapter  V)  assume  that 
all  turbulent  energy  comes  from  the  motion  of  the  whole,  dissociated  and 
broken  up  by  the  obstacles  of  periodic  structure.  But  nature  furnishes 
examples  of  different  turbulences  such  as  the  so-called  thermal  turbu- 
lence, for  example,  where  the  source  of  energy  is  not  the  fan  of  a  wind 
tunnel  but  the  solar  radiation  suitably  transformed  into  kinetic  energy. 
Kolmogorof f ' s  theory  of  local  isotropy  applies  probably  to  thermal  turbu- 
lence, but  the  forms  of  energy  transport  in  the  spectrum  must  be  differ- 
ent from  those  encountered  in  wind  tunnels .   The  same  applies  to  the 
astronomic  turbulent  mediums  alluded  to  previously.   So  the  foregoing 
remarks  limit  the  scope  of  the  recent  theories,  in  a  certain  measure. 

Omitted  entirely  was  the  problem  of  turbulent  boundary  layer  which, 
experimentally,  depends  on  the  same  technique,  but  has  been  approached  by 
different  mathematical  methods.   Furthermore,  it  is  a  complicated  problem 
where  the  turbulence  is  neither  homogeneous  nor  isotropic. 

In  what  measure  are  the  results,  suitably  demarcated,  defined?   The 
statistical  character  of  the  turbulent  velocity  seems  a  clear  and  well- 
established  notion  and  consequently  the  mechanics  of  turbulence  will  be 
a  statistical  mechanics.   As  far  as  the  kinematics  (chapter  II)  are  con- 
cerned, we  are  therefore  on  solid  ground.   The  dynamics  of  turbulence 
(chapter  III)  itself  is  likewise  well  established,  by  means  of  the 
hypothesis  of  the  validity  of  the  Navier  equations.   This  hypothesis, 
generally  adopted  because  it  is  convenient  and,  one  might  say,  necessary, 
has  however  at  times  raised  considerable  doubts.  However,  while  experi- 
mental verifications  do  not  contradict  it,  there  is  yet  no  occasion  to 
reject  it. 

But,  what  should  be  expected  from  experimental  verifications?  First 
of  all,  it  is  found  that  the  theories  are  still  rather  imperfect.   In 
fact,  the  theories  are  usually  limited,  acceptable  in  the  limiting  con- 
ditions which  are  difficult  to  attain  actually  (very  high  Reynolds  num- 
bers, for  example).  No  rigorous  verification  should  be  expected  since 
the  true  conditions  are  too  far  removed  from  the  theoretical  conditions. 

On  the  other  hand,  the  accuracy  of  measurement  is  low.   The  original 
reason  for  it  lies  in  the  very  nature  of  the  turbulent  phenomenon,  and 
its  irregular  and  badly  defined  character.   Furthermore,  the  anemometers 
are  coarse  instruments,  their  operation  not  sufficiently  known  in  the 
presence  of  turbulence,  and  their  interference  with  the  fluid  sometimes 
a  little  mysterious.  All  this  helps  to  lower  the  experimental  precision. 
To  measure  the  "length  of  dissipation,"  for  example,  two  ways  are  open: 
one  is  to  measure  the  correlations  of  the  velocity  at  two  infinitely  close 
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points  which  has  no  practical  sense  since  two  anemometers  cannot  be 
brought  together  indefinitely;  the  other  is  to  use  one  anemometer  for 
measuring  the  mean  square  value  of  the  derivative  of  velocity.   But  the 
operation  of  taking  a  derivative  of  a  function  as  complicated  as  the 
velocity  is  by  itself  inaccurate  and  introduces  serious  scattering  of 
the  points  of  the  derivative  curve.   So  in  both  cases,  the  accuracy  of 
the  results  is  extremely  limited. 

In  consequence,  the  experimental  verifications  can  be  applied  only 
to  the  orders  of  magnitude.   To  illustrate:   to  check  whether  a  parameter, 
according  to  theory,  is  constant,  one  constructs  the  representative  curve 
of  the  parameter  and,  if  this  curve  has  a  sufficiently  extended  maximum, 
one  estimates  that  the  experiment  closely  confirms  the  theory,  conceding 
that,  as  the  experimental  conditions  more  nearly  approach  those  stipu- 
lated by  theory,  the  maximum  flattens  out  more  and  more,  and  one  does 
not  appear  too  severe  in  the  examination  of  the  scatter  of  the  test 
points. 

On  this  assumption,  the  verifications  of  the  experimental  laws  of 
turbulence,  due  in  particular  to  Townsend,  are  encouraging.  Therefore, 
it  is  well  to  retain  the  hypotheses  of  chapters  III,  IV,  and  V,  although 
some  of  them  obviously  have  their  limitations,  and  we  ought  not  hope  for 
results  greater  than  actually  can  be  given.   Take  an  example  drawn  from 
the  theory  of  the  spectrum,  for  instance.  At  the  beginning  of  decay  of 
turbulence  in  a  tunnel  and  at  high  Reynolds  numbers,  a  certain  spectral 
function  F(k)   is  of  the  form  Ck^  for  the  small  values  of  wave  num- 
ber k.  .  It  then  passes  through  a  maximum,  then  becomes  proportional 

to  k  '  ,  and  finally  approaches  zero  as  k  >  °°,  maybe  as  k~'.   But 

the  regions  of  the  axis  of  k  in  which  the  fragments  of  the  laws  to  be 
enumerated  remain  acceptable  are  badly  defined  and  connected  by  zones 
of  which  the  structure  is  not  known.   So  the  future  task  of  the  theorists 
will  be  to  combine  these  partial  results  into  a  single  acceptable  law  at 
least  to  the  extent  that  it  does  not  become  fundamentally  incompatible 
with  the  nature  of  turbulence,  since  k  =  0  up  to  k  =  °°.   Only  then 
will  there  be  a  true  theory  of  turbulence . 
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CHAPTER  I 

GENERAL  ASPECTS  OF  TURBULENCE 
THE  STATISTICAL  METHOD 


1.  Definition  of  turbulence: 

Our  purpose  is  to  discuss  the  operation  which  consists  of  measuring, 
at  a  given  point,  the  velocity  of  a  flow  which,  to  avoid  every  difficulty, 
is  assumed  to  be  steady  and  uniform  at  the  usual  macroscopic  scale.   The 
measuring  instrument  is  an  "anemometer"  with  approximately  determined 
dimensions  and  time  constant,  and  of  sufficiently  ideal  nature  so  as  not 
to  disturb  the  flow  by  its  presence .   If  this  anemometer  is  large  enough, 
it  measures  the  "velocity  of  the  main  flow"  of  the  fluid.   If  it  is  very 
small,  it  can  be  imagined  that  it  operates  in  discontinuous  manner,  never 
undergoing  the  influence  of  more  than  one  molecule  at  a  time.  An  anemom- 
eter sensitive  to  the  individual  action  of  molecules  is,  of  course, 
unattainable,  but  the  idea  of  such  an  instrument  is  convenient  for  repre- 
senting the  extreme  limit  of  fineness  of  kinematic  measurements  in  a  fluid. 
Between  these  two  extremes,  the  indications  of  the  anemometer  depend  upon 
the  structure  of  the  fluid.   It  may  happen,  by  exception,  that,  when  its 
dimensions  are  progressively  reduced,  the  velocity  which  it  indicates 
remains  unchanged  up  to  the  moment  where  the  individual  influence  of  the 
molecules  starts  to  make  itself  felt  and  where  the  indications  lose  all 
statistical  significance.   The  flow  is  then  said  to  be  laminar. 

But,  in  general,  the  matters  are  otherwise.   We  start  with  a  first 
anemometer  which,  through  its  dimensions,  fixes  a  certain  scale  of  meas- 
urement.  This  anemometer  measures  the  mean  speed  of  the  molecules  in  a 
certain  volume  V.   Then  it  is  replaced  successively  by  smaller  anemometers 
in  such  a  way  that  volume  V  decreases  progressively.   It  happens  that, 
from  a  certain  value  V]_  of  V,  the  numerical  indication  supplied  by  the 
anemometer  changes.   If  V  is  decreased  continuously,  the  new  indication 
remains  stable  up  to  a  certain  value  V2,  then  changes  again  and  so  on. 
The  intervals  (V^,   V2) ,    (V2,   V*),  .  .  .  characterize  the  various  scales 

of  turbulence,  and  the  motion  of  the  fluid  is  said  to  be  turbulent.  More 
exactly,  they  are  the  conditions  necessary  for  a  fluid  to  be  turbulent, 
and  which  must  be  defined  and  perfected  to  make  them  sufficient  and 
practical. 

The  last  value  Vn  of  the  series  V]_,   V2,  ...  is  that  from  which 
onward,  the  notion  of  average  loses  its  significance,  the  number  of  mole- 
cules contained  in  the  volume  Vn  not  being  large  enough  any  longer  for 
statistical  purposes.   The  series  V^,   V2,  •  •  •  can  be  discrete  or 
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continuous.   If  it  is  discrete,  it  still  does  not  imply  that  the  critical 
values  V-]_,   Vg,  .  .  .  are  mathematically  defined.   In  the  vicinity  of  Vj_, 
a  rapid  variation  of  the  anemometer  indications  occurs,  which  subsequently 
become  stabilized  in  the  region  (V-,,   Vp),  and  so  forth.   If  this  stabili- 
zation is  not  very  clear,  the  turbulent  scales  are  said  to  succeed  one 
another  in  continuous  fashion. 

The  previous  discussion  ends  with  the  notion  of  the  turbulent  fluid. 
But  the  definition  of  the  laminar  motion  given  above  is  a  little  too 
restrictive  and  the  distinction  between  laminar  and  turbulent  still  not 
precise  enough,  as  proved  by  the  following  example: 

Consider  the  motion  of  air  produced  by  stationary  waves  in  a  sound 
tube.   The  motion  of  the  whole  reduces,  at  rest,  to  large  scale.   But, 
at  each  point  there  exists  a  speed  other  than  zero,  a  periodic  time  func- 
tion which,  at  a  given  instant,  varies  periodically  from  one  point  to 
another.   This  motion,  lying  between  the  system  at  rest  and  the  molecular 
agitation,  has  not  a  turbulent  character. 

Turbulence,  as  shown,  implies  first  the  notion  of  scale.   But  it 
should  be  added  that,  at  a  given  scale,  each  component  of  the  velocity 
at  a  point  is  a  function  of  time  presenting  a  character  of  periodicity 
without  fundamental  period.   This  is  not  a  periodic  function  but  a  sum 
of  harmonics,  the  frequencies  of  which  are  not  multiples  of  an  identical 
fundamental  frequency.   This  irregularity  of  the  turbulent  agitation  is 
essential  and  distinguishes  it  from  sound  agitation,  or  preturbulent 
vortical  motions,  like  the  cellular  vortices  of  Benard.   The  mathematical 
symbol  for  the  turbulent  velocity  is  not  the  ordinary  Fourier  series,  but 
Fourier's  integral.   It  will  be  discussed  later. 

This  concept  of  irregular  agitation  at  a  point  as  function  of  the 
time  is  not  itself  sufficient.   It  makes  it  possible  to  differentiate  the 
turbulent  agitation  from  the  periodic  sound  agitation  (musical  sound), 
but  not  from  the  noise,  which  is  an  agitation  without  definite  period. 
What  distinguishes  the  noise  from  turbulence  is  the  fact  that  it  is  prop- 
agated by  waves  that  exist  on  surfaces  of  equal  phase,  and  consequently 
have  a  regular  spatial  distribution,  notwithstanding  the  irregularity  in 

the  time  of  the  local  velocity. 

A  descriptive  and  purely  kinematic  distinction  is  involved  here. 
Its  cause  (compressibility)  is  not  discussed. 

The  difference  between  turbulence  and  sound  agitation  should  become 
plain  from  the  following  example:   In  a  turbulent  wind  tunnel,  we  select 
at  a  point  the  longitudinal  component  of  the  velocity  with  a  hot  wire  and 
send  the  electric  current  of  the  hot  wire  to  a  loud  speaker  on  the  outside 
of  the  tunnel.   The  atmosphere  becomes  the  source  of  an  irregular  agi- 
tation, which  propagates  by  waves  and  is  not  turbulence,  although,  at 
a  point,  the  internal  turbulent  motion  and  the  external  sound  motion 
have  some  important  kinematic  elements  in  common. 
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For  a  given  scale,  each  component  of  the  turbulent  velocity  is  an 
irregularly  periodic  function  of  both  space  and  time. 

It  seems  that  turbulence  is  well  defined  by  these  kinematic  condi- 
tions and  hence  is  distinguished  from  all  other  more  organized  fluid 
motions . 

2.  Average  values  -  statistics: 

Figure  1  represents  the  record  of  a  component  u  of  the  turbulent 
velocity  as  function  of  the  time.  (Record  of  the  velocity  of  turbulent 
agitation  in  a  20-cm  by  30-cm  tunnel:   airspeed,  20  m/sec;  intensity  of 

turbulence,  5-10  ;  time  of  recording  is  0.03  second.) 

The  most  natural  method  of  measuring  the  mean  velocity  on  the  graph 
consists  in  forming  the  ordinary  integral 


um  =  |  J      u(t)  dt  (2-1) 


extended  over  the  total  duration  T  of  recording.   This  method  is,  in 
general,  not  very  satisfactory  because  the  operation  lacks  precision 
when  the  curve  u(t)   is  complicated. 

A  more  precise  method  consists  in  dividing  the  graph  by  parallels 
to  the  axis  of  t,  suitably  close  together  in  the  ordinates  u-^,   U2, 

u^,  .  .  . ,  in  measuring  the  number  n^  of  points  where  the  line  of  the 
ordinate  —  (u^  +  u^+i)  meets  the  curve  u(t)  and  then  in  computing  the 
quantity 

"  =  Z  ui  y  (2-2) 


n  =  y       n^  is  the  total  number  of  points  met  by  all  parallels, 
i 

For  this  calculation,  a  profitable  first  stage  consists  in  first 
constructing  the  graph  giving  the  corresponding  statistical  frequency 

f ^  =  — i-  for  each  velocity  u^.   Crossing  the  limit  obviously  makes  it 

n 
possible  to  plot  a  curve  of  frequency  f  =  f(u)   (fig.  2)  such  that  the 
proportion  of  the  values  of  the  velocity  comprised  between  u  and  u  +  du 
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p  00  

is  equal  to  f (u)du,  the  integral   /   f (u)du,  which  replaces  y>   — , 

J0  Z n 

being  equal  to  unity.   The  final  expression  of  u  is  then 

/CO 
uf(u)du  (2-3) 


It  replaces  (2-2)  and  should  be  compared  with  (2-1). 

The  practical  operations  enabling  the  replacement  of  v^     by  u 
correspond  to  well-known  mathematical  operations.   The  mode  of  com- 
puting u  is  that  of  a  Lebesgue  integral,  and  um  is  an  integral  of  the 
classical  type  of  Riemann.   If  Ujjj  is  computable,  both  methods  yield 
the  same  result.   But  it  may  happen  that  Riemann 's  integral  does  not 
exist  because  the  function  u(t)   is  too  complicated  mathematically. 
However,  in  general,  Lebesgue 's  integral  exists  (if  u(t)  is  measurable, 
and  naturally  bounded).   This  mathematical  case  corresponds  to  the 
practical  case  where  the  curve  u(t)   is  too  complicated  for  an  accurate 
continuation  of  the  integration.   The  function  f(u),  continuous  and 
differentiable  in  the  current  cases,  is  the  medium  which,  determined  once 
for  all,  replaces  the  calculation  of  Lebesgue 's  integral  by  that  of  an 
ordinary  integral  by  means  of  the  plotting  of  curve  f(u). 

From  the  function  f (u),  other  averages  can  be  computed.   For 
example,  the  amount  of  differences  of  the  speed  with  respect  to  its 
mean  value  can  be  figured  by  computing  the  mean  value  of   (u  -  u)  . 
Rather  than  defining  this  average  by  Riemann1 s  integral 


i  f      f"u(t)  -  uj2dt  (2-4) 

T  Jo 

it  is  simpler  and  more  precise  to  use  the  formula 


_  2     /         2 
(u  -  u)  =   /   (u  -  u)  f(u)du  (2-5) 


which,  once  the  curve  f (u)   is  plotted,  calls  only  for  operations  of  a 
simple  character. 
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Thus,  it  is  seen  how  much  the  construction  of  the  curve  f (u) 
simplifies  the  numerical  calculations  of  turbulence.   Various  experi- 
mental techniques  make  for  direct  attainment  of  this  curve  without  passing 
through  the  numerical  analysis  of  a  velocity  record.   For  measuring  the 

2 
averages  such  as   (u  -  u)  ,  it  is  often  advisable  also  to  use  specialized 

equipment  without  first  plotting  the  curve  f(u). 

3.  Random  variables  and  the  laws  of  probability: 

The  theoretical  significance  of  the  function  f(u)   is  analyzed. 
It  groups  the  statistical  data  contained  in  the  initial  curve  u(t),  with 
this  exception  that  the  chronological  order  in  which  the  velocities  actu- 
ally follow  one  another  does  no  longer  appear.   This  limitation  is  quite 
natural,  though,  and  it  will  be  seen  later  that  this  order  reappears,  in 
a  certain  measure,  by  the  introduction  of  space  and  time  correlations. 

Obviously,  only  statistical  data  can  supply  stable  information  on 
turbulence.   When  the  same  record  of  the  velocity  is  begun  again  several 
times  while  taking  every  reasonable  precaution  so  that  the  conditions 
are  identical,  it  obviously  results  in  curves  u(t)   which  absolutely 
are  not  superposable .   The  function  u(t)  has  not,  therefore,  the  char- 
acter of  permanence  that  is  suitable  for  representing  the  laws  of  a 
physical  phenomenon.   But  this  character  is  relevant  to  the  function  f (u) 

or  to  averages  such  as  u,   (u  -  u)   whose  values  are  characteristic 
numbers  of  the  investigated  flow,  and  which  are  derived  by  simple  mathe- 
matical operations  from  the  function  f(u).   Hence,  we  direct  our  atten- 
tion to  this  function  which  can  be  regarded  as  representing  the  first 
law  of  turbulence . 

To  say  that  the  velocity  is  characterized  by  a  curve  of  statistical 
frequencies  is  to  say,  by  comparing  the  frequencies  with  probabilities 
and  f (u)  with  a  density  of  probability,  that  this  velocity  is  a  statis- 
tical quantity;  f (u)du  is  the  probability  that  the  chosen  component  of 
the  velocity  is  contained  between  u  and  u  +  du. 

A  priori,  such  a  law  of  probability  could  be  dependent  on  the  time. 
That  would  correspond  to  a  turbulent  flow  for  which  the  laws  varied  with 
respect  to  time.   There  is  no  contradiction  to  the  initial  assumptions 
of  permanence  here.   It  is  a  question  of  scale.   In  order  for  the  experi- 
mental operation  by  which  f (u)   is  defined  to  have  any  meaning,  it  is 
necessary  that  two  conditions  be  realized  simultaneously: 

(1)  The  number  of  oscillations  in  the  time  interval  T  involved 
must  be  great  enough  to  furnish  satisfactory  statistics. 

(2)  The  laws  of  turbulence  in  this  time  interval  T  must  be  practic- 
ally permanent. 


NACA  TM  1377  11 


If  it  is  not  so,  it  might  be  difficult  to  reconcile  the  statistical 
theory  with  the  experiment.  Fortunately,  those  "ergodic"  conditions  are 
practically  always  realized  in  the  usual  cases  and  are  therefore  taken 
for  granted  in  the  following: 

The  velocity  has  three  components  Uj_,   u2,  and  u*  which  are  treated 
as  three  random  variables,  components  of  a  random  vector. 

It  should  be  noted  here  that  the  velocity  is  perhaps  not  sufficient 
for  characterizing  the  turbulence.   It  might  appear  useful  to  introduce 
other  quantities,  such  as  the  pressure,  which  should  be  treated  as  a 
statistical  quantity.   But,  owing  to  the  equations  of  motion,  this  then 
will  be  a  function  defined  by  the  velocity  and  its  derivatives.  For  the 
present,  it  is  assumed  that  the  turbulent  motion  is  sufficiently  well 
defined  by  its  velocity  so  that  the  problem  narrows  down  to  the  laws  of 
probability  applied  to  the  velocity. 

The  simplest  of  these  laws,  that  which  immediately  generalizes  the 
experimental  function  f (u),  is  the  law  of  probability  of  the  system  of 
three  components  u^_,   u2,   u*  of  the  velocity.   This  law  may  vary  as 

function  of  the  time  t   (problem  of  spontaneous  decay  of  turbulence)  and 
of  the  space  (variation  of  turbulence  in  terms  of  the  distance  from  the 
walls).   It  is  therefore  a  function  of  t  and  the  ordinates  x^,   x2,   x^ 
of  the  point  of  measurement.   Its  density  is  denoted  by 


f(ul>  u-2,      uy,   xx,   x2,   x5,   t) 


or,  abbreviated,   f(u;  x,  t). 

The  quantity  f  du^  du2  du*  or,  abbreviated,   f  du,   represents  the 
probability  that,  at  the  point  x^,   x2,   x*   (or  x)  and  at  the  instant  t, 
the  three  velocity  components  are  comprised  between  u-,   and  U]_  +  dui, 
u2  and  u2  +  du2,   u*  and  u?  +  du?. 

But  a  single  law  of  probability  defined  in  terms  of  four  parameters 
x-,,   Xp,   x,,  and  t  is  not  adequate  for  characterizing  turbulence.   It 
is  necessary  to  introduce  the  more  profound  concept  of  random  function 
and  to  consider  the  turbulent  velocity  as  a  random  function  of  space  and 
time.   This  is  the  random  velocity  field. 

This  point  is  now  to  be  defined.  An  isolated  random  quantity  is 
defined  by  its  law  of  probability.   But,  to  define  a  system  of  coexistant 
random  quantities  requires  more  than  just  their  laws  of  individual  prob- 
abilities.  The  stochastic  dependencies  or  correlations  between  these 
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quantities  must  also  be  known.   Taking  a  family  of  probability  laws  does 
not  give  the  right  to  speak  of  the  system  of  corresponding  random  quan- 
tities without  completing  the  data  by  those  of  the  correlations.   Sup- 
posing that  the  family  in  question  depends  on  a  continuous  parameter  t. 
We  know  then  an  isolated  random  quantity  U(t)   for  each  value  of  t. 
This  immediately 'suggests  grouping  the  U(t)   corresponding  to  the  var- 
ious values  of  t  in  a  well-defined  system  of  statistical  quantities. 
This  calls  for  the  introduction  of  the  correlations  between  the  U(t]_), 
U(t2),  .  .  .  corresponding  to  an  arbitrary  system  E  of  the  values  t-j_, 
\>2)    ...  of  the  parameter.   To  proceed  thus,  means  to  define  a  random 
function.   Naturally,  this  also  holds  for  laws  of  probabilities  dependent 
on  several  parameters.   Thus,  when  a  turbulent  medium  is  represented  as 
a  velocity  field,  the  system  of  velocities  at  each  point  and  at  each 
instant  precisely  constitutes  a  system  of  coexistant  statistical  quantities, 
of  which  the  physical  interactions  characterizing  the  structure  of  turbu- 
lence have  the  correlations  for  mathematical  description. 

Among  the  systems  E,  the  simplest  are  the  denumerable  systems  and 
even  the  finite  systems  and,  among  the  latter,  the  simplest  one  which  is 
not  trivial  is  that  of  two  elements,  that  is,  of  two  points  of  space  and 
time . 

The  concept  of  random  function  thus  suggests  the  comparison  of  the 
velocity  vectors  at  two  different  points  of  their  field  of  definition, 
that  is,  for  two  different  positions  x  and  x'   (x  represents  the 
point  of  the  coordinates  x-,,   Xp,   x*)  and  for  two  instants  t  and  t'. 

It  concerns  a  statistical  comparison  which  makes  it  possible  to  define 
the  law  of  probability  of  the  two  systems  of  the  velocity  components  at 
points  x,   x'   and  instants  t  and  t'.   This  law  has  a  very  clear 
physical  meaning  and  is  easy  to  define  experimentally  or,  at  the  least, 
to  construct  the  surface  of  (statistical)  frequencies  corresponding  to 
one  velocity  component  at  point  x  and  a  second  component  at  point  x', 
the  measurements  being  spaced  at  a  chosen  time  interval  t. 

In  practice,  the  question  is  frequently  handled  from  a  less  general 
point  of  view.   One  is  not  concerned  with  the  laws  of  probability  them- 
selves as  such  but  only  with  their  most  simple  moments,  those  of  the 
second  order  which  are  associated,  as  will  be  shown,  with  certain 
"physical"  aspects  of  turbulence  and,  possibly,  with  certain  moments  of 
the  third  order  that  play  a  part  in  modern  theories.   The  moments  of  the 
second  order  constitute  the  correlation  tensor  of  the  law  of  probability 
of  the  velocity  field  at  two  points  and  at  two  instants.   They  are  meas- 
ured direct,  without  resorting  to  frequency  curves  or  frequency  surfaces 
or  to  velocity  recordings.  A  detailed  study  follows  later. 
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k.   The  concept  of  random  point  -  velocity  field  -  turbulent  diffusion: 

The  density  of  probability  f(u;  x,  t)   of  the  velocity  field  contains 
the  velocity  of  the  whole  fluid,  but  not  its  density.   With  p(x-]_,  X2,  x^;  t) 
or,  abbreviated,   p(x,  t)  denoting  the  quotient  of  the  density  at  point  x 
by  the  fluid  mass,   p  is  a  normalized  function,  as  a  density  of  probabil- 
ity, which  means  that 

n 

p  dx  =  1  (k-1) 

v 


where  dx  represents  the  element  of  volume  dx^  dxg  dx*  and  the  integral 
is  extended  to  the  volume  V  occupied  by  the  fluid. 

In  all  modern  studies  on  turbulence,  the  fluid  is  naturally  assumed 
incompressible,  so  that  p  is  a  constant,  equal  to  l/V  in  the  volume  V, 
and  zero  at  the  outside,   p  could  be  simply  replaced  by  a  constant;  but 
the  more  general  conclusions  to  be  arrived  at  ultimately  are  more  complete 
if  this  simplification  is  not  made.   On  the  other  hand,  it  is  interesting 
to  foresee,  at  a  certain  stage  of  the  theory,  the  day  when  It  will  be 
possible  to  study  the  turbulent  motions  in  conditions  where  the  compressi- 
bility is  no  longer  negligible.   For  these  reasons,   p  is  treated  here 
as  a  function  of  x  (and  even  of  t,  if  necessary). 

The  product  R(x,  u;  t)  =  pf  is  now  formed.   It  obviously  is 
normalized  with  respect  to  the  system  of  the  six  variables  x,  u 


/R  dx  du  =  1  (k-2) 


R  presents  thus  the  characters  of  a  density  of  probability  with  respect 
to  these  six  variables.   The  quantities  x-j_,  yi^,    x*  are  regarded  as  the 
coordinates  of  a  moving  point,   U]_,  U2,  u*  as  the  velocity  components  of 
this  point.   These  are  six  random  quantities  of  which  the  law  of  probability 
at  instant  t  is  known.   Thus  a  random  point  can  be  associated  with  the 
turbulent  fluid  in  correspondence  with  a  given  scale. 

This  point  is  now  to  be  discussed  as  was  the  velocity  field  in  the 
preceding  paragraph. 

The  position  and  the  velocity  of  this  point  are  random  functions  of 
the  time.   The  theory  of  random  functions  suggests  the  study  of  the  law 
of  probability  of  the  system  of  positions  and  velocities  of  this  point 
for  an  arbitrary  combination  of  instants.   The  initial  analysis  was  on 


14 


NACA  TM  1377 


the  law  of  probability  at  one  instant.   The  generalization  from  one  to 
two  instants  seems  to  us  adequate  for  forming  a  physical  theory  of  turbu- 
lence and,  in  particular,  for  considering  the  statistical  organization 
in  time  of  the  velocity,  and  what  may  be  called  the  "interactions  between 
turbulent  particles." 

The  values  of  the  positions  and  of  the  velocity  must  therefore  be 
associated  to  two  instants  t  and  t'.   This  association  is  a  stochastic 
relationship  defined  by  the  density  of  probability  of  the  following  system 
of  12  statistical  variables: 


*ij  X2^  x^ 


Ul>  u2'  U3 


position 


velocity 


at  instant  t 


xl  }    x2 


ui'>  V  >  v 


position 
velocity 


>  at  instant  t ' 


or,  abbreviated, 


G(x,  x',  u,  u';  t,  t1 ) 


and  it  is  assumed  that  it  defines  the  turbulent  motion, 
can  it  be  measured  ? 


By  what  stages 


According  to  the  theorem  of  compound  probabilities, 
is  the  product  of  two  factors: 


G  dx  dx'  du  du' 


(1)  The  probability  that  at  the  instants  t  and  t'   the  statistical 
image  point  of  the  fluid  might  have  positions  contained  within  the  inter- 
vals  (x,  x+  dx)   and   (x1,  x'  +  dx1). 

(2)  The  probability  (conditional)  that,  these  positions  being  fixed, 
the  velocities  are  contained  in  the  intervals  (u,  u  +  du)  and  (u1,  u'  +  du' ) 

This  last  probability  is  designated  by  H(u,  u';  x,  x',  t,  t')du  du1 . 
According  to  the  theorem  of  compound  probabilities,  the  probability  (l) 
is  the  product  of  the  probability  p(x,  t)dx  that  the  random  point  is 
found,  at  instant  t,  in  the  interval  x,  x  +  dx,  through  the  (conditional) 


probability  that,  the  position  at  instant  t 
at  instant  t1   is  found  in  the  interval  x' , 


being  chosen,  its  position 
x '  +  dx ' . 
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This  last  probability  is  designated  by  p(x';  x,  t,  t')dx'.  We  can 
write 

G(x,  x',  u,  u",  t,  t1)  =  H(u,  u';  x,  x',  t,  t')p(x';  x,  t,  t')p(x;  t) 

CM) 

What  is  the  significance  of  the  three  factors  of  which  ■  G  is  the 
product  1 

We  already  know  p,  which  represents  except  for  a  numerical  factor 
the  density  of  the  fluid. 

The  function  H  is  the  law  of  a  random  velocity  field  at  two  points 
and  at  two  instants.   It  is  natural  to  identify  it  with  the  law  of  the 
field  already  discussed  at  the  end  of  paragraph  3-   It  is  seen  that  it 
gives  no  complete  picture  of  turbulence.   It  does  not  explain  the  func- 
tion p. 

The  function  p  represents  the  turbulent  diffusion  at  a  chosen 
scale.   It  is  the  relative  density  at  point  x'   and  at  instant  t'   of 
fluid  elements  which  have  passed  neighboring  point  x  at  instant  t. 
It  should  be  pointed  out  that  the  diffused  portion  of  the  fluid  is  essen- 
tially compressible  since  the  density  p(x';  x,  t,  t')  decreases  in 
proportion  as  the  point  x'   is  removed  from  the  initial  point  x  where 
it  is  maximum.   The  spread  can  be  materialized  and  p  can  be  measured 
by  introducing  with  the  necessary  precautions  at  point  x  a  dye  that 
spreads  in  the  fluid.   But  it  is  a  rather  ticklish  matter  to  separate 
the  effects  of  the  various  turbulent  scales,  especially  of  the  molecular 
diffusion.   It  is  accomplished  by  adapting  the  particles  of  the  "dye"  to 
the  chosen  scale.   This  way  Kampe  de  Feriet  rendered  the  turbulent  dif- 
fusion of  the  tunnel  flow  visible  in  his  experiments  at  the  Institute  of 
Fluid  Mechanics,  at  Lille,  by  injecting  soap  bubbles  at  a  point  in  place 
of  dyed  particles.   These  soap  bubbles,  because  of  their  size,  were 
sensitive  to  the  turbulent  fluctuations  and  were  used  successfully  for 
measuring  a  density  of  turbulent  diffusion. 

5.  Equations  of  development  of  the  laws  of  probability: 

The  analysis  of  the  turbulent  velocity  field  made  it  possible  to 
represent  a  turbulent  fluid  by  two  random  functions  X(t)  and  U(t) 
playing  the  part  of  the  position  and  of  the  velocity  for  a  random  material 
point.   It  was  shown  how  this  concept  of  random  point  gives  a  very  com- 
plete picture  of  the  turbulence.   But  this  picture  is  qualitative  and 
must  be  made  more  quantitative . 
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At  the  beginning,  no  distinction  is  made  between  the  velocity- 
vector  X(t)  and  the  position  vector  U(t).   The  six  components  of  these 
two  vectors  are  considered  as  those  of  a  unique  vector  X(t)   i-n  a  six- 
dimensional  space,  by  putting  u^  =  x^,  Up=  Xr,  u,  =  Xg,  and  designating 
the  density  of  probability  of  this  vector  by  R(x;  t),  this  notation  being 
the  abbreviation  of  R(x^,  x2,  Xj,  x]^,  Xc,  xg;  t),  that  is,  of 
R(xx,  x2,  x^,  ux,  u2,  u5;  t). 

We  already  had  applied  (4-3),  the  theorem  of  compound  probabilities, 
to  the  law  of  probability  G  of  the  position  and  the  velocity  at  two 
instants  t,  t ' : 

G(x,  x*,  u,  u';  t,  t')  =  fp(x;  t)p(x';  x,  t,  t')|H(u,  u'j  x,  x1,  t,  t1) 

The  question  involved  essentially  the  separation  of  the  velocity, 
which  figures  in  the  factor  H,  from  the  position.   But  there  is  another 
way  of  applying  this  theorem.   It  consists  in  separating  the  two 
instants  t  and  t '   by  writing 

G(x,  x',  u,  u1;  t,  t')  =  R(x',  u';  t')K(x,  u;  x',  u1,  t,  t') 

(5-D 

hence  it  results,  according  to  the  theorm  of  the  total  probabilities, 
that 

R(x,  u;   t)  =   /r(x',  u';  t')K(x,  u;  x',  u',  t,  t')dx'  du' 

(5-2) 

K  is  a  conditional  density  of  probability,  that  of  the  "probability  of 
passage"  of  the  state  of  the  random  point  at  instant  t'   to  its  state  at 
instant  t  >  t ' . 

Abbreviated,  we  get 

R(x;  t)  =   /r(x";  t» )K(x;  x',  t,  t')dx'  (5-3) 


To  exploit  this  equation,  recourse  is  had  to  a  method  patterned  after 
the  concept  of  J.  Moyal  (ref.  33)  •   This  idea  was  to  consider  (5-3)  as  a 
linear  integral  transformation  for  passing  from  the  density  of  proba- 
bility R  at  instant  t'   to  the  same  density  at  instant  t.   The  density 
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of  probability  of  the  passage  K  is  the  kernel  of  the  transformation. 
With  K^t1   indicating  the  linear  operation  which  has  as  kernel  the 
function  K(x;  x1,  t,  t'),  one  may  write  symbolically: 

R(x;  t)  =  Ktt,    R(x';  f )  (5-*0 

This  transformation  has  special  properties  rather  difficult  to 
define.   Suffice  it  to  state  that  it  reduces  to  the  identical  transforma- 
tion for  t  =  t'.   It  is  assumed  that  it  has  an  inverse  K^  i  "■*-  =  K^-i^. 
without,  however,  prejudicing  the  relations  between  K^ i   and  K-^'f 

The  infinitesimal  transformation  applied  to  K^t i   is  now  examined. 
To  this  end,   R  is  assumed  differentiable  with  respect  to  t  and  dR/dt 
calculated 

dR(x,  t)  =  lim  R(x;  t  +  t)  -  R(x;  t) 

dt     x  _^  n  t 


=  lim  Kt+Tt'  ^^  R(X',  f) 

T  -^0         T 


lim  ^if  *t't  -  1  E(X)    t) 


If,  as  assumed,  the  function  R  is  differentiable  with  respect 
to  t,  the  limit  of  the  second  member  exists,  and  it  is  a  function  of  x,t 
independent  of  t1.   It  can  be  computed  by  giving  parameter  t'   (inde- 
pendent of  variable  t)  any  value  not  exceeding  t,  such  as  t'  =  t,  for 
example .   Hence 

f£.  lim  JW^ilR  (5.5) 


bt 


T-^0 


The  operator 


L  =  lim 


Kt+Tt  "  1 
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defines  the  infinitesimal  transformation  of  K,  and  R  confirms  the 
fundamental  functional  equation 

M  =  lr  (5-6) 

bt 


which,  theoretically,  enables  R(x;  t)  to  be  computed  at  instant  t 
when  R(x;  t1)   is  known  at  an  initial  instant  t1. 

The  foregoing  calculation  cannot  be  explained  in  a  simple  manner 
from  formula  (5-3)  because  the  limit  of  the  function  K(x,  x' ,  t,  t') 

does  not  exist  when  t'  >  t;  this  is  a  symbolical  "function  of  Dirac" 

which  expresses  the  identity  K^  =1  in  the  functional  formalism.   It 

is  preferable  to  pass,  as  Moyal  did,  from  the  densities  of  probability 
to  characteristic  functions.  Moyal 's  calculation  follows : 

With  the  function  K  is  associated  the  characteristic  function  of 
the  increment  X(t)  -  X(t'),  the  value  of  X(t')   once  fixed,  or  by 
definition 

*(<x;  x*,  t,  t')  =  LMx-x')  K(x.  x.^  tf   ti)dx      (5_7) 

One  assumes  likewise: 

cp(a,  x,  t)  =  /eiaxR(x;  t)dx  (5-8) 

the  characteristic  function  of  X_(t). 

According  to  the  theorem  of  total  probabilities,  one  has: 

cp(a,  x,  t)  =   /eiaX'o(a;  x',  t,  t')R(x';  t')dx'        (5-9) 

This  relation  replaces  (5-3)- 

The  two  members  of  (5-9)  are  differentiated  with  respect  to  t 

*e  =   iim    £jq*'  »(q;  *',  1 1  ii  *')  -  *(*>  *>  *,  t-) 

at       T_^rW  T  ' 


NACA  TM  1377  19 


Whereas  the  conditions  of  differentiability  under  the   /   sign 

[permutation  of  signs  lim  and  of  /  J  cannot  be  verified  on  (5-3),  they 

can  be  here,  in  all  current  cases.   Since  the  limit  must  not  depend  on  t1, 
let  t1  =  t,  so  that 

^E=  reiaX'R(x',  t)   lim   S(a?  x' ,  t  +  t,  t)  -  1  dx,    (5_1Q) 


because  $ (a;  x,  t,  t)  =1,  for  t'  =  t. 

The  limit  that  figures  under  the   /   sign  is  a  certain  function 
9 (a;  x1,  t'),  whence  follows  the  basic  equation 


^B  =  /eiax,R(x',  t)9(a;  x1,  t)dx' 
ot    J 


equivalent  to  (5-6).  From  that,  two  probability  densities  can  be  recovered 
by  taking  the  Fourier  transforms  of  the  two  members.   Lastly,   dR/dt  is 
expressed  in  form  of  an  integral  transformation  of  R,  equivalent  to  the 
transformation  L,  which,  after  introducing  a  kernel  function  L(x,  x',  t), 
gives 


dR(x'  t}  =  /r(x',  t)L(x,  x',  t)dx'  (5-11) 

St      J 


The  second  member  of  (5-H)  is  represented  in  integral  form  but,  in 
many  cases,  it  can  be  expressed  in  form  of  a  differential  operator  of 
finite  or  infinite  order. 

Examples.-  Supposing  the  probability  of  passage  obeys  the  Laplace- 
Gauss  law  and,  to  avoid  any  confusion  of  the  notations,   to  now  denotes 
the  previous  instant;  the  differences  of  X ftg )   and  X(t)  are  denoted 
by  Sq  =  S  (tQ J   and  S  =  S(t),  the  correlation  coefficient  between  Xrto) 
and  X(t)  by  r  =  r'tQ,  tj.   It  is  known  that  the  law  of  probability 
related  to  X(t),  when  the  value  Xq  of  Xrto)   is  given,  has  for  density 


1 


x  _  j.^ 


W       +   +\        1        2  1-r2  s    s0 
K(x;  x0,  t0,  tj  =-—-—=  e   V    JL 

/2nSVl  -  r2 


(5-12) 
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Using  Moyal's  method,  it  is  shown  that  the  probability  density  R(x,  t 
which  assumes  K  as  kernel  of  the  probability  of  passage  satisfies  the 
equation  of  partial  derivatives 


2R  =  .  fr.  +  Sl\l_(xR)  .  r-s  ¥*  (5-13) 


where  r'   represents  the  value  of  — —  r  (ti,  t)   for  t]_  =  t,  and   S' 

St^ 

denotes  the  derivative  —  S(t) . 

dt 

The  equation  (5-13)  is,  moreover,  demonstrated  very  simply  by  a 
direct  method.   By  definition  of  K 


1 


x  _  r;*o 


,12 


r(x,  t)  =  -      X       f    R(x0,  tQ)e  W-^)lS  soJ  ^      b_lk) 


r2 


this  equation  being,  in  particular,  satisfied  when 

R(x,  t)  =   1  e  2S 
S|/2jt 

When  the  two  members  of  (5-1*0  are  differentiated  with  respect  to 
t  and  x,  equation  (5-13)  is  verified.  Reciprocally,  the  integral  of 
(5-13)  which  is  reduced  to  R(x,  tQJ  for  t  =  tQ,  can  be  put  in  the 

form  (5-lM  ■ 

The  construction  of  the  random  functions  compatible  with  these  laws 
of  probability  is  an  easy  matter. 

To  illustrate: 

Let  h(s)   be  a  variable  random  function  of  the  parameter  s,  obeying 
the  reduced  Laplace-Gauss  law   fh  =0,  h  =  1],  and  so  that  the  increments 
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h(S]\  -  h  fs2)    and  n(s3)  -  n  (sl±)  ,    corresponding  to  two  separate  intervals 
(s±>    s2>)   an<^  (s~5>    stf\  >    are  independent.  The  random  function 


X(t)  =  f    (t  -  s)dh(s)  (5-15) 


is  now  considered. 

It  is  easily  shown2  that  X(t)   obeys  a  Laplace-Gauss  law  having 
for  typical  difference  S(t)  =  ^jL-     and  that,  if  T  =  t  -  tQ 

+  2: 


'('•^--^ds-1'!**---  (5"l6) 

(*  ♦  I)' 


Without  going  into  details,  it  is  simply  recalled  that  the  demonstra- 
tion utilizes  as  intermediary  the  characteristic  function  of  X(t),  of 
which  the  logarithm,  owing  to  the  properties  of  independence  of  the 
dh(s),  is  expressed  by  an  elementary  integral. 

From  these  formulas,  it  follows  that 


r-  .0,  §l-*i 

'  S    2  t 


and  that  R(x,  t)   verifies  the  first-order  partial  differential  equation 


M  +  I  I  JL(xR)  =  o  (5-17) 

dt   2  t  dx 


2See,  for  instance,  Bass,  reference  k,    or  "The  random  functions  and 
their  mechanical  interpretation."  Revue  Scientifique,  No.  32^0,  19^5 • 
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X(t)   is  a  differentiable  random  function.   Its  stochastic  deriva- 
tive is  deduced  from  the  expression  (5-15)  by  operations  of  classical 

rt  3  x 

form  and  written  as   /   dh(s).   It  is  easily  proved  that  is  the 

J  Q  2  t 

related  mean  of  this  derivative,  if  the  value  x  of  X  is  fixed.   The 
equation  (5-17)  has  many  solutions  which  are  densities  of  probability. 
It  does  not  determine  the  function  K. 

Next  consider  the  elementary  random  function 


X(t)  =  f    dh(s)  (5-18) 

Jo 


stochastic  derivative  of  the  function  which  has  been  studied  as  the  first 
illustrative  example.   It  is  shown  that  X(t)   obeys  a  Laplace-Gauss  law 

having  S(t)  =  /t"  for  typical  difference  and  that,  if  tn  <  t: 

r  ft»  to)  =  /5  =  1  ~   -  —  +  •  ■  •  >     r'  =  "  —       (5-19) 
\ '  0/  y  t      2  t0  2t 


Therefore,  R(x,  t)   verifies  the  second-order  partial  derivative, 
equation 


SR   1  ci2R  , 

*  =  I  S2  (5-20) 


This  is  the  equation  of  heat.   The  solution  which  reduces  to  a 
given  function  Rq(x),  for  t  =  t_,  is  given  by  the  classical  formula 

R(x,  t)  =  f    "    J__e"^:^'(X"Xb)2Rfa)dxo 


which  can  also  be  shown  by  using  the  Fourier  transform  of  R,  that  is, 
its  charactistic  function. 
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Thus  equation  (5-20)  defines  here  the  form  of  the  probability  of 
passage  k(x;  Xq,  t,  tQ j  contrary  to  what  occurred  in  equation  (5-17) • 

This  example  presents  an  unusual  peculiarity.   It  is  easily  verified 
that  the  function 

K(Xi  V  *'  V>  "  ^T%  -'*fef™)2  (5-2D 

satisfies  the  functional  equation 

K(x2;  xq,  t2,  t0)  =  /   K(x2;  xx,  t2,    t1)K(x1;  xq,  t±,    t0)dx1    (5-22) 


called  the  Chapman-Kolmogoroff  equation,  which  characterizes  the 
Markoff  processes  (or  more  generally,  the  "pseudo-markovian"  processes), 
a  functional  generalization  of  simple  Markoff  chains.  This  equation 
can  be  written  in  operational  notation  as 


Kt2t0  =  Kt2tL  Kt1t0  (5-23) 


It  expresses  that  the  operations  K  form  a  group.  The  operator 
H+Tt'  H't is  then  vritten  simply  as   t+Tt  — .   It  is  independ- 

T  T 

ent  of  t*.   In  the  general  case  it  is  its  limit  L  only  when  t— >0 
which  must  be  independent  of  t,  but  it  itself  is  not. 

Returning  to  equations  (5-12)  and  (5-13)  it  now  is  assumed  that 
R  is  not  dependent  on  t;   S  is  then  a  constant  and  S'  =  0.   If 
r1  ^  0,  equation  (5-13)  is  written  as 


A(xR)  +  S2^  =  0  (5-21+) 

Sx         Sx2 
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The  only  solution  of  this  differential  equation  which  defines  a 
law  of  probability  is  the  density  of  the  Laplace-Gauss  law 

_  jsi 

r(x)  =  _A_  e   2s2  (5_25) 

This  example  contains,  as  a  special  case,  stationary  random  functions 
for  which  the  function  Wt,  tQ )  depends  solely  on  the  difference 
t  -  tg;  r'   is  then  a  constant,  independent  of  t. 

The  final  example  deals  with  a  vectorial  random  function  of  a  type 
to  be  utilized  later.  Consider  simultaneously  the  random  function 


x(t)  =  f    (t  -  s)dh(s)  (5-26) 

J  0 


and  its  derivative 


X'(t)  =  f     dh(s)  C5-27) 

■J  0 


which  can  play  the  part  of  the  speed. 

It  is  easy  to  form  the  'functions  R  and  K  for  the  vector  X 
having  X  and  X'   as  components.  To  find  the  partial  differential 
equation  verified  by  the  function  R(x,  x',  t)   it  is  not  necessary  to 
first  form  the  kernel  K.   It  is  simpler  to  begin  with  the  expressions 
of  X  and  X',  which,  passing  through  the  intermediary  of  the  charac- 
teristic function  of  X,  X',  gives 


M  +  x-  M  =  A  2I£_  (5_28) 

dt      dx 
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The  operator  L  is  thus  a  differential  operator,  of  the  second 
order  with  respect  to  x' : 


L  =  -x-  A  +  i  -*L  (5-29) 

ax     2  ax*2 


Returning  now  to  the  old  notations  and  adding  the  letter  u  to 
the  velocity  vector,  we  get 


L  =  -ul+^  (5-30) 

dx   2  du^ 


considering  X'(t)   as  the  velocity  of  point  X(t)  .  The  operator  L, 
which  is  the  subject  of  this  example,  appertains  to  the  particular  class 
of  operators  which  make  it  possible  to  define  statistical  kinematics  and 
consequently  hydrodynamics,  as  will  be  proved. 

6.  Random  velocity.  Hydrodynamic  equations  (refs.  3,  ^) ' 

The  problem  involves  the  separation  of  what  is  position  and  what 
is  velocity  in  the  vectorial  random  function  with  six  dimensions  X(t). 
In  classical  mechanics  the  velocity  U(t)   is  the  derivative  of  the 
position  X(t)  .   In  the  present  case,  the  assumption  is  made  that  the 
velocity  is  the  stochastic  derivative  of  the  mean  square  of  the  position. 
By  theory  of  random  functions,  it  follows  that,  if  the  random  function 
X(t)   is  completely  or  at  least  known  locally,  the  function  U(t)   can 
be  defined  by  operations  comparable  to  those  used  in  classical  mechanics 
to  deduce  the  velocity  from  the  position. 

The  preceding  paragraph  contained  an  example  (equations  (26)  and  (27)) 
of  two  random  functions  X(t)   and  X'(t)  =  U(t)   linked  by  a  relation 
of  this  kind,  proving  in  a  general  way  that  there  is  an  operator  that 
satisfies  the  necessary  relation: 


ARdu  =  0  (6-1) 
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and  is  such  that  the  basic  equation  (5-6)  (with  original  notations) 
takes  the  form 

M  +  2~  uk£2_  =  AR  (6-2) 

St    k    6xk 


If  A  is  a  differential  operator,  it  can  be  expressed  by  an  expan- 
sion in  series  (infinite  or  finite)  in  terms  of  the  symbolical  powers 

of  the  partial  derivatives  of  which  the  coefficients  are  functions 

auk 

of  x  and  u.   The  vectorial  operator  having  for  components  occurs 

cbck 

only  in  the  first  member,  by  its  scalar  product  with  the  velocity. 

The  product  G  =  RK  can  then  also  be  decomposed  in  G  =  ppH  (com- 
pare formulas  (4-3)  and  (5-1)).  The  function  H  defines  the  correlations 
in  the  velocity  field  and  function  p  is  the  mathematical  representation 
of  turbulent  diffusion. 

Equation  (6-2)  must  therefore  play  the  fundamental  part  in  the  theory 
of  turbulence.   No  attempt  is  made  here  at  particularization;  it  simply 
is  shown  that  it  contains  the  general  equations  of  hydrodynamics. 

First,  the  two  members  of  (6-2)  are  integrated  with  respect  to  uk, 
with  due  regard  to  (6-1).  We  introduce  the  density  of  the  fluid 


p(x,  t)  =  /r(x,  u;  t)du  (6-3) 


(or  more  accurately,  a  quantity  p  normalized  and  numerically  propor- 
tional to  the  density)  and  the  relative  mean  of  the  velocity  for  a 
given  position,  the  components  of  which  are 


ui  =  -  /uiR(x>  u>  t)du  (6-M 


We  end  with  an  equation  of  continuity 


£♦££(*) -°  (6-5) 


NACA  TM  1377  27 

which  simply  expresses  the  fact  that  the  velocity  is  the  (stochastic) 
derivative  of  the  position.  Thus  it  is  seen  that  the  velocity  of  the 
whole  of  the  turbulent  motion,  defined  qualitatively  in  paragraph  h,   has 
the  first  of  the  qualitative  properties  of  the  hydrodynamic  velocity  in 
general.   At  the  chosen  turbulent  scale  (which  is  arbitrary)  the  mass 
is  conserved. 

Now  it  will  be  seen  that  the  velocity  also  satisfies  the  equations 
of  motion.  Multiplying  equation  (6-2)  by  the  velocity  component  u. 

and  integrating  with  respect  to  u,  gives  for  the  first  member 


^ +  ¥  4riUk 


Introducing  the  speed  of  fluctuation 


u'i  ■  H  -   "i 


and  considering  the  equation  of  continuity,  we  put 


Tik  =  "Pu'iu'k  (6"6) 


The  first  member  becomes 


^i   ^-  -H   1  ^-  ^ik 
at  k    dxR   p  k  oxR 


the  second  member  is  written 


u-j_AR  du  =  p7i 

7.   being  a  certain  function  of  x  which  is  deduced  from  the  operator  A 
and  the  probability  density  R.   Hence  the  equations  of  motion 
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du-j^       ,_ —    du.j_  ]_    ^- —    dTjk  ., 

dt  k        kdxk  i        p     k     dxjj 

The   7.   play  the  part  of  the  components  of  the  density  of  the  field 
of  external  forces  (gravity,  for  example).   The  Tj_^  are  the  components 
of  the  stress  tensor  (with  fixed  turbulent  scale )3.   They  are,  except  for 
the  factor  -p,  the  components  of  the  correlation  tensors  of  the  velocity 
components  at  a  point.   They  can  be  measured  directly  by  statistical 
methods,  and  the  equations  of  hydrodynamics  for  a  given  turbulent  scale 
can  thus  be  verified. 

It  should  be  noted  that  in  this  case  the  equations  of  motion  of  the 
ensemble  (average)  are  involved.   For  the  present,  nothing  about  the 
behavior  of  the  rate  of  fluctuation  has  been  assumed.   Later  on  the  usual 
assumption  will  be  made  that  the  velocity  satisfies  the  hydrodynamic 
equations  and,  more  precisely,  the  Navier  equations.   It  should  be  remem- 
bered that,  based  upon  this  hypothesis,  Reynolds  was  able  to  establish 
equations  similar  to  equations  (6-7)  for  the  mean  turbulent  motion.   It 
is  apparent  that,  without  it  being  necessary  to  repeat  Reynolds'  calcu- 
lations, the  statistical  theory  in  question  here  is  entirely  different 
from  that  of  Reynolds.   It  is  probably  more  complete,  but  it  still  has 
not  been  pushed  far  enough  to  be  verified  by  experiment,  due  to  a  lack 
of  suitable  hypothesis. 

As  simple  example  of  the  fundamental  equation  (6-2),  the  case  is 
chosen  in  which  the  function  U(t)  is  itself  differentiable  (in  mean 
squares,  ref s .  15,  16).  In  this  case,  the  statistical  image  point  of 
the  fluid  has  an  instantaneous  acceleration.   If  rv(x,  u,  t)   represents 

the  relevant  mean  of  the  acceleration,  that  is,  of  the  velocity  derivative, 
when  the  position  and  the  velocity  are  fixed,  it  proves  that  the  density 
of  probability  of  X(t)   and  U(t)   satisfies  the  equation 

dR   <?-        dR    v2-   d 


dt   k=l   oxk   k=l  ouk 

The  operator  A  =  -  5_  - — (rkR)   ^s  therefore  linear,  differential 

k=l  duk 

and  of  the  first  order.   The  field  of  the  external  forces  to  which  the 
fluid  is  subjected  has,  necessarily,  for  components,  the  quantities 


^This  expression  of  stresses  was  originally  given  by  Reynolds  (On 
the  dynamical  of  incompressible  viscous  fluids  and  the  determination  of 
the  criterion.   Phil,  trans.  Roy.  Soc .  CLXXXVI,  part  I,  123,  1895), 
proceeding  from  the  Navier  equations.   However,  the  exact  meaning  of  the 
Reynolds  stresses  is  different  from  that  of  Tj^. 
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*  -  -i/-i  x:  4Ma-  ■  lJv *■ -  h         (6-9) 

averages  formed  of  the  acceleration  components  when  the  position  is  fixed, 
provided  only  that  Uj_r/R  tends  toward  zero  when  the  velocity  increases 

to  infinity.   If,  in  particular,  the   Fj_   are  not  dependent  on  u,   T^   is 
identical  with  y^,   and  R  satisfies  the  simple  equation 

t?  +  H  "Jr- +  YL  yj*r-  =  °  (6-10) 


Once  the  7^  are  given,  this  equation  defines  R(x,  u,  t)   from 
R(x,  u,  tQ  \.      But  the  form  of  the  linear  transformation  from  R(x,  u,  tgl 
to  R(x,  u,  t),  that  is,  the  probability  of  transition,  depends  upon  the 
form  of  the  original  probability  R'x,  u,  tg;,  contrary  to  what  happens, 
say  for  the  equation  of  heat  (type  (5-19));  this  illustration  does 
not  appear  to  rest  on  hypotheses  sufficiently  inspired  by  reality  to 
serve  as  basis  of  a  turbulence  theory.   First  of  all,  the  starting  point 
must  be  modified,  as  will  be  done  in  the  following  paragraph: 

7-  Systems  of  molecules  (refs.  13  and  Ik): 

The  statistical  quantities  to  which  the  analysis  of  velocity  records 
leads,  represent  only  certain  scales  of  turbulence,  those  which  correspond 
to  the  ensemble  of  "vortices"  whose  dimensions  are  superior  to  a  limit 
approximately  fixed  by  the  employed  anemometer .   Can  the  theory  be  changed 
so  that  all  the  possible  scales  can  be  represented  simultaneously?  It 
seems  that  it  suffices  for  this  purpose  to  start  from  the  finest  scale, 
that  is,  the  molecular  scale.   The  gas  is  therefore  considered  as  a  system 
of  N  molecules,  and,  to  explain  the  method  with  as  much  simplicity  as 
possible,  the  assumption,  which  is  not  verified  for  air,  is  made,  that 
the  molecules  are  identical,  monatomic,  comparable  to  material  points 
subjected  to  central  interactions.   V(r)   denotes  the  potential  of  the 
force  of  interaction  of  two  molecules  separated  by  the  distance  r.   If 
r  is  great  (with  respect  to  the  diameter  of  the  molecules,  which  will 
not  be  introduced  explicitly),   V(r)   is  negligible.   For  the  small  values 
of  r,   V(r)   expresses  the  repulsion  of  the  molecules,  generalized  form 
of  shocks.   No  other  information  about  V(r)   is  needed  beforehand,  at 
least  in  a  general  theory. 

The  motion  of  this  system  of  molecules  is  controlled  by  the  equations 
of  dynamics .   But  the  extreme  complication  of  the  trajectories  of  the 
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molecules  prompted  the  replacement  of  rational  mechanics  by  statistical 
mechanics.  We  prefer  to  introduce  random  mechanics  where  each  molecule 
is  a  random  point.   By  means  of  an  ergodic  hypothesis  the  random  motion 
of  a  molecule  can  be  considered  as  being  statistically  equivalent  to  the 
motion  of  the  ensemble  of  the  fluid,  in  quasi-steady  conditions,  the 
ensemble  of  successive  states  of  the  visualized  molecule  replacing  the 
ensemble  of  the  simultaneous  states  of  all  the  molecules .  The  molecule 
would  therefore  be  the  concrete  image  of  the  abstract  random  point  serving 
up  to  now  for  representing  the  fluid.   The  molecular  scale  is  a  true  ulti- 
mate scale  of  turbulence,  separated,  however,  from  the  actual  turbulent 
scales  by  a  poorly  defined  but  finite  interval.  While  for  experimental 
turbulence  the  concept  of  a  random  image  point  is  a  mathematical  abstrac- 
tion, it  is  a  natural  idea  and  a  starting  point  for  the  molecules. 

This  idea  is  now  explored  but  by  a  method  slightly  different  from 
that  discussed  in  the  preceding  paragraphs.   The  probability  density 
%(xl>  •  ■  •>   XN'  ut  >    •    '    •>   Uy  >  ^)      °f  the  positions  and  the  velocities 
of  N  molecules  simultaneously,  rather  than  singly,  is  introduced.   The 
gas  appears  then  as  a  random  point  with  three  N  dimensions,  in  a  space 
of  configuration,  and  no  longer  as  a  random  point  of  ordinary  and  physi- 


cal space.   For  the  time  being,  the  notations 


xl >   x2> 


^N 


shall 


have  a  vectorial  character  and  represent  the  system  of  the  three  coordi- 
nates of  the  molecules^"  numbered  1,  2,  .  .  .,  N. 

The  potential  of  interaction  of  the  molecules  of  rank  i  and  j 
is  indicated  by  Vij,  and  the  external  force  to  which  the  molecule  of 
rank  i  is  subjected,  by  y^.  If  m  is  the  mass  of  a  molecule,  the 
equations  of  motion  of  the  molecules  read 


dxi 

dt 


=  u. 


du. 


m- 


dt 


7±  + 


(7-1) 


the  summation  applying  to  all  values  of  j   from  1  to  N,  when  V. .   is 
assumed  to  be  zero.   Naturally  V^ i  =  vii* 

When  these  equations  are  written,  the  hypothesis  is  made  that  the 
velocity  of  each  molecule  is  dif f erentiable .   It  follows ,  that  the  random 


f(x,  t)   represented 

of  components  x-j_,  Xn,   Xx- 


In  the  preceding  paragraphs  the  notation 
already  the  probability  density  of  a  vector  x, 

What  is  used  only  temporarily  is  the  meaning  of  the  subscripts,  particular 
to  this  paragraph.  A  change  in  notation  could  be  avoided  only  by  compli- 
cations in  writing  which  would  be  more  harmful  than  useful. 
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point  of  3N  dimensions  having  for  coordinates  the  ensemble  of  the  coordi- 
nates of  the  molecules  is  represented  by  a  random  vector  function  of  the 
configuration  space  of  3N  dimensions  doubly  differentiable.   This  prop- 
erty which,  as  stated  before,  is  not  necessarily  true  for  the  three- 
dimensional  random  point  used  previously,  is  the  reason,  according  to  the 
theory  of  statistical  functions,  why  fjj  verifies  the  partial  derivative 


equation 


*N   f   ^H   if      ^  SVijVfN   0        in  * 
+  >   u. +  -  }_    7i  +  /   > — -  =  °        (7-2) 


which  expresses  that  fjj  is  an  integral  of  the  equations  of  motion.   It 

has  the  same  form  as  Liouville's  equation  of  statistical  mechanics.   But 
it  should  be  remembered  that  its  original  meaning  is  a  little  different. 
This  equation  replaces  and  defines,  in  the  space  of  3N  dimensions,  the 
equation  (6-2),  and  must  now  be  exploited. 

The  turbulence  involves  "particles"  or  "eddies"  formed,  at  a  given 
scale,  by  groups  of  s  molecules,   s  being  a  very  large  number,  but 
at  the  same  time  very  small  with  respect  to  N.   Let  us  investigate  what 
functional  equation  is  satisfied  by  the  probability  density  fs  relative 

to  the  molecules  of  rank  1,  2,  .  .  .,  s.   fs  is  defined  by 


1 


fs    -    /%   *<8+l   dVl    •     •    •    *<N   dUN  (7-3) 


The  fundamental  equation  (7-2)  must  be  integrated  with  respect  to 
xs+l>   us+l»  ■  •  •*  XN>   UN-  ^e   °Pera"ti°n  1S  obvious,  except  for  the 

Svij 
terms  in  -;      V^  .   being  function  solely  of  il.   and  x . ,  the  integra- 
te.      J  J 
J 

tion  gives  a  zero  result  if  i  and  j  are  both  superior  to  s . 

1  ^vii  dfs 
If  i  and  j  are  both  inferior  to  s,  then . 

m  Sxj   di^ 

If  J  <  s,   i  >  s,  the  integration  with  respect  to  u^  shows  that 
the  result  still  is  zero. 

If,  finally,   i  S  s,   j  >  s,  we  first  can  integrate  with  respect 
to  x^,   u^,  for  k  ^  i. 
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1  avii  ofs+1 

The  result  is —  ,  where  fs+]_  represents  the  probability 

m  ox  .   ou . 
J    i 

density  relative  to  s  +  1  molecules  of  rank  1,  2,    .    .    .,      s,      j.   The 
last  integration 


-  /  t- — i dx<   du^   =  -  - —  /  - — if     ■,    dx,   du; 

mj    axj        Bui  J        J        m  oi^  J    Sxj      s+1        J        J 

Svij 

cannot  be  extended  farther,  because  both  ^ and  f-.i   are  functions 

0x4         t)"t~x 

of   Xj. 

Now  it  will  be  noted  that  in  the  sum 


I  ±  ±, 


>oV, 
- —  /  — — fg.-i  dx,  du, 
du,  J  ox.   S+1   J    J 


j=s+l  oui 


the  terms  are  identical  because  the  molecule  j  plays  an  anonymous 
part. 

The  final  result  is 

Sfs       v^        Sfs        If-        Sfs       N   -   s   ^      o      ^Vls+1^ 

T—  +  2_   ui^~  +  -  Z_    7i—-  =  Z_    —  / f  s+i   <3xs+i  dus+1 

at    i=i   dx^       m  i=i   du1     m   1=1  du^_J    °xs+i 

(7-M 


This  equation  seems  capable  of  serving  as  basis  for  a  theory  of 

turbulence.   The  scale  there  appears  explicitly  for  the  numbers  —     and 

s 

—  and  the  problem  (not  taken  up  here)  consists  in  formulating  a  reasonable 

hypothesis  which  enables  ?s+±     "to  be  expressed  with  the  aid  of  f„,  so 
that  (j-h)   becomes  a  functional  equation  in  f  . 


NACA  TM  1377  33 


For  s  =  1,   (7-M  is,  in  a  certain  measure,  comparable  to  (6-2) 
and  reads 


Sfx     Bfx   1  ^   N  -  1  d   /  ^V12 


+  Ul_i  +  -7l_i  =  /  _^f  dx  du      (7-5) 

ot      ox-j_   m  ou-j_     m   ou-j_^  0X0 

Equation  (7-5)  differs  from  (6-2)  by  the  presence  of  the  terms  in 
fp.  The  presence  of  the  second  member  relates  (7-5)  to  (6-2),  but  its 

form  makes  it  different  from  it.   Equation  (7-5)  reduces  to  a  functional 
equation  in  f^  only  if  fo  is  tied  to  f^  by  a  suitable  assumption. 

Such  an  assumption  has  been  made  by  J.  Yvon  who  demonstrated  with  its 
aid  Boltzmann's  fundamental  equation  on  which  the  kinetic  theory  of  gases 
is  based.  More  recently  it  was  taken  up  again  by  Born  and  Green,  who 
also  studied  the  case  of  s  =  2,  and  applied  it  to  the  case  of  liquids  on 
the  basis  of  an  assumption  by  Kirkwood  binding  f,  to  f  . 

It  seems  that  the  case  of  large  values  of  s  has  never  been  studied. 

Equation  (7-5)  like  (6-2)  and  despite  the  not  necessarily  linear 
character  of  the  second  member  with  respect  to  fs,  involves  equations 

of  hydrodynamics  which  are  easy  to  write,  if  the  method  indicated  in 
paragraph  6  is  applied.   Since  these  equations  have  so  far  not  been  used 
in  concrete  applications  to  turbulence,  only  these  summary  indications 
are  given  here.   The  application  of  equation  (7-M  to  turbulence  raises 
difficulties  which  are  far  from  being  solved,  or  even  stated  but  it  is 
important  to  know  that  this  equation  exists. 


CHAPTER  II 


CORRELATIONS  AND  SPECTRAL  FUNCTIONS 


8.  Introduction  -  correlations  in  space  -  homogeneity,  isotropy: 

In  this  chapter  we  deal  no  longer  with  the  general  laws  of  proba- 
bility of  turbulence .  We  limit  ourselves  to  the  study  of  the  correla- 
tions of  the  law  of  the  field  H(u,u'j  x,x',t,t').  Following  are  some 
preliminary  remarks  on  the  subject: 

Experiments  furnish  the  correlations  in  form  of  temporal  averages . 
Chapter  I  shows  how  it  was  possible  to  change  their  Interpretation  for 
converting  them  into  stochastic  averages .   The  assumption  is  made  that 
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this  operation  is  always  possible.   But  for  starting,  there  is  no  neces- 
sity for  knowing  whether  it  has  been  effected,  since  only  properties  of 
symmetry,  of  tensional  character,  are  involved  here.   On  the  other  hand, 
arriving  at  the  dynamics  of  turbulence  we  shall  see  that  the  use  of  tem- 
poral averages  leads  to  serious  difficulties,  which  disappear  when  they 
are  transformed  into  stochastic  averages,  and  our  calculations  will  deal 
only  with  stochastic  averages .   Temporal  averages  are  resorted  to  only 
when  the  calculations  have  reached  the  laws  which  must  be  compared  with 
physical  reality. 

The  components  of  the  speed  f luctuations-^  at  point  X]_,X2,x*  are 
designated  by  u^_,   Ug,  and  u*.   They  are  statistical  quantities  whose 

mean  value  is  zero.   Time  plays  no  part  at  present  and  the  parameter  t 
in  the  formulas  is  disregarded.   The  correlations  of  the  velocity  between 
the  two  points  x  and  x '  =  x  +  (■   are  defined  by  a  tensor  having  for 
components 


Rap(x,x')  =  ua(x)up(x')  (8-1) 


The  scalar  is  introduced  also: 


=  ^IW  =  zZua(x)ua<x')  (8-2) 


a 

The  turbulence  is  said  to  be  homogeneous  (within  a  certain  field  of 
space)  when  the  R  □  are  not  separately  dependent  on  the  two  points  x,x' 

but  only  on  their  relative  position,  that  is,  on  the  vector  |  =  x'  -  x. 
Then  Rap(|)   replaces  Rap(x,x'). 

The  turbulence  is  said  to  be  isotropic  at  point  x  if  it  is  homo- 
geneous in  the  vicinity  of  x  and  the  tensor  R  „   is  invariant  to  any 

rotation  of  the  axes  and  any  reflection.   The  form  which  isotropy  imposes 
on  the  components  of  this  tensor  is  discussed  later. 

These  two  definitions  concern  only  the  second-rank  tensor  R<xp. 
They  are  "second-rank  properties"  of  the  statistical  vector  Uq,.   They 
must  be  extended  to  certain  third-rank  tensors .   Homogeneity  and  isotropy 
can  be  given  a  more  complete  and  also  more  restrictive  definition  by 
extension  to  all  possible  moments  of  velocity  components  taken  at  any 
number  of  points.   It  is  then  more  convenient  to  define  the  homogeneity 

-^They  were  called  u'-^u'gjU' *   in  chapter  I.   But  it  is  advantageous 
to  use  thereafter  the  notation  u'   for  other  purposes. 
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as  the  invariance  to  translations  of  the  law  of  probability  of  the 
ensemble  of  velocity  vectors  having  for  origin  a  certain  number  of 
arbitrarily  distributed  points,  and  the  isotropy,  once  homogeneity  is 
achieved,  as  invariance  to  rotations  and  reflections. 

Without  restrictive  hypotheses  the  tensor  R  n  has  nine  distinct 

components,  functions  of  the  seven  variables  X]_,X2,x^,x'-l,x'2,x'  *,t . 

If  the  turbulence  is  homogeneous,  the  nine  functions  remain,  but  of  four 
variables  only  |-.,£p>£vt.   The  isotropy  is  now  written. 

To  the  tensor  RaR(i)  we  associate  the  scalar  bilinear  form 


where  Xa  and  Y   are  two  arbitrary  vectors .   This  form  is  the  mean 


value  of  the  product  of  two  scalar  products  2. uo(x)Xcl  an(^ 


ua(x')Ya.   If  it  is  invariant  to  rotations,  it  is  an  algebraic  com- 
bimation,  separately  linear  with  respect  to  Xq,  and  Ya,  of  invariants  of 
vectors  Xa,Ya,|a  with  respect  to  rotations.  These  invariants  are 


a 


Y~  X  Y       V5X       Y~  i   Y  r2=Y~i 

^ —  a  a      /        a  a      Z. a  a  Z —  c 

There  are  therefore  two  scalars  A(r),  B(r),  functions  of  r,  so  that 


YL   Ra3XaYp  =  A  51  XaYa  +  B   XI  ^aXa  )   H  hYfi)  (8-3) 


The  identification  shows  that,  if  8ag  represents  the  classical 
symbol  of  Kronecker,  zero  if  p  ^  a,  equal  to  unity  when  p  =  a: 

Rap  =  A5ap  +  B^p  (8"i+) 
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The  tensor  RaR,  in  the  isotropy  hypothesis,  depends  therefore  solely 

on  two  distinct  functions  of  two  variables  r  and  t.   The  notations  are 
changed  as  usual  and  one  puts  with  Karman  and  Howarth  (reference  30): 


RaP  "  U02 


■  \i0   +  g5 


L 


^2^abp     a(3 


(8-5) 


f  and  g  being  functions  of  r,   t  becoming  unity  for  r  =  0,  and 
Uq^  a  simple  function  of  t. 

If  r  =  0,  R  R  is  reduced  to 


up 


V°)  =  U02^°)5ap  (8"6) 


It  is  seen  that 


U]_  =  ^2  =  ^  =  V      U2U3  =  U3U1  =  \U2   =  °      ^8"7^ 


Also,  owing  to  the  isotropy,  the  tensor  R  R  is  symmetrical. 


Lastly,  it  should  be  noted  that  the  quadratic  form    _  R  „X  X„   is 

<xp  a,  p 

the  first  member  of  the  equation  of  an  ellipsoid  of  revolution.   This 

shows  that,  if  one  takes,  no  matter  how,  a  symmetrical  table  of  numbers 

R  R,  they  are  not,  in  general,  components  of  a  correlation  tensor.   These 

numbers  must  verify  the  inequalities  which  state  that  they  are  coefficients 
of  the  first  member  of  the  equation  of  an  ellipsoid.   In  other  words,  the 
roots  of  the  equation  of  the  third  degree   |RaR  -  S6aR|  =  0   (equation  in 

S  classical)  must  be  positive. 

A  similar  method  permits  the  reducing  of  the  components  of  the  tensor 


Tap7  "  ua(x)uR(x)u7(x  +  i)  (8-8) 
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which  will  be  needed  later.   This  tensor  is  symmetrical  with  respect  to 
a  and  (3,  and  the  trilinear  form 


TaP7XaYpZ7 

is  invariant  to  rotations.   By  applying  the  classical  notations  of  Karman 
and  Howarth,  it  is  found,  after  a  few  calculations,  that 


T<xB7  ~  "0 


a>c             a>    c             h      <-             c-b-   2ar    ,    „ 
^a6B7  +  ^357a  +  p  75ap  +  rf ^^7 


(8-9) 


a,  b,  and  c  being  three  functions  of  r  and  t. 

9.  Properties  of  the  functions  f,   g,  a,  b,   c.   Incompressibility  : 

Other  forms  of  symmetry  less  particular  than  isotropy  can  be  vis- 
ualized, as  for  instance,  axial  symmetry,  or  invariance  to  rotations  about 
a  given  axis,  instead  of  about  a  point.   The  case  of  isotropy  is  that  in 
which  the  Ra|3  depend  on  the  smallest  number  of  separate  functions. 

The  functions  f  and  g  are  correlation  coefficients.   For  example, 
f  can  be  defined  by  taking  two  velocity  components  along  the  axis  of  x-, 

(direction  of  velocity  of  the  ensemble)  at  two  points  situated  on  a 
parallel  to  this  axis  (fig.  3) 


(r)  .  ul(xl'x2'x?)ul(xl  +  r."2.*?)  (91) 


"o2 


g  is  defined  by  taking  two  components  still  parallel  to  the  axis  of 
the  x-,,  but  at  two  points  located  on  a  line  perpendicular  to  this  axis. 


^The  authors  (reference  30)  use  q,  I,      k  where  we  use  a,  b,   c. 
But  the  letter  k  is  also  used  in  a  just  as  classical  although  more 
recent  fashion  to  designate  the  spectral  frequency  and  will  be  used 
with  this  meaning  in  the  present  report. 
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For  example 
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,(r)  =  Ul(Xl-X2^X5)Ul(X1^2  +  r>X^) 

%2 


(9-2) 


Experience  indicates  that  near  r  =  0,  the  functions  f  and  g  are 
continuous  and  twice  diffe rent i able  and  allow  a  tangent  to  the  "horizontal" 
origin  (fig.  k) . 

Theory  confirms  experience,  which  unfortunately  lacks  precision  when 
the  distance  r  becomes  small.  Much  greater  accuracy  is  obtained  by 
correlation  measurements  with  difference  in  time  (experiments  by  Favre, 
Report  to  the  VII"^n  Intern.  Cong,  of  Appl.  Mech.,  London,  19^8),  because 
only  one  anemometer  is  used  and  the  timing  can  be  reduced  as  much  as 
desired. 

Hence,  one  may  write  developments  of  the  form 


f(r)  =  1  +  ^£f"(0)  + 


;(r)  =  1  +  E-g"(0)  + 


(9-3) 


f"(0)  and  g"(0)  are  negative  quantities,  possibly  time  functions, 
having  the  dimensions  of  the  inverse  of  the  square  of  a  length.  The 
length 


A  = 


g"(0) 


(9-h) 


is  called  length  of  dissipation. 

The  triple -correlation  functions  a,b,c  have  interpretations  simi- 
lar to  those  of  f  and  g .   For  example : 


c(r)  =  ^-u12(x1,x2,x5W1Jx1  +  r,x2,x5) 
uq3 


(9-5) 
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It  follows  from  formula  (8-9)  that,  if  Tao   has  a  well  defined 

limit  when  | — >0,  this  limit  can  only  be  zero.   Hence,  a(0) ,b(0) ,c(0) 
are  zero.   In  a  general  way,   Tagy  is  an  odd  function  of  Za>ln>t.y 

Consequently,  considered  as  function  of  r,   c(r)  is  an  odd  function. 

o 
Its  development  has  no  term  in  r  .   Lastly,  it  is  shown  that  it  has  no 

term  in  r,  either.   The  coefficient  of  this  term  is  the  mean  of 

2aul   1  5   5   1  „    ul5(xl  +  r>x2>x3)  "  u15(x1,x2,x5 

u-l  = u-j_^  =  -  lim  * — s 

Sxx   3  5xx      3  r^0  r 


or 


i  lim  I 
3  r   \n  r 


0 


ul5(xl  +  ri>x2'x3)  "  ul5(xl>x2>x3) 


But  on  account  of  the  homogeneity  extended  to  the  averages  of  the 


third  order,   ui  (xi  +  r>x2>x^)     ^s   eclual  to  ul  (xl^x2^x3J-   Hence  the 
limit  is  zero.   Lastly,   c(r)   is  an  infinitesimal  of  the  third  order 
with  respect  to  r.   The  same  holds  for  b(r)  and  a(r). 

If  the  flow  is  incompressible,  the  various  functions  that  charac- 
terize the  correlations  of  isotropic  turbulence  are  not  independent.   From 

\ — ^a 
the  relation  2_ =  0  it  is,  in  fact,  immediately  deduced  that 

dx„ 


ZM  =  0  (9-6) 


This  equation  is  general.   In  the  particular  case  of  isotropic 
turbulence  it  leads  to  Karman's  equation 


r  bf 


g  =  :f  +  ^S  (9-7) 

2  or 
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In  the  same  manner,  one  then  finds  that: 


c  =  -2b 


a  =  -d  -  £  ■& 

2  dr 


(9-8) 


So,  of  the  five  functions  f,g,a,b,c  which  define  the  double  and 
triple  correlations  of  isotropic  turbulence,  only  two  remain  distinct  if 
the  flow  is  incompressible. 


The  incompressibility  results,  in  particular,  in  the 
relation  g"(0)  =  2f"(0).  Therefore,  in  terms  of  dissipation  length  A: 


f(r)  =  1  - 


2A2 


g(r)  =  1  -  4+ 
A2 


(9-9) 


Together  with  the  length  A,  two  other  numerical  parameters  (time 
functions)  are  frequently  used  to  give  an  idea  of  the  turbulence. 

First,  the  correlation  length 


L  =    f (r)  dr 

Jo 


(9-10) 


is  introduced;  it  depends  upon  all  values  of  f  (r)  for  0  <  r  <  °°,  while 
A  depends  only  on  the  form  of  the  function  f(r)  at  the  origin.  It  may 
be  pointed  out  that  L  is  expressed  by  means  of  the  scalar 


R 


=  y      ua(x)ua(x  +  £)  =  Uq  (f  +  2g).  Owing  to  the  incompressibility, 


,of 


f  +  2g  =  3f  +  r2i-  =  2f  +  -°-  (rf).   If  rf(r)  approaches  zero  when  r- 
dr        dr 


(f  +  2g)  dr  =  2L 


(9-11) 


NACA  TM  1377  ^1 

The  intensity  of  turbulence  in  one  direction,  that  of  the  axis  of  x-,, 
for  example,  is  the  dimensionless  quantity 

^~  (9-12) 

U 


where  U  is  the  mean  velocity.   The  total  intensity  of  turbulence  is  the 
ratio 


W^2  _  1  >/i(Q)  (9_15) 

u     3   u 


If  the  turbulence  is  isotropic,  its  intensity  is  the  same  in  every 

u  ' 


Up, 

direction.   It  is  equal  to  the  total  intensity  and  to 


In  the  case  of  isotropic  turbulence  the  mean  square  values  of  the 
rotational  components  are  directly  associated  with  the  quantities  Uq 

and  A. 

Introducing  the  components 

bu-z        dug  du-]_   duv  Sup   du-^ 

o>l  = -  av>  =  -  ou   =   -  

0X2   dx*  dxi   ax^  ox-^   Bxo 


O        p        p 

of  the  rotational,  the  averages  o^  ,     oy>   >     ^   ^OT   homogeneous  and 
and  isotropic  turbulence  are  computed: 


2   (^u^\  fin2\  ^u3  ^u2 
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We  note  that,  for  example: 

c*2   ^2^°  62 


'3u, 

=  lim  _i 

aX2/    &2-*°  ^22 


u32(xl>x2  +  ^2'x3)  +  ^^l'^'5^) 


2u5(x1,x2,x5)u3(x1,x2  +  62,x5)] 


The  first  tvo  bracketed  terms  are  equal  to  uq  .   The  third  is  equal 
to  R^v(0,|2j°)  =  Uq  6>  or>  except  for  the  third  order,   Uq  1 


Lastly 


Likewise 


'^1\   =  I__o! 

dxj     A2 


Sx5 J     *2 
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^3 


Lastly: 


lim- 


6263 


1*25(0, |2,  -  i3j  -  1^(0,0,13)  -  r25(o,i2,o\  +  r12(o,o,cm 


"o 


2  lim- 


t2h  [ 


f  -  gF  f 


2AC 


Finally 


2   ,,2 


^  „  2^0_  +  2^0_  +  ^0_  =  £0 


A2     A2     A2     A2 


2         2  2 

Naturally,  uy>       and  0)5   are  equal  to  a>^_    , 


10.  Spectral  decomposition  of  the  velocity: 

We  know  that  the  representative  curve  of  a  turbulent  velocity  compo- 
nent u(t)   in  terms  of  time  suggests  the  idea  of  an  irregular  periodic 
phenomenon.   The  analytical  representation  of  u(t)   is  not  a  periodic 
Fourier  series,  but  rather  an  almost  periodic  series  or  a  Fourier  integral. 

In  the  first  case,   u(t)   is  a  sum  of  harmonics  without  common  base 
period;  by  adopting  the  complex  notation 


u(t)  =  /    ^e 
n=-°° 


iofct 


(10-1) 


the  pulsations  <%  forming  a  succession  of  real  numbers  increasing 
with  n.   It  can  always  be  supposed  that  a>_n  =  -ci^.   u(t)   is  a  real 
quantity  if  A   =  \i* >    *^e  notation  A*  designating  the  conjugate 
imaginary  of  An.   The  series  An  must  be  convergent. 
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In  the  second  case: 


u(t)  =    A(tD)eia)tdo)  (10-2) 


the  function  A(cn)  being  absolutely  suramable  and  such  that  A(-o)  =  A*(a>) 

The  numbers  An  or  the  function  A(m)   depend  upon  the  position  of 
the  anemometer. 

This  representation  is  practical  when  the  turbulence  is  steady  in 
time,  which  enables  a  vigorous  application  of  the  ergodic  principle  and 
the  calculation  of  the  averages  connected  with  u(t)   in  a  time  interval 
as  great  as  desired.   But,  among  the  simplest  problems  involved,  there 
is,  first  of  all,  that  of  a  turbulence  homogeneous  in  space  (at  least  in 
a  sufficiently  restricted  range)  and  developing  in  time.   The  pseudo- 
periodic  character  of  the  speed  is  then  manifested  in  space  rather  than 
in  time.   (A  detailed  study  follows.) 

The  component  ua  of  the  velocity  of  fluctuation  at  point  x-^XgjX? 

at  a  given  instant  can  be  developed  in  series  or  by  Fourier  integral, 
depending  upon  whether  the  spectrum  is  discrete  or  continuous 

ua(x,t)  =         Z^       5'(t)e  V11      22      33    ' 


ni*°2'n3 


Ua(x,t)  =  fz^^yty^^^^)^    dA2  dA3 


(10-3) 


(10-4.) 


n-|_,   ng,  and  Hi  are  three  integers  that  vary  independently  from 

(W"5) 

-oo  to  +0=  in  formula  (10-3),  and  on  which  the  quantities  Z^  , 


a 


(n-^)  (ng)     [n,] 

^1    ,   A2    ,  ~\y      '      are  dependent.   In  formula  (10-4),   u(X  is  a 

triple  integral  extended  over  the  entire  space  A  of  wave  numbers  (The 
wave  number,  inverse  of  a  length,  is  the  equivalent  in  space  to  the  fre- 
quency, inverse  of  time.   It  is  a  vector).   ua  is  a  real  quantity  if 

Za(-\L>-A2>"A3>t)  =  Za*(Al'A2'A3't, 


NACA  TM  1377  ^5 


These  notations  are  of  classical  form.   But,  since  the  u   are 

steady,  random  functions  of  space,  it  is  preferable  to  represent  them  as 
stochastic  Fourier  integrals 

ujxjt)  =  I  ei(AlXl+V2+A3x3)dha(A,t)  (10-5) 


In  this  formula,  the  ha(A,t)  are  random  functions  of  A-^,^*^ 

with  orthogonal  increments  (or  noncorrelated) .   Or,  in  other  words,  if 
A  and  A'   are  two  different  points  in  the  space  A,  the  mean  of  the  pro- 
duct of  the  two  increments  dha(A,t)  and  dho(A',t)   is  zero.   On  account 

of  the  complex  notations  it  naturally  is  a  question  of  a  product  of 
hermitian  symmetry.   Hence 


dha*(A,t)  dh3(A*,t)  =0   if    A'  f   A  (10-6) 


If  A1  =  A  at  the  same  point,  this  average  is,  in  general,  no  longer 
zero.   It  is  infinitely  small  of  the  order  of  dA.  We  put 


dha*(A,t)  dhp(A,t)  =  cpap(A,t)  dA  (10-7) 


Now  the  spectrum  of  turbulence  in  the  different  cases  (10-3),  {lO-K) , 
and  (10-5 )  is  defined. 


11.  Spectral  tensor  and  correlation  tensor: 

Assuming  spatial  homogeneity  it  is  now  attempted  to  form,  from  the 
spectral  decomposition  of  the  velocity,  the  expression  of  the  compo- 
nents RaR  (^1^^2'^3,t)  or  abbreviated,  RaB^)  °^  the  correlation  tensor 
of  the  velocity  at  a  given  instant. 

If  the  formulas  (10-3)  °r  (10-4)  are  utilized,  the  spatial  averages 
in  a  very  great  volume,  physically  limited  but  practically  infinite  with 
respect  to  microscopic  turbulent  lengths,  must  be  used.   It  involves  a 
cube  with  edges  parallel  to  the  axes,  of  arbitrarily  large  length  2a,  of 
which  the  center,  which  can  be  placed  anywhere  because  of  the  homogeneity, 
is  placed  in  the  origin  of  the  axes  of  the  coordinates. 
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First,  take  the  case  of  the  Fourier  series  (10-3).   To  form  Raa, 
involves,  first,  the  product  ua(x)ug(x  +  |),  by  associating  an  arbitrary 
term  of  the  series  which  represents  u  (x),  or  rather  u  *(x),  to  any 

one  term  of  the  series  that  represents  Up(x  +  i) .      To  form  the  average 

in  the  volume  V  =  8a3,  involves  division  by  V  and  integration  of  xa 

in  space.   The  exponentials  are  restricted,  and  likewise  their  integrals 
in  V,  so  that  the  quotients  by  V  approach  zero  when  a — >°°.   There 
is  an  exception  for  the  terms  of  the  series  u a(x)   and  ua(x  +  | )  which 

correspond  to  the  same  wave  numbers,  and  which  give  the  products 

z#(n1,n2,n3)e   p=1  ^  ^       frl'^'^)  e  p=l  V  T  VJV         . 
a 

i  £  ipAp(np) 

(n1,n2,n5)z  (n1,n2,n5)e  p=l  (ll_l} 

a  3 


independent  of  xa.   They  are  therefore  equal  to  their  mean  value  and 
consequently 


3 


1 


iv\ 


(°p) 


Z_    Za  Z3V        e  (11-2) 


nx,n2,n5 


This  is  none  other  than  the  Parseval  formula  for  the  almost  periodic 
Fourier  series  of  three  variables.   It  naturally  assumes  conditions  of 
uniform  convergence  which  need  not  be  defined  here. 
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In  the  case  of  the  Fourier  integral,   ua(x)un(x  +  |)   is  a  sextuple 
integral 

3  3 

n  -1!-  Vp+1XI  *'p(*p+*p) 

ua(x)up(x  +  O   =    /       Za*(A)Z   (A')e       P  P  dAdA1 

Jaxa 

(n-3) 

extended  over  all  values  of  A^_,   A2,   A^,   ^'i>   ^*2>  an(^  ^*3' 

Integrating  under  the  sign  /  in  the  finite  volume  V  and  dividing 

by  V,  we  find  that,  mathematically  speaking,  the  mean  value  which  is 
being  sought  is  the  limit  for  an  infinite  of 


Za*(A)Z  (V)e  P=1      dAdA'   /  e  P_1         dx  = 
v  'AXA       P  Jv 


3 

A   /               fe   P  P      sin  A*-,  -  A,)  a 
2   /   Z*(A)Zq(A')e  P        dA  dA'  i-1 ±1— 


sin  (A1 2  -  A2)  a  sin  (  A'x  -  A-z 
A  o  -  ^2         ^'  3  ~  ^3 


(11-4) 


Suppose  Za(A)  satisfies  "Dirichlet's  conditions"  or,  more  gener- 
ally, is  a  limited  variation;  then  the  integral  which  figures  in  (11-4) 
has  a  finite  limit  which,  when  Za(A)   is  continuous,  has  for  value 

3 

p  1  X"  Vp 

8«3   /  Za*(A)Zp(A)e  P"1     dA  (ll-5) 
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This  is  Dirichlet's  theorem. 


The  mean  value  of  ua(x)un(x  +  |),  the  quotient  of  this  integral  by 

V  =  8a5,  thus  approaches  zero  when  a — ><*>,  and  the  formalism  of  the 
Fourier  integral  does  not  furnish  the  spectrum.  This  conclusion  is 
paradoxical . 

But  every  difficulty  disappears  when  it  is  assumed  that  the  turbu- 
lence might  have  only  a  line  spectrum,  incompatible  with  the  so-called 
Fourier  integral.   But  such  an  assumption  does  not  seem  reasonable. 

It  seems  more  satisfactory  to  concede  that  the  functions  representing 
the  turbulent  velocity  are  too  complicated  for  applying  Dirichlet's  theorem 
and,  more  precisely,  that  their  oscillations  are  too  crowded  to  be  repre- 
sented by  functions  with  bounded  variation.  A  mathematical  process 
avoiding  this  difficulty  is  to  follow. 

But  an  approximate  argument  can  also  be  made.   In  reality,   a  is 
very  great  (with  respect  to  turbulent  wave  lengths ) ,  but  finite .   Hence 
the  mean  value  of  uQ^(x)un(x  +  i)      is  given  by  formula  (11-4),  without 

transition  to  the  limit.   However,  a  transition  to  the  approximate  limit 
can  be  made  and  assumed  that,  with  suitable  accuracy 

3 

3  r  1 II  Vp 

8*     -7  W.W  hu   P=1 


So 


ua(x)up(x  +  I)  :  ^-  /  Za*(A)Zp(?0e  *  x     dA      (11-6) 


,  if  a  spectral  function  cpag(A)   is  defined  by 


Vp 


Rap(l)  =  ua(x)up(x  +  I)  =  /cpa(3(A)e  P  1  dA      (11 -7) 


we  get 


cpap(A)  =  2^a*(A)zp(A)  (11-8) 
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In  a  more  accurate  way,  the  9an(A)   constitute  the  components  of  a 

tensor,  the  spectral  tensor,  transformed  from  Fourier's  correlation  tensor. 
But  the  approximate  formula  (ll-8)  has  the  drawback  of  yielding  a  spectral 
tensor  of  a  too  restrictive  form.   In  fact,  there  is  no  reason  that  the 
cp  _  be  the  product  of  a  function  of  subscript  a  by  a  function  of  sub- 

script  3,  that  is,  the  general  product  of  the  vector    / Z.  (A)  by 

V  V 

itself.   This  new  defect  can  be  corrected  by  superposing  a  temporal  average 
on  the  spatial  average  (ll-8).  But  it  is  much  preferred  to  turn  elsewhere 
and  have  recourse  to  the  representation  of  the  velocity  by  Fourier's  sto- 
chastic integral  (10-5).  The  product  u  (x)u„(x  +  i)     is  a  double  sto- 
chastic integral 

3        3 

p    -1  5Z  Vp+1  ^  *'P(xP+£p) 

ua(x)up(x  +  I)  =  /   e   P=1  dha*(A)dhp(V) 

(11-9) 


On  averaging,    it   is   seen  that  the   quantity     dh  *(A)dh   (A')      is 

involved  which  is,    as  explained  earlier,    zero  when     A'    /  A,   and  of  the 
form     cpao(A)dA     if     A'    =  A.      Hence   simply 


1  XI  Vp 

Ra3^)   :=  UJX)VX  +   5)    :      /e     Pl  <Pap(A)dA  (11-10) 


The  correlation  tensor  Ran(l)   Is  "the  Fourier  transform  of  the 
spectral  tensor  <pag(A).   The  <PaQ(^)  here  have  the  generalty  desirable, 
They  are  subject  to  the  two  following  conditions: 

(a)  cp  R(A)  has  hermitian  symmetry 


V  =  <pPa; 


* 
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(b)  Whatever  the  complex  numbers  X  ,  the  hermitlan  form 

3 
a,p=l 

which  is  real  according  to  the  first  condition,  cannot  be  negative. 

Any  system  of  numbers  cp  R  constituting  a  tensor  and  satisfying 
these  two  conditions  can,  a  priori,  be  used  as  spectral  tensor. 


12.  Spectral  tensor  of  isotropic,  incompressible  turbulence: 

I.  Suppose,  with  Heisenberg  (reference  23),  that  the  turbulence  is 
isotropic.   Then,  as  supplementary  condition,  the  flow  is  incompressible. 

The  method  used  to  express  the  isotropy  of  the  correlation  tensor 
is  applied  to  all  tensors,  whatever  its  meaning.   The  spatial  tensor  is 
isotropic  when  two  functions  A(k),   B(k)   exist  depending,  moreover,  on 
t  such  that 


<Pap00  =  A(k)AaAp  +  B(k)5ap  (12-1) 

2  2     2     2 

k   denoting  the  quantity  A^_  +  \^     +  A^  • 

If  the  flow  is  incompressible,  the  functions  A  and  B  are  formed 
by  a  relation  equivalent  to  Herman's  relation  (9-7)  between  the  functions 

<c —  du_ 
of  the  correlation  f  and  g.  We  write,  in  fact,  that  >   — =  °- 

dx_ 


a 


Since 


3 

N. 

n1  S_  Vp 

a  _ 
a 

1  J  a     ^            ^~  Aadha 

(12-2) 
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we  must  have   identically 

IW  =  0  (12-3) 

hence,   with     p     being  fixed 


± 


.  VPapM    =  ° 
a=l 


When  this  condition  is  applied  to  the  components  (12-1)  of  the 
isotropic  spectral  tensor,  it  is  seen  that 

Ak2  +  B  =  0 


We  put 


A  =  --E —     B  =  -£- 


jtk4        knk2 


hence 


F(k)  /     A  A3\ 


F(k)   is  the  spectral  function  introduced  by  Heisenberg.   It  will  be 
interpreted  later. 

II.  With  Kampe  de  Feriet  (reference  26)  we  first  express  incompressi- 
bility  without  disturbing  the  isotropy. 

The  cp  n  are  subject  to  the  condition  that  the  hermitian  form 
Tap 


+  =  ZW  ^aP  (12"5) 
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is  positive  or  zero.   It  is  positive,  except  when  the  vector  X   is 


proportional  to  vector  A^,  since 


5Z  V^o-P  =  ° 
a(3 


(12-6) 


on  account  of  the  incompressibility  condition  (12-3)- 

The  hermitian  form  \|r  is  now  reduced.   According  to  classical  theory 
and  the  equation  in  S,   three  real  numbers  S^,S2^S^  and  three  linear 

forms 


Yn  = 


\2*  Y2  =  YL   bccXa     Y3  =  H  caXa      (12"T) 


can  be  found  with  which  \|r  can  be  expressed  in  the  form 


*  =  S1 


Yl   +  S2  Y2 


+  S, 


y. 


the  three  vectors  of  components  aajtla,ca  being  "two  DY  "two  orthogonal. 
If  \|r  >  0,  the  numbers  S-^^So^S^  are  positive  or  zero. 

But,  there  exists  a  particular  vector  7  ,  transformed  from  vector  A 

by  (12-7),  so  that  ^  =  0.   This  can  happen  only  when  the  three  terms  of 
which  \|r  is  the  sum,  are  zero  at  the  same  time.  Hence  two  separate 
hypotheses 


(a)  So  =  0 


s5  =  0 


V\x  =  ° 


t  =  S, 


5~~aaXa. 


Pa(3  =  Slaa*a(3 


(12-8) 
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a       being  the   function  of     A.      This   form  of     m  a,    similar  to   (11-8),    is 
too  restrictive; 

(b)     S5   =  0  Z>a\*  =  °  XIVa  =  ° 

*   =  Sl  IZ^^I2  +  S2  |l>aXa|2  (12"9) 


We  write  rather 


*  =  si2_  Wxp  +  s2  2_  VbpVxp 


and  identify  with  the  original  form  (12-5)  of  i|r.   It  is  seen  that 


<PaP  =  Slaa*a3  +  S2ba*bp  (12"10) 


This  is  the  most  general  form  of  a  spectral  tensor  in  an  incompress- 
ible medium.   In  more  simple  form 


S,a  =  a'       ,/S0b  =  b' 
la     a      V  2  a     a 


There  are  two  vectors  a'  ,     b'   forming  with  vector  ~h       a  trirec- 
tangular  trihedral  and  such  that 


<Po3  =  a'a*a'3  +  b'a*b'P  i12-11) 
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By  Pythagorean  theorem 


ZW2  =  ^rva 


2  +  -1- 


.2 


Za'aXa 


2  +  ^ 


b" 


Lb'aXa 


(12-12) 


with 


.2 


■  H    '".' 


b'2  -    V~  b'   2 

Z —         R 


By  this  formula,  the  b1   can  be  eliminated  from  the  expression 


Of   CD 


op 


^M-SK-v^L-f 


(12-15) 


Finally  we  put 


(.2  =  F(k) 

hnk2 


'l  "  *|  *a  -  ca 


This  is  the  canonical  form  which  Kampe  de  Feriet  has  given  for  the 
spectral  tensor  of  an  incompressible  fluid 


Vx) 


=  c  *(A)ca(A)  +  n±l[b 


i±)L 


~<x  x  '  3 


c 

4rtk2\ 


ap  "  v2 


(12-U) 


with 


Zca(^)\c  =  ° 
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The  isotropy  stipulates  that  ca  =  0,  so  that  formula  (12-^)  is 
obtained  again. 

15.  Energy  interpretation  of  the  spectral  function  F(k): 

Except  when  stated  otherwise,  the  isotropic  turbulence  is  assumed 
to  be  in  the  incompressible  state . 

Hence 


2     2     2     2 

ul  =  "2   =  u3   =  ^ 


The  total  energy  of  the  turbulent  fluctuations  per  unit  mass  is 


(15-1) 


But 


9mw  -  m  (15-2) 

2«k^ 


Hence 


E   =    1  1M.  dA  (15-3) 

k«  ,  A     k2 


In  passing  to  the  polar  coordinates  in  the  space  of  the  wave  numbers 
one  finds  that 


E  =    F(k)dk  (13-M 

P 
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Which,  after  minor  changes,  gives  the  amount  of  the  energy  of  the 
turbulent  fluctuations  per  unit  mass  in  the  "sphere  of  the  wave  numbers" 


Z 


Aa2  <  k2  as 


Ek  :   I  F(k')dk' 


(13-5) 


k '   being  an  integration  variable . 

The  "small  values  of  k"  correspond,  of  course,  to  the  large  scales 
(large  vortices). 

The  next  problem  is  to  find  the  expression  of  the  energy  e  dissi- 
pated as  heat  per  unit  mass.  While  the  preceding  calculations  were  by 
nature  strictly  kinematic,  the  dissipative  function  here  must  be  so  intro- 
duced that  it  yields  equations  of  motion,  that  is,  Navier  equations; 
hence,  an  assumption  associated  with  the  dynamics  of  turbulence,  which 
is  discussed  in  chapter  III. 

With   v  denoting  the  coefficient  of  kinematic  viscosity,  the  energy 
dissipated  in  heat  by  the  molecular  motion  is 


'duA2 


>oX 


chiz   dug 


(13-6) 


This  is  to  be  expressed  in  spectral  terms  and  averaged.   If 


a 


1 YZ  Vp 


P=1     ^aM 


we  get 


du, 


a 


c*r 


i£ 


i   V  P-l  dha(A) 
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and 


c>ua       dug 

+  — -  =1     /  e 

dxp        Sxa  .  /A 


1  2_  Vp 

'*p  ^a  +  V  top 


which,    squared,    reads 


Sua       ^p' 
dxp        axa/ 


1 5  M* 


a3  dha*(>,)  + 


Aadhp*(A)      A'p   dha(A')   +  A'adhp(A') 


(13-7) 


To  obtain  the  averages,    involves   operation  under  the    j   sign.      The 

result   is   zero,   except  when  the  points     A     and     A'      coincide   in  the 
space     A     of  the  wave  numbers .      It  leaves 


Sxp     av 


A, 


dh 


a 


+  Vp 


dha*    dhp  +   dh3*   dha 


+  A 


a 


dhr 


/  [Ap2cPaa  +  Vp(<Pap  +  ^Pa)   +  \x2<Ppp]d* 


(13-8) 


Consequently,   if     e     is  the  mean  dissipation  of  energy  per  unit  mass 


e  =  v 


Y.    V  +  JZ  Vp<PaP 


a 


aP 


dA 


(13-9) 
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This  formula  supposes  neither  incompressibility  nor  isotropy.   If, 

now,  the  turbulence  is  incompressible,  2. ^c^R^aB  =  ^"   Owing  "to  "the 

ap 
F 


isotropy,   2__  V 

aa       2«k 

Hence,  by  integrating  in  polar  coordinates  in  space  A 


e  =  2v  I  k2F(k)dk  (13-10) 

JO 

e  being  a  physical  quantity,  hence  finite,  it  is  seen  that  the 
preceding  integral  is  finite,  which  gives  a  first  limitation  of  the 
possible  forms  of  the  functions  F(k). 


Ik.   Relations  between  spectral  function  F(k)   and  correlation  func- 
tions f(r)   and  g(r): 

The  spectral  tensor  q>  „  was  defined  by  the  formulas 


Rap(h)  :;  ua(x)up(x  +  h)  :  hafiMe  ^ 


A 


which  assume  homogeneity  only.   Incompressibility  and  isotropy  are  to  be 
added. 


According  to  section  8 


W6)  =  uo2 


Sia£a   +  go, 


"Tp^R  -  B^ap 
rc 


g  =  f  +  -  — 
2  dr 


and  according  to  section  12 


m      /,,    F(k)  L  AaAp\ 
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By  Fourier's  reciprocity  formulas 

3 

n  -1  Yl  Vp 

^R(A)=^-T        R.R(l)e         P=1  ^  <U-1) 


the  integral  being  extended  over  the  entire  space .  The  general  formulas 
for  defining  F  in  terms  of  f  or  g,  and  f ,  g  in  terms  of  F  must  be 
particularized. 

F    <\ — 
It  is  pointed  out  that  =  >  cp    is  the  scalar  invariant  of  the 

spectral  tensor  cp  „   ("trace  of  tensor"  or  contracted  tensor).   Hence  it 

seems  logical  to  compare  with  F  the  analogous  invariant  of  the  correla- 
tion vector  R  _ ,  defined  by 

R(r)  =  Y_   Raa  =  "o^*"  +  2g)  =  uQ2(3f  +  rf '  > 


Thus  the  following  problems  must  be  solved: 

(a)  express  F(k)  by  means  of  R(r),  and  conversely 

(b)  express  f (r)   and  g(r)  by  means  of  R(r) 

(c)  express  f(r)   and  g(r)  by  means  of  F(k) 

Also  to  be  defined  in  precise  manner  is  what  is  called  the  trans- 
versal (lateral)  and  the  longitudinal  spectrum  of  turbulence. 

(a)  Relation  between  F(k)   and  R(r): 


3 


i   >       A  % 

il     p  p 

R(r)   =  IcpcGacLA=     |e     P"  ^1  dh  (U-2) 

2rtk2 
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To  reduce  the  triple  integral  to  a  simple  integral,  simply  select 

a  new  axis  of  the  ?w  perpendicular  to  the  plane  }      A  £   =  0;  it 
results  in  the  integral 


R 


(r)  =  [eirA3  !M  dA  =  ±    Teirk  COS  9F(k)sin  9  dk  d6  dq> 
JA     2nk2      2jt  Ja 


Which,  after  a  short  calculation,  yields 


Conversely 


or 


jL/r3f\    =  & 


(U-3) 


B(r)=2         H^F(k)dk  (14-4) 


^1  =  -         r  sin  rkR(r)dr  (l4-5) 

k  it  Jo 

(b)  Relations  between  f(r),   g(r),   and  R(r): 
The  integration  of  the  differential  equation  yields 


rf •  +  3f  =  JL.  (14-6) 

uo2 
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with  the  condition  f(0)  =  0.   It  is 


UQ2f(r)  =  —    r'^r'Jdr' 
r5  Jo 


(14-7) 


R(r)  having  for  limit  3^  t  when  r  tends  toward  zero,  it  con- 
firms that  the  integral  is  equivalent  to  Uq  r*,  and  that  f (r)  tends 
rather  toward  unity. 

From  that,  it  is  deduced,  by  formula  (9-7)  that 


u02g(r)  =  iR(r)  -  -i-  I  r,2R(r')drI 
2       2r5  Jo 


(1^-8) 


(c)  Relations  between  f(r),   g(r),  and  F(k) 
A  simple  calculation  gives 


Uo2f(r)  =  2    F(k) 


'sin  rk   cos  rk \ 


'3k3     r2k2  I 


dk 


un2g(r)  = 


l^W 


sin  rk   sin  rk\,, 
+  dk 


~3k3 


rk 


(14-9) 


From  that  the  first  terms  of  the  limited  developments  of  R(r),   f(r), 
and  g(r)  for  small  values  of  r  are  derived  in  terms  of  the  dissipated 
energy  e 
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R(r)  =  3u^(l  -  ^r' 


1J2.    e 


^0 


vuq 


f(r)=:l   ^r2_±_  + 

VUqC 


30   „-2 


g(r)  =  1  -  -Lr2— S—  + 
15     2 


(14-10) 


Comparison  with  (9-9)  indicates  that  the  dissipation  length  A  is 
related  to  the  dissipated  energy  e,  the  viscosity  v  and  the  kinetic 

energy  of  the  turbulent  fluctuations  E  =  ?^q       through  the  formula 


F(k)dk 


A2  =  10v-  =  5- 

e 


k^F(k)dk 


(14-11) 


15.  Lateral  and  longitudinal  spectrum: 


The  spectral  measurements  do  not  provide  the  spectral  function  F 
directly,  but  merely  the  spectral  terms  corresponding  to  certain  simple 
associations  of  velocity  components,  for  the  observers  placed  in  the 
particular  relative  positions .   It  concerns  longitudinal  components  of 
the  velocity  (parallel  to  the  velocity  of  the  main  flow),  at  points 
placed  either  on  a  parallel  to  the  axis  or  on  a  perpendicular  to  the 

The  corresponding  correlations  are  UQ2f(r)  an^  UQ2g(r). 


axis, 
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According  to  the  original  notations  of  G.  I.  Taylor 

\ 

J "i  00 
cos  asrA(co)da) 
0 


un2g(r)  =  J   cos  asrB(a))da> 
0 


"0 


(15-1) 


where  A(u>)   and  B(cjd)  are  termed  the  longitudinal  and  lateral  spectral 
functions.  They  now  must  be  expressed  by  means  of  F(k). 

Since  g  =  f  +  — f ' ,  elementary  calculation  shows  that 


B(cjq)  =  ^[a(o))  -  ooA'(u))] 


(15-2) 


and  at  the  same  time  that 


^00 


R(r)    =  2 


Sin  r^F(k) 


dk 


0  rk 


loo 

cos  ojrfA(o))  +  2B(o))J  dm       (15-3) 


Elementary  Fourier  transformations  yield 


A(o>)   +  2B(cd)  =  2    ^Uk 
Jo      k 


(15-4) 


A  and  B  are  determined  by  the  two  equations 


2B  =  A  -  aA.' 


A  +  2B  =  2 


F{k] 


dk 


(15-5) 
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The   result   is 


AM .  riM  l  .  & 


dk 


ia 


"-■J^H 


dk 


(15-6) 


CHAPTER  III 


DYNAMICS  OF  TURBULENCE 


l6.  Introduction: 

To  construct  a  dynamics  of  turbulence,  that  is,  to  set  up  the  dif- 
ferential or  finite  laws  which  govern  the  development  of  the  statistical 
quantities  characterizing  the  turbulence  in  time  and  space,  it  is  neces- 
sary to  start  from  elementary  laws  and  apply  the  statistical  methods  to 
them.   The  most  natural  idea,  the  only  one  which  actually  produces  con- 
crete results,  consists  in  utilizing  the  Navier  equations.   To  what 
extent  are  they  applicable  to  turbulence?  The  turbulent  motion  is  always 
a  macroscopic  motion  with  respect  to  a  finer  scale  motion,  and,  at  the 
limit,  with  respect  to  the  molecular  disturbance.   Therefore  it  is  rea- 
sonable to  believe  that  it  satisfies  the  equations  of  the  mechanics  of 
fluids.   The  next  step  is  to  find  the  solution  of  the  Navier  equations 
which,  for  certain  limiting  conditions,  have  the  turbulent  aspect,  and 
calculate  the  particular  averages  from  these  solutions.   Unfortunately, 
it  is  rather  difficult  to  define  these  limiting  conditions.   So,  the 
remaining  resource  is  to  examine  whether,  among  all  the  possible  solutions 
of  the  Navier  equations,  there  exist  any  sufficiently  complicated  for  rep- 
resenting the  turbulence,  without  attempting  to  determine  them  logically 
by  the  limiting  conditions.   But  then  a  new  difficulty  arises.   Solutions 
for  the  Navier  equations  are  known  only  for  simple  conditions  which  are 
far  from  resembling  turbulence.   In  other  words,  while  conceding  their 
validity,  we  practically  do  not  know  how  to  solve  them. 

Since  it  is  not  acceptable  to  take  the  averages  on  the  solutions  of 
the  Navier  equations,  it  is  attempted  to  take  the  averages  on  the  differ- 
ential equations  themselves  and  to  write  the  differential  equations  which 
verify  the  statistical  quantities.   Since  this  method  has  produced  results 
it  is  set  forth  in  the  simplest  case,  that  of  the  Karman  correlation 
tensor. 
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17.  Fundamental  equation  of  turbulent  dynamics: 

It  is  assumed  that  the  total  velocity  U  is  uniform  and  along  the 
axis  of  x-,  .   The  Navler  equations  read  then 


^a  +  ^  +  £"_  Up^  =  _1  ^_  +  v  AUa         (17.1} 


St     Sx-l    p    dxp    P  dx 


a 


in  the  absence  of  external  forces,   p  is  the  density,  usually  assumed 
constant,  and  p  the  pressure.   These  equations  must  be  supplemented 
by  the  equation  of  continuity 


a  =  0  (17-2) 


3xa 


The  notations  are  abbreviated  by  representing  the  velocity?  at 
point  x'   =  x  +|   by  u'   =  u  (x' )  =  u  (x  +  £ ).  Multiplication  of 

(17-1)  by  u1   and  summation  with  respect  to  a,  gives 


—    du      ^ —    du_   ^ du„ 

_  u'  — £  +  U  ^   u'  — a  +  )   u'auQ — a 
a     dt      '  a    aaXl    ap     pdxp 


.       Sp  r d      U  . 


1 


a^xn.  as         \R2 


ca  a3  dxp£ 


'•In  chapter  I,  section  6,   u^,u^,u7  represented  the  velocity  compo- 
nents of  the  whole  and  u'-^u^u1*,  those  of  the  velocity  of  fluctuations. 

The  over-all  velocity  having  then  U,o,o,  as  components,  the  notations  are 
changed;  Uj_  represents  the  velocity  fluctuations  at  point  x,  and  u'j_ 

the  velocity  fluctuations  at  point  x ' . 
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Permutation  of  the  points     x     and     x'      yields   a  similar  equation, 
which  is  added  to  the  preceding  one 

d     t —  ^ —    /       Su„  du'   \       v —    /  Su„  du'    \ 

—  Y^  u„u'    +  u  >       u'— a  +  u — a    +  >       u'  uR— a  +  u  u'R— -a    = 


cR2       axy 

(17-4) 


lr-f.^p               dp'     \            V-          ,5\             ^u'^ 
—     >        u  (xr-^— +   u  — +    V    >         u'a 2  +  Uqt a 


Taking  the  averages  of  the  two  members  it  is  then  assumed  that  the 


turbulence  is  homogeneous,  nothing  more.   The  scalar  product  2 uau'a 

is  a  function  of  £a,  and  is  not  individually  dependent  on  x   and  x'  . 
Consequently 


a   v~  11  n«    =     a   v~  „  ,,1  a v     ,.  ,,i    _   a 


\    ii  ii  i   —    B   \    11  ii'  B   \    ii  n  '   _   o   \    ,,  „i 

Z_  uau  a  -  "  TT   2_  uau  a      77  ^—  UaU  a  ~  TT "  2—  uau  a 
dx-L  d^  ax'x  di1 


The  term  containing  the  over-all  velocity  U  disappears.   The 
following  term  can  be  written 


u  a7-^aup  +  Z_  ua7-7^  au  p  "IrZ  Uau  a 
op     dxp       ap    dx  p         3"  axp  -a" 


>   —     u  u 

z >  ,   < a  a 

3  °x  3  a 
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The  last  two  terms   disappear  as   a  result  of  the  equation  of  continuity. 


YL    ua%TuauP      is  bitten:8 


Ldx~  a  P 


P     dxp     a  0     o|d     a 


Likewise,   as   a  result  of  the   homogeneity,       /      u u1    u'„     becomes 

ox  p 

H  ^T  H  V^V^  =  +   II  ^-  H  W  "5  ) 
p     ^'p   V  P  V  ^p     a       ^ 

Since  T aft-y(-£)  =  "^aB7^)'  ^e  system  °^  these  two  terms  is  equal  to 

P   ^3  a 
The  viscosity  term  reads 

P   d|    a  a 

where  A  is  the  Laplacian  symbol  in  three-dimensional  space. 

Lastly,  it  is  shown  that  the  pressure  terms  disappear.  This  is  done 
by  extending  the  definition  of  the  homogeneity  to  the  averages  containing 
the  pressure,  with  due  regard  to  the  incompressibility.  The  averages  such 

as  pu'   are  not  functions  separately  of  x  and  x'   but  merely  of  the 

difference  £  =  x'  -  x,  so  that,  on  these  functions,  t-2—  =  -  °— . 

'  dx'     dx 

See  the  definition  and  the  properties  of  the  tensors  R  R  and 
T     in  chapter  II,  sections  8  and  9- 


JJC  -    z 


-    -  -    - 


---- 


-  7-    -    :      -  ?o=   '    -  —  :->- 


-_. - 


.    -         i  -        :;:;r-L::    i 

—  =  __•-_:  _ "  - : : 


-     -  _  _:  :  _  _  _ 

=    —    -    ""I 
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According  to  the  properties  of   r_o,   se;xion  9): 

— 


-5  =-;%=-^   =  Vj£  -  +^ 


(18-1) 


":.r  iivargeroe  ::  this  vector  is  er_ 


-2  cr 


r- 


-1 


(16-2) 


.    r    fundaaental  equation     -   -1C      la   now  transcribed; 


ri    rr|_     at      '       _  1    _r  :,     ;r[     \  r 


In  fact,   it   is   iuu  latlacian  of  a  function     <p(r)     in 

-~1~er.sisr.al   s:i:e   is 


-:-  -    =  ~ 


b  '  "  *') ' 


(16-3) 


Multiplica'.izr-    by      r~      is  follovei    by   integration  with  resres.    to      r, 
I"    is    eas:._.    verifiei   '.a.',    ill    _r.e    ter rs    :ar-sel   rut    for     r  =  C.      Searing 


^r?  =  u_2_£~_  _i_  r::-   •  =  _^r;  f  +  if 


rr 


-:-) 


cr.e  arrives  at  _r.e  equation 


|r  '^2-  =  --/-  -;-')  *  -o5(*'  -  f)       0--5) 
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This  is  the  Karman-Howarth  equation  (ref .  30)  in  its  classical  form. 

The  general  equation  (17-IO)  has  the  form  of  the  equation  of  heat 
propagation  in  ordinary  three-dimensional  space,  the  quantity  being 
propagated  is  the  scalar  R,  the  invariant  of  the  double  correlation 

tensor.   The  "second  member"  div  T  is  tied  to  the  triple  correlations 
as  R  is  to  the  double  correlations .  The  result  is  an  equation  of 

partial  derivatives  between  two  unknown  functions  R  and  div  T.   For 
reasons  which  appear  later  on,  R  is  usually  considered  as  the  principal 
unknown.   The  problem  of  solving  (17-10)  consists  in  defining  R  by  means 
of  suitable  physical  hypotheses  on  the  triple  correlations.  The  equation 
with  two  unknowns  finally  arrived  at  is  the  result  of  the  nonlinear  char- 
acter of  the  Navier  equations.   Its  mode  of  operation  will  be  explained 
later. 

The  "isotropic"  equation  (18-5)  leaves  no  trace  of  the  original 
properties  of  tensorial  symmetry.   Its  purpose  is  to  connect  the  func- 
tions f  and  c,  which  can  be  measured  directly.   If  it  is  borne  in  mind 
that 


f"(r)  +  ^f(r)  =  ^f(r) 


■♦*-?  !?(**•) 


(18-6) 


then  one  may  write : 


It  is  an  "equation  of  heat,"  isotropic  naturally,  but  in  a  fictitious 
five -dimensional  space.   This  remark  is  interesting,  but  it  seems  arti- 
ficial and  without  real  physical  significance. 

Nevertheless,  an  unusual  property  of  the  functions  UQ(t)   and  f(r,t) 

can  be  derived  from  it.  Multiplying  by  r^dr  followed  by  integration 
from  0  to  00  yields 
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^in2  A  dr  =  0  (18-8) 

St  °   Jo 

k    '  k 

provided  that     r  f       and     re     tend  toward  zero  when     r — »°°,   which  seems 

reasonable,   and  difficult  to  check.      Consequently 


X, 


uo^(t)         rMr,t)dr  (18-9) 


is  a  numerical  constant,  independent  of  t,  during  the  development  of  the 
turbulence.   It  is  the  Loitsiansky  invariant,  which  connects  the  energy 

of  fluctuation  E  =  ^-uq   to  the  correlations  f . 

It  also  should  be  noted  that,  since  the  second  term  of  (17-IO)  is 

the  divergence  of  2v  grad  R  +  T,  its  integral,  on  a  sphere  of  center  0 

and  infinitely  large  radius,  is  zero.  Consequently,  the  integral  of 

dR   . 

^—  is  zero,  and 

dt 


R  dxx  dx2  dx^  =  Cte  (18-10) 


the  triple  integral  being  extended  over  the  entire  space .   If  the  turbu- 
lence is  isotropic,  it  leaves 


/  r2R(r,t)dr  =  Cte  (18-11) 


But,  according  to  the  expression  of  R(r),  remembered  at  the 
beginning  of  this  paragraph: 


I  I-2R  dr  =  Uq2  r5f   =  Uq2  lira  (r^f  ) 
10  0       r — >oo 
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If  r'^f ' — >0,   r^f  tends  toward  zero  also,  and  this  limit  is  zero 

"loo 

r2R  dr  =  0  ( 18-12) 

Jo 

p 

This  proves  that  the  function  R,  which  is  equal  to  3Uq   for 

r  =  0,  and  to  zero  for  r  =  00,  takes  negative  values  for  sufficiently 
high  values  of  r.   If  formula  (18-12)  is  approximated  to  formula  (ik- 5) 
which  expresses  F(h)  by  means  of  R(r),  it  will  be  found  that,  pro- 
viding R(r)   tends  sufficiently  fast  toward  zero  at  infinity,  the 
development  of  F(k)   in  powers  of  k  starts  with  a  term  in  k^;  near 
k  =  0   (large  size  eddies) 

F(k)  =  Ck^  (18-13) 

with 

J~>00 
r^B('r)dr 
0 

The  constant  C  can  be  expressed  by  f,  because,  according  to  (l4-7) 


■oo  noo  Or 

Uq2     I     r*"f(r)dr  =  r  dr  rl2R(r')dr* 

Jo  .0  Jo 


In  order  to  change  the  order  of  integration,  the  integrals  of  the 
second  member  are  written  in  the  form  (fig.  5) 

;nr"  pT  Pr' 

r"   dr"  r,2R(r')dr»    =     lim  r|2R(r')dr'  r  dr 

Jo  r-"mj0  Jr" 


ir 


i     lim     r2    /     r,2R(r,)dr' 

Jo 


2  r — >oo 


i  r'^Ci-'JcLr' 


2  do 
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Now: 


r2  I  r,2R(r')dr'  =  u^r^r) 


So,  if  f(r)   tends  toward  zero  at  infinity  fast  enough  so  that 
not  only  r^f  but  even  r^f •  approaches  zero,  then 


2  P°° 
C  =  i£_    r4f  or  (18-14) 

5*  Jo 

To  a  factor  — ,   C  is  therefore  identical  with  Loitsiansky 's 
3* 
invariant  (18-9). 

These  results  could  be  extended  to  include  the  nonisotropic 
turbulence . 

Lastly  it  should  be  noted  (ref .  2)  that,  in  the  case  of  isotropic 

turbulence,  the  vector  T  defined  by  (l8-l)  is  radial,  that  is  to  say, 
that  the  TQ  are  proportional  to  the  |_.   In  effect: 


Uq5(c  +  2a)  /ox 


Hence  one  may  put : 


-4 
T  =  T  grad  r  (l8-l6) 


T(r)  =  Uq  (c  +  2a)  being  a  scalar. 


19.  Local  form  of  the  fundamental  equation: 

The  problem  is  to  ascertain  what  the  equation  (17-10)  becomes  when 
r— *0. 
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It  is  knovn  that  f(r) — >1.   As  f  =  1  -  -£—  +  .  .  .,   f"  +  ±f « 

2A2  r 

5 
approaches  — — . 

A2 
Lastly,  as  c  is  of  the  third  order  with  respect  to  r,   c'  +  — 


r 


tends   toward  zero. 

The   final  equation,    due   to  G.    I.   Taylor,    reads 

2 
d„  2  _      nn,u0 


-9^  =   -10v-^-  (19-1) 

itO  A2 


Recalling  that,  according  to  paragraph  Ik: 


u-,2  =  §-E     A2  =  10v£ 
^    3  e 

The  formula  (19-I)  becomes  therefore^ 


It  establishes  an  elementary  relationship  between  the  dissipation 
of  energy  e  by  viscosity  and  the  decay  as  a  function  of  time  of  the 
energy  of  fluctuation  E;  the  turbulent  energy  of  fluctuation  is  totally 
dissipated  as  heat  by  viscosity,  this  dissipation  following  a  process  to 
be  analyzed  in  chapter  V. 

The  foregoing  demonstration  stipulates  isotropic  turbulence.   But 
the  result  (19-2)  is  valid  in  more  general  cases.   In  fact,  (19-2)  can  be 
derived  from  (17-10).  First  it  is  known  that  R(0)  =  2E.   The  connection 
between  A  R  and  e  is  established  by  the  intermediary  of  the  spectral 
functions  cp  R(A).   If  the  flow  is  incompressible,  formula  (13-9)  reads 

9  r* 

E  being  function  of  the  sole  variable   t,  ^-  can  be  replaced  by 

d 

— ,  a  notation  which  represents  an  ordinary  derivative  with  respect  to 
dt 

time  (not  a  derivative  with  respect  to  motion  which  would  be  without  sense 

here ) . 
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e  =   v 


\  k2  5Z  U^)^  (x9-3) 


but  on  the   other  hand, 

3 


p  i  II  Vp  v- 

!(0   =  Z_  WD    =    /  e  >       qw(*)** 


consequently: 

3 


Vp 


AR  =    -     /  k2e     P  X  5Z  WA)dA  C1^) 


and  in  particular,  when  f    approaches  zero,   vAR  is  reduced  to 
-v  Ik2  ^  cpQa(A)dA  =  -e. 

To  prove  that  div  T  =  0,  at  the  limit,  the  simplest  way  is  to  revert 
to  the  Navier  equations . 

It  is  plain  that  equation  (19-2),  limit  of  (17-IO)  when  r — >0  can 
be  obtained  without  going  by  way  of  the  finite  values  of  r.   Simply  mul- 
tiply the  Navier  equations  by  u-^u^u,,  add  them  up  and  average.   Assuming 

the  disappearance  of  the  pressure  terms  to  have  been  proved,  it  is 


~   —    ^    Urr2    +     Y~  U«Uo^    =     V    V~    U    — -Jfifc  (l9"5) 

2  dt  /L~     a  z—    a  pSxo  ^x  2 
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The  term  i.  -^-  \   u    is  written  ^L  The  second  member, 
2  dt  ^—  a  dt 

according  to  the  foregoing  proof,  is  equal  to  -e.   Hence  it  must  be 


diL 


v —    ou 

shown  that   >   u  u — ^  =  0.   This  quantity  reads 


2  V  PaxP      a    2  ^3  a  3   2  a   a  V  dx3 


(19-6) 


p 
By  virtue  of  the  homogeneity,   ua  ur  is  not  dependent  on  x,  and 

the  first  term  is  zero . 

Owing  to  the  incompressibility,  the  second  term  is  zero.   The  first 
member  is  therefore  zero,  which  provides  the  proof. 


20.  Solution  of  the  fundamental  equation,  when  the  triple  correlations 
are  disregarded: 

When  the  triple  correlations  are  discounted,  equation  (17-10)  becomes 
the  equation  of  heat 


^2.  =  2v  £R  (20-1) 

at 


This  hypothesis,  which  mathematically  is  convenient,  is  physically 
quite  difficult  to  justify.   Its  accuracy  increases  as  the  viscous 

effects  become  greater  (prevalence  of  term  2v  AR  over  the  term  div  T  ), 
or  as  the  Reynolds  numbers  decrease  (whatever  their  definition).   It  also 
would  be  verified  if  the  velocity  components  followed  the  Laplace- 
Gaussian  law.   And,  what  is  more  interesting,  it  may  be  added  that  it 
corresponds  to  a  state  of  turbulence  in  which  the  forces  of  inertia  are 
negligible  against  the  forces  of  viscosity. 


Equation  (20-1)  must  thus  be  resolved  knowing  R  for  t  =  0,  and 
I   aware  that,  if  r — >0,   R 
temporarily  considered  as  known. 


being  aware  that,  if  r — >0,   R  has  for  limit  a  finite  quantity  3uq  (t) 
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In  isotropic  turbulence,   R  depends  only  on  r,  t,    and  (20-1)  reads 


dR   _  /d2R   2  dR\ 

^—  =  2v  +  -  t—  • 

dt     U.2   r  dr. 


(20-2) 


To  satisfy  the  limiting  conditions,  elementary  solutions  of  the  form 


R(r,t)  =  R1(t)R2(s)      s 


V8^t 


(20-3) 


are  necessary.   Hence 


tR'-Jt)      1 


R,(t)    2R9(s) 


R,2(s)    /    l\    ,  , 
— - +  s  +  -R'2(s) 


(20-3) 


Q 
The  two  terms  of  this  equation  have  a  constant  value  -  — .   There- 
in 
fore,   Ro  satisfies  the  linear  differential  equation 


sR 


"2  +  2(s2  +  l)R'2  +  QsR2 


(20-4) 


or,  reduced  to  classical  form,  by  the  transformation, 


sR2  =  e  S  H 


(20-5) 


Quantity  H  confirms  in  fact  the  Hermitian  differential  equation 


H"  -  2sH'  +  (Q  -  k)H   =  0 


(20-6) 


78  NACA  TM  1377 


Since  R(r,t)  must  cancel  out  for  r  infinite,  hence  also  H(s) 
for  s  infinite,  Q  -  k  is  necessarily  an  even  integer  2n.  In  that 
case  H  is  the  nth  Hermitian  polynomial : 


Hn(s)  =  es2-^e"s2  (20-7) 

as 


and 


R2(s)  =  ±e-s\(s)  (20-8) 

S 


But  Rg  must  have  a  finite  limit  when  s — X),  which  requires  that 
H  (s)  "be  an  odd  polynomial,  hence  that  n  be  an  odd  integer:   n  =  2p  -  1 
Q  =  4p  +  2.   The  function  Rj_(t)   is  then  equal  to  — = — .   Hence  the  ele- 

mentary  solution: 

r2 


*.*>  -  £  ^P-fe)  (--9) 


The  general  solution  is  a  superposition  of  elementary  solutions : 

r2 

P=l  tP     \\/8vt/ 

the  constants   A   being  so  chosen  that  the  series  converges  and  R(0,t) 
is  equal  to   3uq2 • 

Examples:   1.  Limited  to  the  term  in  p  =  2  it  is  found  that,  by 
interpreting  the  Hermitian  polynomial,   Hj(s),  and  denoting  an  inde- 
terminate constant  by  A: 
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R(r,t)  =  A3  _  £l_e  0vt  (20-11) 


hence 


-A  5 

f(r,t)  =  e  8vt     U(,2(t)  =  At  2 


This  example  is  rational.   It  is  compatible  with  the  hypotheses 

made  at  that  time.   In  particular,   r^  approaches  zero  when  r — ><», 
whatever  the  non-negative  number  m  may  be.   However,   f  =  0  for 
t  =  0,  except  when  r  =  0.   At  the  initial  instant  the  turbulence  is 
concentrated  at  one  point  from  which  it  ultimately  spreads  throughout 
the  entire  fluid.   Such  a  turbulent  structure  seems  rather  difficult 
to  conceive  when  assuming  it  to  start  from  t  =  0.   An  interpretation 
will  be  given  later  (section  ~^k)   in  connection  with  the  final  phase  of 
decay  of  turbulence  behind  a  grid. 

For  p  =  1,  one  would  have 


r2 

r'2 

A      8vt 

f  = 

.3      rr.2e'8vtar 

"I6 

r5  Jo 

t2 

R(r,t)  =Aeuvl     f  =  2-    \     r^e   Ul/bdr'        (20-12) 


f  would  not  converge  toward  zero  when  r — ><»,  which  proves  that  this 
example  cannot  be  suitable  for  a  real  motion. 

2.  (Karman-Howarth  (ref.  30).   Suppose  the  correlation  function  f 
does  not  depend  separately  on  r  and  t  but  only  on  the  vari- 
able -  -   -£- 


u02(3f+r|i 
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it   then  yields 

R  =  uQ2(t)[3f(s)   +   sf(sj]  (20-13) 


so  that  R  is  the  product  of  a  function  of  t  "by  a  function  of  s. 
Thus  R  is  an  elementary  solution  of  the  problem  treated  in  the  first 
example,  corresponding  to  a  chosen  Hermitian  function  Hp  i- 

The  function  f(s)   satisfies  a  differential  equation  similar  to 
(20-4)  written  by  Karman  and  Howarth.   But  a  great  advantage  accrues 
from  the  use  of  function  R  and  reduction  to  Hermitian  polynomials . 

Remark.   According  to  (20-10),   R(r,t)   is  a  sum  of  "elementary 
solutions"  corresponding  to  various  values  of  the  integer  p,  starting 
from  p  =  2.   For  each  one  of  these  solutions,  except  for  p  =  2,  it 
is  verified  that  the  invariant  of  Loitsiansky  (18-10)  is  zero.   For 
p  =  2  it  has  a  finite  value . 

21.  Solutions  involving  a  similarity  hypothesis  (see  ref.  7): 

If  the  triple  correlations  are  no  longer  negligible,  the  fundamental 
equation  cannot  be  solved  completely.  It  is  assumed  that  the  correlation 
functions  f(r,t)   and  c(r,t)   are  not  separately  dependent  on  the  tvo 

variables  r,t,  but  solely  on  i|/  =  — f —  where   Z(t)   is  a  length  in 

Z  (t ) 

terms  of  time.   That  is  to  say,  that  at  each  instant  t  the  correlation 

curves  are  superposable ,  by  means  of  a  unit  change  on  the  axis  of  r. 

However,  it  is  recognized  that  this  hypothesis,  called  "total  similarity" 

is  often  too  specific,  and  it  is  therefore  replaced  by  a  "partial 

similarity"  which  is  verified  only  in  a  finite  interval   Z^<  r  <  Zg. 

For  total  similarity,  the  Loitsiansky  invariant  reads 


r^f(r,t)dr  =  uQ2(t)Z5(t)    i|Af(*)di|r       (21-1) 

0  Jo 

Consequently,   Uq  (t)Z-)(t)   is  a  constant  independent  of  t  in  the 

course  of  the  motion.   This  is  no  longer  true  in  the  case  of  partial 
similarity. 
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To  transform  the  fundamental  equation,  it  is  possible  to  put 

2/ 


R(r,t)  =  3V  (t)a(i|r) 


T(r,t)  =  3uQ3(t)p(t)  > 


*l   = 


UqZ 


(21-2) 


a(t)   is  a  function  that  takes  the  value  unity  for  \|r  =  0;  (3(\(f)   is 
a  function  that  cancels  out  with  i|r .  R7   is  a  Reynolds  number  in  terms 

of  the  time,  and  T  is  the  scalar  which,  according  to  (l8-l4),  defines 
the  triple  correlations . 

Considering  the  relation  (19-1)  which  affords  — Uq2  in  terms  of 

o 

Uq   and  of  the  length  A,  the  fundamental  equati 


on 


dR  -* 

—  =  2v  &R  +  div  T 

dt 


is  easily  transformed  to 


-l_(z2Wa'(t)   +  5^(*)   +  a"(t)   +  ia'U)  +  ±R, 
4V  v      '  ^2  t  2— ' 


P'OlO   +  tP(+) 


=  0 


(21-3) 


or,  schematically  arranged: 


a1(t)31(jf)  +  a2(t)p2(\|f)  +  a3(t)35(v|r)  +  aj+(t)p1+(t)  =  0    (21-4) 


Equation  (21-3)  is  a  sum  of  four  terms  each  of  which  is  the  product 
of  a  function  of  t  by  a  function  of  \|r .   How  can  an  equation  such  as 
(21-4)  be  resolved? 
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The  a^(t)   are  considered  as  the  coordinates  of  a  point  A  in 

four -dimensional  space,  this  point  describing  a  curve  (a)  in  terms  of 
the  parameter  t.   The  same  applies  to  the  (3-j_(v).   It  results  in  two 

curves  (fig.  6)  such  that  the  straight  lines  joining  the  origin  to  any 
point  A  of  the  first  curve  and  to  any  point  B  of  the  second  are 
always  perpendicular.  These  curves  must  fit  two  orthogonal  complemen- 
tary subspaces  in  four-dimensional  space.   Hence  there  are,  a  priori, 
three  possible  cases: 

1.  (A)  is  a  line  passing  through  the  origin.   (B)  is  in  the  com- 
plementary three-dimensional  space  orthogonal  to  this  line. 

2.  (A)  and  (b)  are  in  two  completely  orthogonal  planes  passing 
through  the  origin. 

3-  (A)  is  in  a  three-dimensional  linear  space.   (B)  is  the  orthog- 
onal line . 

In  the  first  case,  the  a^(t)  are  proportional  to  constant  num- 
bers nij_ . 

In  returning  to  the  notations  of  equation  (21-3),  the  quantities 

O  72 

are  proportional  to  constant  numbers .   In  other  words ,   Z ,   — ,  and  FL 

A2       ~l 

are  constant.   Following  this,  the  functions  <x(ii)     and  p(t)  are 

connected  by  a  unique  differential  equation.   Obviously  one  may  assume 

o 
Z  =  A.   Then  A   is  a  linear  function  of  time  and  UqA  is  constant. 

Considering  equation  (19-1),  which  we  recall  here 

Au  2  =  _10vVl  (19-1) 

dt  0        a2 
we  obtain  of  necessity: 


2  =  A2  =  lOvt     uQ2  =  5-f.     RA  =  C^         (21-5) 
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The  Loitsiansky  relation  is  not  verified,  so  that  no  total  simi- 
larity can  prevail;  a(t)  and  p(\|/)   are  joined  by  the  relation 


a"  +  (%  +  2  la'  +  5a  +  -%  (p '  +  ^A   =  0  (21-6) 

12      *r  2  V   V  j 


which  can  be  subjected  to  experimental  check. 

A  complete  discussion  of  the  second  and  third  possibility  is  fore- 
gone, in  favor  of  the  case  where  I   =   A  is  prescribed.   Then  (21-3) 
becomes 


i(A2)Va'(t)  +  an(lr)  +  ^a'(t)  +  5a(t)  +  fe> 


and  is  represented  schematically  by 


(21-7) 


a1(t)p1(|0  +  a2(t)p2(^)  +  a5(t)p3('^)  =  0 


(21-7) 


The  discussion  then  deals  with  three-dimensional  rather  than  four- 
dimensional  space. 

Two  cases  are  possible: 

(1)  The  point  a^(t)   describes  a  straight  line  (issuing  from  0), 
and  point  P^OlO  remains  in  the  plane  perpendicular  to  this  line 
(passing  through  0). 

(2)  The  roles  of  the  two  points  are  permutated. 

The  first  case  is  not  distinct  from  that  which  we  have  studied. 
In  the  second  case,  there  exist  three  constants  m-,  ,m?,m,,  so  that 
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2„, 


a,^)        cc"(->)  +  |a'(t)   +  5a(>)        P'(^)   +  |p(+) 


-vku 


mc 


2m, 


ml(  ^  )     +  m2  +  m3^\  =  ° 


(21-8) 


They  axe  the  conditions  proposed  by  Sedov.   They  can  be  checked 
directly  by  experiment;  m^  cannot  be  zero,  because  K-,   would  be  =  C   , 

and  at  the  same  time  aOlO  =  C   ,  which  is  not  possible;  nh  and  m, 
cannot  be  zero  simultaneously. 

If  m-z  ^  0,  the  similarity  cannot  be  total,  because  one  would  have 

Uq2A5  =  Cte   and  by  (19-1) 


A^  =  kvt 


"6 


2t5/2  =  cte     ^t3A  =  cte 


which  is  incompatible  with  the  second  equation  (21-8). 

Following  this,   Uq  and  A  are  solutions  of  the  system  of  differ- 
ential equations : 


m1(A 


:j'  +m2 


+   — unA  =  0 
v  0 


/   2>        l 


(21-9) 


a(r)      verifies   the   second-order  differential  equation: 


(2.         '"P 

a"  +   !  -  -  r-^if)a.'   +  5a  =  0 
\\|r        ^vm^ 


(21-10) 
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after  which  3  is  determined  by  the  first -order  differential  equation 


p»  +  8p  =  — 2_^a«  (21-11) 

t    2vmj_ 


The  system  (21-9)  can  be  reduced  to  quadratures.   If  the  solution 
cannot  be  written  in  finite  terms,  it  still  is  desirable  to  integrate 
as  far  as  possible. 

When  the  first  equation  (21-9)  is  set  in  the  form 


J_=!^JS(A2V  +k  (2l.l2) 

Rq    IOv  \   / 


where  Rp,  and  k  are  constants  replacing  m-j_,  mo,  and  nw,  it  is 
found  that 


A  =  C^Wr.  -  Ro)k  (21-13) 


Whence  it  is  deduced  that 


m  =  cte   *-%    L  _  5oY"; 

dt  R  -   (l  -  kjRjrjl  R  J 


=  cte = =0 -l    _  =0  (21-U) 


The  problem  is  thus  reduced  to  a  quadrature.   By  (21-13),   A 
follows  from  R,  and  Uq  =  v=. 

The  equation  (21-10)  generalizes  equation  (20-4)  to  a  certain  extent, 
but  it  is  less  simple  and  will  not  be  discussed. 

Other  types  of  solutions  satisfying  a  similarity  hypothesis  can  be 
examined,  such  as  the  case  where  A  is  very  small,  that  is,  where  the 
dissipation  due  to  viscosity  is  low.   In  that  case  a  length  other  than 

Poo 

A  must  be  chosen  for  length  Z,  such  as  1=1  f(r)dr.  The  discussion 

1  0 
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proceeds  on  the  assumption  that  the  viscosity  terms  of  the  fundamental 
equation  are  negligible  against  the  other  two.   It  leaves  then 


^.=  div  T  (21-15) 

dt 


The  calculations  are  not  developed  (for  further  details,  consult 
Batchelor's  report  (ref.  7))-   The  sole  purpose  in  this  paragraph  was  to 
give  an  idea  of  the  various  methods  that  can  be  applied  to  solve  the 
fundamental  equation  on  the  basis  of  the  similarity  hypothesis.   The 
physical  study  of  the  problem  is  deferred  until  later. 


22 .   Transformation  of  the  fundamental  equation  in  spectral  terms : 

The  formulas  of  reciprocity  between  the  correlations  and  the  spec- 
tral functions  made  it  seem  interesting  to  transcribe  the  fundamental 
equation  (17-10)  in  spectral  terms.   (Compare  spectral  intermedium, 
section  19- ) 

By  (11-10): 


R 


with 


<P(W*3)  =  H^aa 


Likewise 


3 


1  A  Vp 

AR  =  -/  e  p_±     k2cpdA  (22-2) 
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Hence 


^  -  2v  AR  = 


P1  2_  Vpk  \ 

e     P=1  p+  2vk2cp  dA  (22-3) 


— >  -> 

This   quantity  is  equal  to     div  T,   where     T     is   the   vector  having 


for  components      /      Tn    .     We  put 


a 


1  h   Aplp 
div  T         /  ^f(A)e     p  dA 


The  fundamental  equation  becomes 


-^  +  2vk2cv  =  t  (22-1+) 

St 


In  this  equation  cp  and  \|r  are  unknown  functions  of  A,   t,  and 


k2  =  ^AP2. 


P=l 

The  total  energy  of  the  fluctuation  is 


E  .  i  1^ 


(it  was   shown  in  section  9  that       /  i[rdA  =  0) 

,/A 


For  isotropic  turbulences,  the  spectral  function  F(k)  =  2nk2cp, 
which  does  not  depend  separately  on  A^,A2>A^  as  in  the  general  case 

for  cp,  but  solely  on  k,  is  introduced. 
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The  factor  >|/   can  also  be  expressed  by  means  of  the  scalar 

T  =  u03(c  +  2a)  =  ^r(c"  +  M  =  ^!fcr^' 

<-   V      /    Pr 


It  is  actually  known  that 

div  T  =  T'  +  -T  =  —  (t2tY  (22-5) 

r    r2  V    ' 


which,  after  elementary  calculations,  gives 


.j,(k)  =  _1_     sin  rk  ^r2T\ '  ^  (22_6) 

2rt2  Jo     rk  V 

Equation  (22-4)  reads 

|£  +  2vk2F  =  2«k2t  (22-7) 

at 


Putting:   2itk2i>  =  T'(k,t),  where  $   is  a  function  connected  with 
triple  correlations  and  such  that 


$   dk  =  0 
0 


yields 


—  +  2vk2F  =  ^  (22-8) 

at 

By  the  use  of  the  formula  (22-5)  and  the  limited  development  of 
sin  rk,  it  is  readily  proved  that  the  function  $(k)   is  infinitely 

small  as  k  ,  provided  that  the  correlation  function  c  approaches 
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zero  faster  at  infinity  than  — .   But,  with  the  relation   J  |  dk  =  0, 

r^  JO 

nothing  more  is  known  a  priori  about  its  form. 

Equation  (22-8)  is  not  essentially  different  from  (17-10).   It 
raises  the  same  difficulties  and  is  always  a  single  equation  between 
two  unknowns  F  and  J.   However,  it  is  much  easier  to  make  mathemati- 
cally or  physically  reasonable  assumptions  about  the  form  of  $  than 
about  the  form  of  the  triple  correlations,  as  will  be  shown  in  the  fol- 
lowing two  methods,  both  due  to  Heisenberg,  which  take  advantage  of 
(22-8);  the  first  involves  the  introduction  of  the  correlations  in  time 
and  will  go  beyond  the  purpose  originally  assigned  to  it,  the  second 
leads  effectively  to  certain  possible  forms  of  the  spectral  function  F. 


23.  First  theory  of  Heisenberg  (ref.  23): 

The  hypotheses  to  be  examined  primarily  are  those  which  form  the 
subject  of  the  last  part  of  Heisenberg 's  report.  We  shall  utilize  the 
ideas  and  follow  the  calculations  as  closely  as  possible,  but  change 

the  notations  and  substitute  stochastic  averages1^  for  the  spatial  and 
temporal  averages. 

The  expressions  giving  the  statistical  velocity  and  pressure  of 
homogeneous  turbulence  read: 

3  3 


p   i  >_  Apxp  i  X  ApXp 

ua  =  I  e  P=1     dha(A)      p  =  /  e  P=1     dq(A)    (23-1) 


JA 


The  equations  of  motion  are  first  written  in  spectral  terms.   It 
can  be  assumed  that  h a(A)   is  derivable  on  each  test  with  respect  to 

the  time,  and  its  derivative,  which  corresponds  to  — — ,  is  designated 


at 


by  ha(A). 


10Bass,  reference  5.  The  calculations  are  rather  difficult.   The 
reader  who  wants  to  avoid  reading  them  will  find  a  summary  of  the 
results  at  the  end  of  paragraph  2h. 
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The  expressions  (23-1)  must  be  introduced  in  the  equations  (23-2) 


at 


(X 


Z_ub:— 


o"u„    du 


+  U- 


a 


bXn  Sx 


1  op      y  ^2ua 

+  v  2_  — - 


1     P  dxa 


P   ox,2 


Sxa 


=  0 


(23-2) 


where  U  is  the  constant  over-all  velocity. 

The  only  term  of  which  the  equation  merits  a  more  detailed  dis- 
cussion is 


-i     A'pe 
,/AXA 


A' 


p*p  1 

e 


^pxp 


dha^(A')dhp(u) 


(23-3) 

Putting  p.  -  A'p  =  ?Vp  and  considering  the  equation  of  continuity 
yields 


Y~A  dh  (A)  =  0 


which  had  already  been  used  in  the  form 


01%       / 
^P     P  JW 


^pXp 


dV^n  -  A)dhp(p)      (23-M 


the  subscript  (u)  specifying  that  the  integral  is  extended  over  the  values 
of  p..   This  quantity  is  the  Fourier  transform  with  respect  to  A  of 


iAp     dha*(u  -  A)dhp(u) 

J(m-) 
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Ey  means  of  some  elementary  transformations  two  equations  are 
obtained: 


dh 


V   PJ(u)   a  P        1    a 


-JAa  dq(A)  -  vk2dha(A) 


>    (23-5) 


^=4rvp  / *v^  -  a)^  (u) 

p      k^  ap      ,/p.  p 


where  h  (A)   designates  the  derivative  of  h  (A)  with  respect  to  A 


a 


av 


The  demonstration  of  paragraph  17  from  equations  (23-5)  is  then 
resumed.   The  first  equation  (23-5)  is  multiplied  by  dha*(A),  added  up 

with  respect  to  a,   i   is  then  changed  to   -i   and  the  two  equations 
obtained  are  added  member  by  member. 


Patting 


Z  =  -2i  Yl   Ap  <%*(*)  f     dha*(u  -  A)dhp(u)        (23-6) 
a3  i/(n) 


while  considering  the  equation  of  continuity,  gives 


3_ 


dha(A)  2  +  2vk2  ^>~ 


dh  (A) 


=  R(Z)         (23-7) 


R(Z)   denotes  the  real  part  of  Z. 

After  averaging,  it  is  known  that  when  the  turbulence  is  isotropic 


|dh  (A)  I2  =  -1—   dA 
■   a,  o 

1    2itk2 
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Consequently 


2rtk2\?t  /         " 


What  was  called  $  in  paragraph  22  is  given  by  the  equation 


2«k2R(Z)  =  t  dA  (23-9) 


In  order  to  find  a  useful  expression  of  Z  several  hypotheses 
originating  with  Heisenberg  are  made . 

First  of  all,  it  is  known  that  the  statistical  functions  dha(A) 

are  of  orthogonal  increments,  that  is,  that  the  averages  of  the  pro- 
ducts dh  (A)dh  *(A' )  are  zero,  if  the  two  points  A, A'   of  the  wave 

number  space  are  distinct.   The  first  hypothesis  is  that  the  statistical 
variables  dha(A),   dha*(7\')   are  not  only  orthogonal  but  independent 

(in  probability). 


Recalling  that  the  average  values  dh  (A)   are  zero  and  that 


a" 


dha*(A)  =  dha(-A),  it  is  seen  that  the  average  value  of 

clha*(X)dV(n  -  A)dhp(^) 
is  zero,  except  when  at  the  same  time 

A  =  A  -  u     u  -  A  =  p     A  =  p 


that  is,   A  =  u  =  0.   If  it  is  assumed  that  the  probabilities  of  sudden 
jumps  of  h  (A)   are  zero,  then  Z  =  0. 

Paragraph  20  dealt  with  the  case  where  Z  =  0.   It  is  certain  that, 
in  general,   Z  is  not  rigorously  zero,  but  it  is  likely  that  Z  is  small 

compared  to  the  viscosity  term,  or  the  term  of  development  — .   In  other 

ot 
words,  the  hypothesis  of  orthogonality  is  too  general,  and  that  of  the 
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independence  too  restrictive,  yet  it  is  difficult  to  formulate  an  inter- 
mediary hypothesis.   But,  starting  from  the  independence  it  is  possible 
to  obtain,  by  means  of  certain  interesting  transformations,  a  significant 
expression  for  Z.   This  expression  is  defective  since  there  is  contra- 
diction.  But  it  seems  likely  that  it  represents  a  good  approximation 
of  Z.  We  shall  therefore  study  it. 

2k.   First  theory  of  Heisenberg  (continued).   Space-time  correlations: 

If  the  average  of  R(Z)   disappears,  it  is  because  Z  contains  an 
odd  number  of  factors  dh  (A);  the  equations  of  motion  are  not  linear. 

Replace  dhn(u)  by  the  time  integral  of  its  derivative  dhR(u)  and 
express  dh-,(u)  by  means  of  the  equations  of  motion.   Certain  terms 

linear  in  dh  are  replaced  by  quadratic  terms;  products  of  four  fac- 
tors dh  appear  in  Z,  and  Z  is  not  zero.   A  kind  of  technique  of 
solving  the  equations  of  motion  is  involved  here,  which  reduces  them 
to  integral  equations  (bearing  on  simple  averages). 

If  T  is  a  very  great  positive  number 

Pt  Pt+T 

dhp(u,t)  -  dhp(n,-T)  =  /   dhp(u,T)dT  =  I    dhp(u,t  -  T)dT 


(24-1) 


If  T  is  sufficiently  large  it  may  be  conceded  that  hn(u,-T) 

tends  toward  zero,  that  is,  that  the  fluid  has  "started  from  rest." 
Hence 


dhp(n,t)  =    dhp(n,t  -  x)dr  (24-2) 


To  lighten  the  notations,   h'g(u)   is  to  represent  the  func- 
tion hg(p,t  -  t).   Then 


Z  =  -2i  ,    AB  dh  *(A)    dr     dha*(p.  -  A)dh'  (u)    (24-3) 
a3  Jo    J(u) 
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and  dh'R(u)   is  eliminated  by  means  of  the  equations  of  motion  from 
which  we  removed  previously  the  term  of  the  pressure  in  dq: 


^L 


dh'p(ti)  =  -i  >   txy  jf        dh'R*(u'  -u)dh'7(u')  + 


i^|^uaVp'  /    dh,a,.*(^'  -  n)cUi'p.(^)  -  (vu2  +  imU)dh'p(ti) 


(2k-k) 


with  p.2  =  ^_  ^ 


Now  a  new  assumption  is  made  according  to  which  the  original  assump- 
tion of  independence  (for  t  =  0),  which  is  certainly  verified  for  the 
great  values  of  r,    remains  valid  for  the  whole  time  interval  t ,    leaving 
only  the  averages  with  four  factors 

r 


Z  =  -2Z-H         dT  ti7dha*(A)dJia*(H  -  TOdh'p*^'    -  vi)dh'    (n')   + 

agr         uo         J(n,n') 


y-       V^aVV     pdT      P  dh*(A)dh*(^   .   X)dh-     ,*(n'    -   n)dh'B,(n')  (24-5) 


2 

aPa'p'  ^2  Jo  J(ii,n') 

The  dh  being  independent  in  distinct  points  of  the  frequency  space, | 
it  is  necessary  to  associate  the  four  points  which  correspond  to  the 
four  dh  by  pairs,  not  forgetting  that  dh  *(A)  =  dh_(-A).   Considering 

the  transformation  of  the  first  integral,  the  only  admissible  combinations 
are 

(l)st 

A  =  A  -  \i  u'  —  M-  =  H1 

hence  p.  =  0 .   No  contribution  to  the  integral . 
(2)nd 

A  =  u  -  u'      u-  -  A  =  n' 
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conditions  which  are  reduced  to  \i'    =  \x  -  ~h,   and  leave  an  integral  other 
than  zero 

/   u7dh^*(A)dh'p(A)  dha*(u  -  A)dh'7(u  -  A)        (24-6) 

J(n) 

(3)rd 

A=u'      u-A=u-u' 

conditions  which  are  reduced  to     ^'    =  A,   and  leave  "behind  the   integral 

/      u7dha*(A)dh'7(A)    dha*(u   -  A)dh'p(n   -  A)  (24-7) 

J(^) 


After  some  changes  in  the  denomination  of  the  subscripts,  the 
first  term  of  Z  assumes  the  form 


n 


-2/>      dr    (Apu7  +  A7up)dha*(A)dh'p(A)  dha*(u  -  A)dh'7(u  -  A) 
a(37  JO    Ku) 

(24-8) 

Putting  u  -  A  =  A',  while  considering  the  incompressibility  condi- 
tion, equation  (24-8)  becomes 


-2^_  A7A'p    dr      dha*(A)dh'p(A)  dha*(A')dh«7(A') 

aP7       JO    J(V)  (^_g) 

A  similar  calculation  which  need  not  be  given  in  detail  gives  the 
expression  of  the  second  term  of  Z.  Finally 

Z  =  2    dT  /    ^ — ^-^ 


0    J(A')  o     ,0     ^~ 

;k2  +  k'2  +  2  2_  ARA 


6A  5 


^A7A'3dha*(A)dh'p(A)  dha*(A1)dh'7(AI)         (24-10) 
a37 
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Now  the  components  of  the  correlation  tensor  in  space  and  time  are 
introduced: 


Sap(s,t,T)  =  ua(x,t)u3(x  +  i,t   -  t)  (24-11) 


The  ua  are  expressed  in  spectral  terms.  Putting 


dh^*(70dh'B(A)  =  dha*(A,t)dhB(A,t  -  t)  =  6a3(A,t,T)dA 

(24-12) 


it  is  seen  that 


y  ap^i 


Sop(l,t,T)  =  l^2—    P5P  eap(A,t,T)dA        (24-13) 

a  formula  which  ties  the  correlations  SaB  to  the  generalized  spectral 
tensor  9ag.  Now,   Z  is  expressed  by  means  of  0aR: 


Z  =  2  dA    dT 


k2  -  k'2 


10    J(A')k2  +  k,2  +  2^--  AsA, 


ZZ.V,3eap(^t'T)ea7^,'t'T)d^'  (24-14) 

aP7 


For  easier  reading  of  this  formula,  the  term  under  the  sign   /  represents 
the  scalar  product  of  the  two  vectors  having  for  components 

Z_  A'peap^'t'T)    and    y~Ay9ay(A',t,T) 


The  real  part  of  Z  is  obtained  by  transforming  the   9  3   into 
correlations  Sa0  by  (24-13)  and  taking  the  real  parts  of  the 
exponentials . 
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All  these  calculations  do  not  necessarily  assume  isotropic  turbulence, 
But  in  this  particular  case  they  are  a  little  simplified  and  will  be  termi- 
nated.  If  there  is  isotropy  and  incompressibility: 


\,AC 


3aj3 


(A,t,r)  =  A(k,t,T)  6^0  -  -SLfi 


a(3 


(24-15) 


A  =  —  2_ 9aa.  being  a  real  quantity  to  be  interpreted  right  away. 


A  being  real,   Z  is  real  and  equal  to  R(z).  Therefore 


R(Z)  =  2  dA    dT 


„  r 


,2  _  ,..2 


J(A') 


k2  +  k'2  +  2 


I 


¥'6 


v 


apr 


_  A7A'pA(k,t,T)A(k',t,T)(-HLl  _  S/vQ)[   a  y  . 


Ja3j  v   .2 
k 


JaB 


(24-16) 

To  terminate  the  calculation,  the  last  factors  are  developed,  and 
with  the  vector  A  taken  as  support  of  the  third  axis  one  passes  to 
polar  coordinates.   The  angle  6  of  vectors  A  and  A'   is  introduced, 
defined  by 


A&A'5  =  kk'  cos  9 


Finally,  one  puts  cos  9  =  £  and  introduces  the  elementary  function 


C(l   -   C2) 


J(k,k')   =  -4n(k2   -  k'2)k3k'3    I      ^X  " 

'  J-1  k2  +  k-2  + 

=  *(k2    -  k'2Vk^  +  k'^    -  SkV2) 


^^ 


2kk't; 


i  (^ 


+  k-2Vk2 

kk' 


log 


k  +  k' 
k   -  k' 


(24-17) 
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Summary  and  conclusion.-  By  the  successive  assumptions  of  the  pre- 
sent paragraph,  the  fundamental  equation  (22-8)  reads 


..  loo    noo 

|£  +  2vk2F  =    dT    j(k,k')A(k,t,T)A(k',t,T)dk'     (24-18) 

o*  UO    Jo 

The  function  j(k,k')  given  by  (24-17)  changes  sign  when  k  and 
k'   are  permuted.   The  integral  of  the  second  member  of  (24-18)  with 
respect  to  k,  from  0  to  oo^  is  therefore  zero  and  (24-18)  contains, 
as  it  should  be,  the  known  equation  (19-2): 

—  +  e  =  0 

St 

as  it  follows  from  formulas  (13-4)  and  (13-IO). 

The  function  A(k,t,f)   is  tied  to  the  space-time  correlations: 


sa,p(£>t>T)  =  ua(x^t)u(3(x  +  i^t  -  T) 
by  the  relation 


Yi  8^(1, t,r)  =  2  /V  51  AP^P  A(k,t,r)dA 


If  the  turbulence  is  isotropic  it  gives,  more  simply: 


51SaaU>t,T)    =   8*     j     k2   sir^rkA(kjt;T)dic  (24_19) 


JO      rk 


and  conversely 


A(k^>^l/2^(^S>       ^^ 


SQja  being  naturally  a  function  of  r  alone. 
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If  t  =  0,  the  general  tensor  S  n(|,t,T)   is  reduced  to  the  COrre- 
a0' 


lation  tensor  in  space  RaR(|,t)   and  y S    becomes 


ZRoa  =  R' 


F(k,t)  =  4jtk2A(k,t,0) 


(24-21) 


is  verified  whereby  the  formulas  (24-19)  and  (24-20)  contain  the  formu- 
las (l4-4)  and  (14-5).   Therefore,  (24-18)  is  a  nonlinear  integral  equa- 
tion, in  which  appear  at  the  same  time  A(k,t,-r)   and  A(k,t,0)  =  F(k,t) 
and  which  makes  it  possible  in  principle  to  determine  the  general 
correlations,  in  time  and  space,  of  isotropic  turbulence  by  means  of 
their  spectral  function  A(k,t,-r). 


25.  Second  theory  of  Heisenberg: 

This  "second  theory"  came  first  in  Heisenberg's  report  (ref.  23). 
By  means  of  postulates  of  physical  character,  it  effectively  enables  the 
form  of  the  spectral  function  F(k)   to  be  determined  in  certain  conditions, 
But  the  formula  set  up  still  contains  an  indeterminate  constant.   It  seems 
that  Heisenberg's  idea  was  to  carry  the  expression  obtained  for  F  in  the 
nonlinear  equations  deduced  from  the  equations  of  motion,  which  would  have 
made  it  possible  to  calculate  the  constant.   The  theory  of  sections  23 
and  24  does  not  exactly  come  up  to  these  expectations.   It  results  in  a 
nonlinear  integral  equation  in  which  not  the  spectral  function  F(k,t)  but 
its  space-time  extension  A(k,t,-r)   is  involved.   If  A(k,t,0)  =  F(k,t) 
are  known,  by  the  theory  which  is  to  be  discussed,  it  seems  theoretically 
possible  to  deduce  from  it  A(k,t,i")  for  every  time  interval  t.   But 
the  mathematical  problem  thus  posed  appears  difficult  and  does  not  fit 
into  any  class  of  known  problems,  and  no  attempt  is  made  to  solve  it. 

The  dissipation  of  energy  per  unit  mass  in  the  ensemble  of  the  spec- 
trum is ,  as  is  known 


e  =  2v 


k2F(k)dk 


(25-1) 
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The  dissipated  energy  per  unit  mass  in  the  spectral  range 

\- —   2    2  ~i  " 

J> A  <   k  ,  that  is,  by  the  eddies  of  "dimensions  greater  than  — 

and  of  relatively  large  scale  are  designated  by  ek.   If  k — =><»^ 
ek  tends  toward  e,  and  if  k  >0,   ek  tends  toward  zero.   Heisenberg 
reduced  ek  to  a  form  similar  to  e  by  introducing  an  apparent  coef- 
ficient of  viscosity  v^  so  that 

ek  =  2^v  +  vk)  /  k,2F(k')dk'  (25-2) 

vk  must  tend  toward  zero  when  k — )°°. 

Now  this  first  hypothesis  must  be  perfected  in  such  a  way  that   v^ 
assumes  a  usable  form. 


Since   lim  v^  =  0,   vk  can  be  represented  by  an  integral   / 

poo  Jk 

What  about  the  function  under  the  sign  1      ?  It  is  assumed  that  this 

Jk 
function  depends  on  the  values  k'   ranging  between  k  and  <»  directly 
and  by  means  of  F(k'),  but  that  it  does  not  depend  on  k.   In  that  case, 
it  can  be  determined  by  dimensional  consideration.   Its  product  by  dk' 
having  the  dimensions  of  a  viscosity  v,  it  is  necessarily  of  the  form 

C  / — ,  where  C   is  a  purely  numerical  constant. 

V  k-5 

Hence 


ek  =  2(v  +  C  /   /liii-1  dk")   /  k'2F(k')dk'        (25-3) 


ek  is  the  sum  of  two  terms.   The  first  corresponds  to  the  part  of 
the  energy  of  the  fluctuations  of  the  spectral  range  (0,k)  which  is 
transformed  directly  into  heat.   The  second  is  the  part  of  that  energy 
which,  before  being  transformed  into  heat,  serves,  first,  to  maintain 
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the  energy  of  the  fluctuations  of  the  spectral  range  (k,oo)  of  the  "small 
eddies ." 

Now,  the  fundamental  equation  (22-8)  can  be  decomposed  precisely  in 
this  fashion.   Integration  member  by  member  from  0  to  k  yields 


—  /  F(k')dk'  +  2v    k'2F(k')dk'  =    i]/(k')dk'      (25-4) 

dt  Jo  Jo  Jo 


The  first  member  manifests,  first,  the  energy  loss  of  the  fluctua- 
tions of  the  range  (0,k),  then  the  portion  of  this  energy  dissipated 
directly  by  viscosity.  The  second  member  is  therefore  the  portion  which 
is  transformed  into  energy  of  the  fluctuations  of  the  range  (k,oo),  and 
the  formula  (25-3)  gives  an  expression  of  the  function  f.      Ultimately 


|-    F(k')dk'  +  2[v  +  C     M^H   cLki]     k'2F(k')dk'  =  0 

^  Jo  \     Jkjk'3    jjo 

(25-5) 

This  is  Heisenberg's  fundamental  equation-*-!.  For  the  present,  it 
is  not  discussed  in  its  most  general  form,  but  only  in  a  special  case. 

Heisenberg's  hypotheses  are  incompatible  with  the  concept  of  abso- 
lute steadiness,  because  (25-5)  has  no  solution  independent  of  t.   In 
chapter  V,  it  will  be  shown  that  Heisenberg's  equation  can  be  applied 
to  the  problem  of  spontaneous  decay  of  turbulence  behind  a  grid,  an 
essentially  unsteady  problem.   But,  an  assumption  of  partial  steadiness 
can  be  made  as  will  be  shown  in  chapter  IV.  For  the  values  of  k 
greater  than  a  fixed  approximate  value  kQ,  it  can  be  assumed  that  the 

energy  dissipated  by  the  large  eddies  (wave  numbers  lower  than  k)  is  a 
constant.   Only  the  manner  of  distribution  between  the  turbulent  fluc- 
tuations of  the  small  eddies  (wave  numbers  higher  than  k)  and  the 

It  appears  in  this  form  in  the  second  of  Heisenberg's  reports 
(reference  2k    ),  but  the  particular  case  in  question  is  treated  in  his 
first  report  (reference  23  ) . 
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dissipation  by  viscosity  (transformation  into  heat,  molecular  fluctua- 
tions) depend  en  k,  the  proportion  dissipated  as  heat  increasing  with 
k  according  to  formula  (25-3). 

It  will  be  shown  (in  chapter  IV)  that  this  assumption  of  partial 
steadiness  is  equivalent  to  an  assumption  of  statistical  equilibrium  or 
similarity,  which  is  the  more  accurate  the  higher  the  Reynolds  number  of 
the  turbulence.   On  these  premises  the  hypothesis  is  discussed  from  the 
mathematical  point  of  view.   Equation  (25-5)  must  be  replaced  by 


2   v  +  C 


feli-  dkj     J     kl2F(kf)dk'   ==  Cte  =   e0         if     k  >  Isq     (25-6) 


To  solve   it,    simply  put 


(k)   :     /     k,2F(k')dk'  F(k)   =  ie'(k) 


(25-7) 


The  function  6(k)   cancels  out  for  k  =  0  and  verifies  the  equation 


v  +  C 


fO'(k') 


,5 


dk' 


£0 


29(k) 


(25-8) 


which  by  differentiation  becomes  the  elementary  differential  equation 

,2 


Qk  "   \eoJ  k5 


(25-9) 


*,  h   eo 
Putting  RJ  -  —  —  and  introducing  an  integration  constant  kj_  it 

is  seen  that 


(k)  =  Hk5 


L 


k 

3 


(25-10) 
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and  finally 


F(k)  = 


.2 

\3 

;  *  (tfl 

i 

r(^ 

k 
'3 


(25-11) 


This  formula  defines  the  spectrum  of  isotropic  turbulence  for  k  >  Icq, 

in  terms  of  an  indeterminate  parameter  k]_,  and  by  means  of  Heisenberg's 

assumptions,  themselves  suggested  by  Weizsacker's  hypothesis  of  similarity 
(reference  k2) ,   which  will  be  examined  later. 

Discussion.-  If  k  is  small  against  kj_,   F(k)   is  reasonably  pro- 
_5 

portional  to  k   .   If  k  is  large,   F(k)  is  proportional  to  k"?   (fig- 
ure 7)-   These  two  laws  have  been  indicated  by  Heisenberg.  The  first  is 
of  special  interest  and  used  for  deriving  the  corresponding  form  of  the 
correlation  functions . 

The  formula  (25-11 )  has  been  established  only  when  k  >  kg  and, 

actually,  has  no  physical  significance  if  k — ^0,  because  F(k)  cannot 
become  infinite.  It  is  supposed  here,  in  first  approximation,  that,  if 
k  <  kQ,  the  true  form  of  F(k)   is  F(k)  =  0,  and,  moreover,  that  k 

decidedly  is  smaller  than  k-^.   To  interpret  this  hypothesis  mathematically 
is  the  same  as  replacing  k^  by  00  in  (25-11);  hence,  one  uses  the 
physical  limiting  conditions  in  which  k]_  is  infinite  (very  high  Rey- 
nolds numbers ) . 

Putting 


F(k)  =  0 


=  F(k0)1/iLN 


if  k  1  kQ 


>         (25-12) 


if  k  >  k 


0 


the  correlation  function  R(r)   can  be  computed. 
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According  to  (1J+—I4-) : 

2  Poo       _5 
R(r)  =  2F(kn)kQ3     sln  rkk  3dk  (25-13) 

J*o   rk 


This  integral  is  convergent  when  Kq  >  0;    it  is  not  derivable  twice 

under  the  integration  sign.  An  attempt  is  made  to  find  the  form  of  func- 
tion R(r)  for  small  values  of  r.  It  is  shown  that  there  exists  a  num- 
ber a  such  that 


R(r)  =  3uq2  (l  -  ArA 


A     being  a  constant, 


When  r — >0,  it  is  easily  verified  that  R(r)   has  a  finite  limit, 
equal  to 


3u02  =  2F(k0)k03  /   k  3dk  =  3F(k0)ko         (25-U) 
Jk0 


Forming 


3^2  -  R(r)  =  2F(k0)k03  /  (l   -   £l^k"3dk       (25-15) 


and  making  the  change  of  the  variable  rk  =  x  in  the  integral,  leaves 

r*f"  (l  -^Mx"3dx  (25-16) 
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Posting 


I  = 


sin  x 


dx 


The  integral   I     differs  from  I  by  infinitesimally  small 

quantities  with  r.   Hence,  in  first  approximation  when  r  is  not  too 
great, 


R(r)  =  3uo2ll  -  Ar3 


(25-17) 


with 


A  =  ^f'l   =  0.801  kQ 


The  exponent  looked  for,  a,  is  equal  to  — .  This  formula  is  impor- 

3 
tant,  and  we  shall  give  later  on  demonstrations  for  it,  of  wholly  differ- 
ent appearances.  The  formula  that  gives  the  correlation  function  f(r) 

2 
is  similar  to  (25-17)-   It  again  is  a  "law  in  r^." 

It  is  a  limiting  law  for  infinite  kj_.   If  k-j_  is  finite,  form- 
ula (25 -11)  yields  for  R(r)   a  law  of  the  ordinary  type 


R(r)  =  3UQ2  1 


2A2 


which  proceeds,  for  r  =  0,  with  a  horizontal  tangent  and  a  finite  radius 
of  curvature  equal  to  A.  When  k^ — ^°°,  this  radius  of  curvature  tends 

toward  zero  and  the  limiting  curve  has  the  form  (25-17)  (figure  8). 


106  NACA  TM  1377 

Note.-  Comparison  of  Heisenberg's  equations  (25-5)  and  (24-18). 
Equation  (25-5)  is  an  equation  in  F,  which,  theoretically,  identifies 
F  and  can  he  used  to  reduce  it  to  a  (complicated)  differential  equa- 
tion. The  calculation  was  made  in  the  "steady"  state.   It  leaves,  in  the 
expression  of  F,  the  constant  C,  which  equation  (25-5)  can  obviously 
not  determine.   If  (24-18)  had  been  an  equation  in  F  too,  it  could  have 
been  used  to  determine  C.  But  the  procedure  that  results  in  this  equa- 
tion and  whose  original  purpose  is  to  eliminate  the  triple  correlations, 
meets  the  expectations  only  partially,  because  it  replaces  in  some  way 
the  triple  correlations  by  the  correlations  of  space  and  of  time .   Equa- 
tion (24-18)  serves  therefore,  as  far  as  the  mathematical  problem  is 
accessible,  to  compute  the  form  of  the  space -time  relations,  once  the 
space  correlations  are  determined  by  (25-5)-  But  it  does  not  appear  that 
the  numerical  solution  of  (24-18)  has  ever  been  attempted  up  to  now. 


CHAPTER  IV 


THEORY  OF  LOCAL  IS0TR0FY  AND  STATISTICAL  EQUILIBRIUM 

26.  Introduction: 

The  importance  attached  to  the  study  of  isotropic  turbulence  was,  up 
to  the  last  few  years,  justified  by  considerations  of,  seemingly,  essen- 
tially practical  value. 

The  development  of  turbulence  toward  isotropy  had  in  its  favor  a 
rather  feeble  theoretical  argument,  namely,  that,  far  from  the  sources  of 
turbulence,  the  Reynolds  stresses  could  be  assumed  zero  (see  section  6), 
which  is  a  condition  necessary,  but  far  from  being  sufficient,  for 
isotropy. 

Batchelor  (reference  12)  proved  in  1948  that  this  argument  is  of  no 
value j  the  elements  of  the  anisotropy  of  turbulence  during  its  formation 
could  be  recovered  in  the  final  phase  of  its  evolution. 

But,  from  the  experimental  point  of  view,  the  isotropy  of  turbulence 
behind  a  grid  seems  well  established,  by  numerous  experimental  works,  as 
will  be  shown  later.  But  none  of  these  evolved  a  general  concept;  it 
seemed  that  a  particular  case  was  studied  here  solely  by  reason  of  its 
simplicity. 

The  theories  to  be  examined  here  form,  by  way  of  contrast,  an  unusual 
ensemble,  not  only  because  of  the  successful  experimental  verifications 
but  particularly  because  of  the  scope  in  which  they  contribute  to  further 
studies  of  turbulence.  They  bring  to  light  the  characters  of  turbulent 
motion  which  are  of  general  significance. 
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The  most  outstanding  result  is  perhaps  the  better  understanding  of 
the  role  played  by  the  viscosity  in  turbulence.  The  viscosity  is 
involved  at  the  birth  of  turbulence  and  at  its  death;  at  its  birth  it 
is  the  condition  necessary  for  the  creation  of  the  rotation  and  the 
diffusion  of  large  turbulent  masses  in  the  fluid;  it  intervenes  at  its 
death  (decay)  since  the  turbulent  energy  is  ultimately  transformed 
into  heat;  but  the  complexity  of  the  turbulent  motion  stems  precisely 
from  the  fact  that  the  viscosity  plays  a  negligible  part  during  the 
major  part  of  the  transformations  of  turbulent  energy. 

The  fundamental  outline  is  as  follows :   relatively  large  turbulent 
masses  with  relatively  low  velocity  gradients  are  formed  in  the  creative 
regions  of  turbulence.  The  sense  of  the  word  relative  can  be  defined  by 

p 

Uf-vLp,  u-^ 

introducing  the  Reynolds  number   u   ,  where   p-^—  is  the  mean  kinetic 

v  2 

energy  and  Lq  the  mean  quantity  of  these  turbulent  masses.  2  In  order 
to  have  turbulence  in  the  proper  sense  of  the  word,  the  ratio 


Uq      V 

Lo 

must  be  high.   If  it  is  small,  we  find  ourselves  at  once  in  what  Batchelor 
has  called  the  final  period  of  turbulence,  in  which  the  rotation  decays, 
so  to  speak,  on  the  spot.  The  forces  of  inertia  play  a  secondary  part. 
The  flow  is  quasi -laminar.  This  is  the  phase  through  which  all  turbulence 
passes .   It  has  been  experimentally  and  theoretically  studied  by  Batchelor 
and  Townsend  (references  10,  12),  and  we  shall  come  back  to  it  later  on. 

If  the  Reynolds  number  is  high,  the  vortices  undergo  a  double  evolu- 
tion, according  to  Heisenberg's  representation  (reference  23);  they  get 
bigger  on  the  one  hand  and  smaller  on  the  other,  that  is,  their  spatial 
dimensions  increase  on  the  one  hand,  and  break  up  into  smaller  vortices 
on  the  other.   Leaving  aside  the  first  aspect  of  this  evolution,  for  the 
time  being,  the  theory  of  similarity  studies  the  development  of  the  energy 
passed  on  to  the  small  vortices;  the  latter  undergo  a  development  similar 
to  that  of  the  generating  vortices,  they  pass  on  their  energy  to  smaller 
vortices  and  so  forth  up  to  the  moment  where  the  size  of  the  vortices  is 
such  that  their  Reynolds  number  is  small.   A  bigger  and  bigger  part  of  the 
energy  is  then  changed  into  heat.   This  somewhat  vague  representation  is 
defined  a  little  more  accurately  further  on,  in  particular,  by  the  analysis 
of  the  rotational  development. 

This  rather  vague  expression  is  specified  to  some  extent  by  the 
study  of  dynamic  equilibrium  (chapter  V). 
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It  is  reasonable  to  think  -  and  this  assumption  will  be  backed  up 
by  its  consequences  -  that  the  development  of  vortices  is  the  more  rapid 
the  smaller  their  size.   On  the  other  hand,  the  statistical  character  of 
the  velocity  of  the  fluctuations,  which  is  prescribed  by  experience,  is 
assumed. 

Hence  it  may  be  imagined  that,  when  the  energy  has  been  transmitted 
to  vortices  sufficiently  smaller  than  the  initial  vortices,  the  structure 
of  these  vortices  is,  on  the  average,  independent  of  that  of  the  generating 
vortices,  and  that  the  relatively  slow  development  of  the  latter  has  no 
effect  on  the  development  of  the  smaller  size  vortices.   The  latter,  there- 
fore, find  themselves  in  a  quasi -balanced  state,  one  determining  factor  of 
which  is  the  energy  per  unit  mass  and  time  which  is  supplied  by  the  bigger 
to  the  smaller  vortices. 


27.  Definition  of  local  homogeneity  and  local  isotropy: 

At  first  we  shall  follow  Kolmogorof f ' s  method  (reference  31.   For 
more  details  see  Batchelor,  reference  6)  which  seems  the  most  direct  as 
well  as  the  easiest,  at  least,  to  begin  with.  Kolmogorof f  introduces 
quantities  which  can  be  completely  defined  in  a  limited  spatial  domain  G 
and  a  limited  temporal  field  T.   Assume  that  M  and  M'   are  two  points 
of  G,  that  t  and  t'   are  two  instants  of  t;  we  call  r  the  dis- 
tance MM ' ,  and  put  t  =  t  -  t ' . 

The  difference  w  of  the  velocities  at  M  and  M'   at  instants  t 
and  t'   is  introduced 


wi 


u1(M,t)  -  u^MSf)  (27-1) 


and  defined  by  its  law  of  probability  in  the  space-time  domain  (G,T). 
One  notices  immediately  that  w  is  little  dependent  on  the  spectral 
components  of  the  velocity  which  develop  slowly  in  the  spatial  domain  G 
and  in  the  time  interval  T. 

By  its  definition  w  depends  a  priori  on  M  and  t,  as  well  as  on 
M '   and  t ' . 

Kolmogoroff  defines  what  he  calls  the  local  homogeneity  and  local 
isotropy.   A  condensed  version  of  these  definitions  reads  as  follows: 

The  turbulence  is  said  to  be  locally  homogeneous  when  the  laws  of 
probability  for  the  statistical  function  w  depend,  in  the  domain  (G,T), 
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only  on  t  and  the  vector  MM1   (actually,  instead  of  this  vector,  the 
vector  e  defined  later  on  is  involved) . 

The  turbulence  is  termed  locally  isotropic  when  it  is  locally  homo- 
geneous and  when  these  probability  laws  depend  neither  on  the  orientation 
nor  on  the  sense  of  the  coordinate  axes . 

According  to  the  foregoing,  it  may  be  assumed  that  in  every  turbu- 
lence having  a  rather  high  Reynolds  number,  there  exists  a  spatial 
domain  G  which  is  small  compared  to  the  size  of  the  big  vortices,  and 
a  time  interval  T,  short  compared  to  the  time  necessary  for  their  energy 
to  change  appreciably,  and  such  that  the  turbulence  is  locally  homogeneous 
and  isotropic  in  the  domain  (G,T). 

But  these  definitions  can  have  no  physical  significance  unless  the 
reference  system  used  by  the  observer  is  itself  physically  coordinated 
with  the  turbulence.   The  domain  G  must,  in  some  fashion,  follow  the 
turbulence  with  the  mean  flow  velocity.  Kolmogoroff,  in  his  original 
report  (ref .  31)  defines  the  velocity  w  about  a  reference  point   (Mqj^o) 

in  terms  of  the  components  g^  of  the  vector  MqM  and  of  t  =  t  -  tQ 
by  the  formulas 


w,- 


±   =  u±(M,t)  -  ujM^to)  (27-2) 


then  he  introduces  a  vector  e  having  as  component 


-i  =  S±  -  t%(mo^o)  (27-3) 


This  is  a  random  vector  whose  physical  significance  is  evident  in 
figure  9.   Since  the  essential  point  of  the  vector  Tu(MQ,t0)  arises  from 

the  "mean  velocity"  in  the  domain  G,  the  use  of  the  coordinates  e. 
is  tantamount  to  following  the  turbulence  in  its  entire  motion. 

The  correlation  measurements  in  time  relative  to  the  laws  of  proba- 
bility of  w  are  few  and  far  between.   Some  measurements  have  been  made 
in  water  and  in  air,  but  none  in  the  wind  tunnel,  as  far  as  we  know. 
Obviously,  in  such  measurements,  time -space  intervals  must  be  combined 
and  it  must  be  ascertained  whether  correlations  close  to  1  can  be 

MM' 
obtained  for  a  finite  value  of  the  ratio  which,  if  it  exists,  define 

T 

the  mean  velocity  at  the  scale  involved. 
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But,  at  the  present  state  of  research  the  study  is  limited,  to  simul- 
taneous quantities  for  defining  the  statistical  values  w.   As  in  ordinary, 
isotropic  turbulence,  it  is  usually  limited  to  the  moments  of  the  second 
and  third  order.  The  tensor  of  the  moments  of  the  second  order 


WifMjMoWMSMQ 


simultaneously  introduces  three  points,  but  can  be  expressed,  on  account 
of  the  isotropy,  homogeneity  and  incompressibility,  as  a  sum  of  moments 

relative  to  the  couples   (M,Mq  i,   (M',Mq),   (m,M' 


w^M^ojwj^M'^oj    =  i  w^MjMoJvj^Moj 


r±  (M  '  ,Mq)  v j  ^M  '  ,Mq)    -  v±  (M,M  'J  w j  (u ,M ') 


(27-4) 


In  this   case,   the  tensor 


Bij   =  wi(M,M')Wj(M,M)'    =    (u'±   -  ujju'j    -  Uj) 


plays  a  role  parallel  to  that  of  the  tensor  R-j.   in  isotropic  turbulence, 
and  by  the  same  argument  it  is  apparent  that  the  isotropy  leads  to  writing 


Bij 


Mill 


LBdd(r)    -  Bnn(r)J    +   5ijBnn(r)  (2T"5) 


where 


*Wr)   -■    fu'd  -  ud^2 


)2  Bnn(^)   "    (*'»  -  ^ 
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ud  and  1^  representing  the  velocity  components  along  MM'   or  normal 
to  MM' . 

A  total  moment  can  be  defined  too: 


»<r)  =  JZ\±  =    IKu-  -u±)2 


The  incompressibility  involves  the  relation 


dB 
Bnn  =  Bdd  +  |  -^  (27 -6) 


equivalent  to  Karman's  relation  in  isotropic  turbulence  (9-7)-   Other 
relations  can  be  established,  in  particular,  the  expression  of  the 
dissipation  of  energy  by  viscosity 


.  .  i5v«l  (27.7) 

2  V  **  Jr-0 

Likewise,  the  moments  of  the  third  order  can  be  expressed  by  the 
single  function 


Bddd(r)  =  (u'd  "  ud)5 


28.  Similarity  hypotheses.   Statistical  equilibrium: 

The  definition  of  these  local  quantities  enables  the  notion  of  sta- 
tistical equilibrium  of  small  vortices  to  be  defined.   According  to  our 
hypotheses,  this  equilibrium  is  not  dependent  on  the  particular  charac- 
teristics of  the  flow .  The  only  quantities  which  can  be  of  influence 
are: 

(1)  The  energy  supplied  per  unit  mass  and  time  to  the  small  vortices 

by  the  larger  ones; 

(2)  The  factor  governing  the  kinematic  viscosity. 
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If  the  fluid  is  assumed  incompressible,  the  specific  mass  can  be 
discounted. 

Hence  Kolmogorof f ' s  postulate  of  similarity:  in  a  domain  where  the 
turbulence  is  locally  homogeneous  and  isotropic,  the  laws  of  probability 
depend  only  on  r,   e,   and   v. 

This  postulate  justifies  the  dimensional  analysis  which  results  in 
setting 


Vr>=   V^Pdd(f)  (28-1) 


In  this  formula,   (3^  is  a  universal  function,  and  I   =  I  — 

represents  a  length,  called  the  local  scale  of  turbulence  by  Kolmogorof f. 
This  length  is  much  smaller  than  the  length  of  dissipation  A  of  the 
isotropic  turbulence.   It  is,  in  effect  (compare  ( 1^-11 )) 


e  =  15-^- 
A2 


whence 


A 

I 


^K.kS 


unA 
where  the  Reynolds  number  E-,    =  — ^—  must  be  at  least  some  tens  If  the 

— A     y 

conditions  of  similarity  are  to  be  realized. 

The  nondimensional  statistical  quantities  must  be  universal  func- 

r 

tions  of  — .   Consequently 

Z 


NACA  TM  1377 


113 


/r\     fu' 
(coefficient  of  dissymmetry)    S, f— 1  =  — — 


(coefficient  of  flatness) 


d  "  ud 


(u'd  "  ud)' 


-.)' 


ud)' 


(28-2) 


and  similar  quantities  Sn,  A^,   with  the  normal  velocity  component  un. 

They  are  universal  functions  in  a  certain  domain  2  «  r  «  L,  in 
which  L  is  a  length  characterizing  the  size  of  G.   It  is  to  be  noted 
that,  accoruing  to  the  hypotheses,   2  and  L  can  be  slowly  variable 
functions  of  time.  The  precise  definition  of  the  range  of  validity  of 
the  formulas  is  ticklish,  but  in  any  case,  if  r  is  made  =  0,  the 
values  of 


S,,(0)  = 


cV3 


bx 


and 


Afl(0) 


(28-3) 


must  be  universal  constants,  which  affords  a  first  confirmation  of  the 
similarity  hypothesis . 

Townsend  (ref.  37  >  38),  at  Cambridge,  developed  a  method  of  measuring 
S(0)  and  A(0).   The  measurements  were  first  made  in  isotropic  turbulence 
produced  by  a  grid.   These  coefficients  are  constant  over  a  wide  range  of 
Reynolds  numbers  (ratios  varying  from  1  to  5)  and  at  varying  distances 
from  the  grid  (likewise  in  a  ratio  near  to  5);  that  is,  at  different 
moments  of  the  development  of  the  turbulence . 

Townsend  quotes  the  value  3.k9     t  0 .Ok   for  A  and  -O.38  for  S,  with 
a  little  lower  accuracy.   It  should  be  recalled  that  for  Gaussian  distri- 
butions the  values  of  A  and  S  are,  respectively,  3  and  0. 
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To  check  the  theory  of  locally  isotropic  turbulence,  by  which  the 
structure  of  the  small  vortices  is  independent  of  the  nature  of  the  large 
vortices  and  the  formation  of  turbulence,  Townsend  made  also  measurements 
in  the  wake  of  a  cylinder  (refs.  39  >  ^0,  and.  4l),  that  is,  in  turbulence 
neither  homogeneous  nor  isotropic  (in  the  original  sense).   He  determined 
thus  Sd(0),   Sn(0),   Ai(°)  and  An(°)   over  "the  entire  extent  of  the 
wake,  at  distances  of  more  than  80  diameters  from  the  cylinder.   The 
values  obtained  are  the  same,  with  sufficient  accuracy,  as  those  obtained 
for  isotropic  turbulence  behind  a  grid.   Furthermore,  the  measurement  of 
the  quantities 


u-]_  being  parallel  to  the  overall  velocity  (axis  of  x^_),  shows  that,  as 

for  isotropic  turbulence  and  throughout  the  entire  wake,  these  quantities 
satisfy  the  equations  of  isotropy 


'^ul\  l/^A  l/^u3\ 


Sx   /  2\dx    /         2\3x 


\ 


By  the   same  argument,   the  expression  of  the  energy  dissipation  is 


13  falY 

— v( J  ,  as  for  isotropic  turbulence,  and  as  the  equations  (27-5)  and 

2  \dx  / 

(27-7)  niay  lead  to  presume. 

These  formulas  are  very  important  because  they  show  the  possibility 
of  defining  the  dissipation  of  energy  in  any  turbulent  flow  by  measure- 
ments bearing  on  a  single  velocity  component. 

29.  Case  of  high  Reynolds  numbers  (reference  6): 

The  study  of  the  statistical  equilibrium  of  vortices  can  be  extended 
further  when  the  Reynolds  number  is  high  enough  so  that  the  viscous  dis- 
sipation in  the  largest  vortices  of  the  statistical  equilibrium  is  rela- 
tively low;  it  can  be  assumed  then  that  the  statistical  equilibrium  of 
these  large  vortices  is  not  dependent  on  v. 
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This  is  Kolmogorof f ' s  second  hypothesis: 

If   Z  «  r  «  L,  the  laws  of  probability  regarding  w  are  then 
dependent  on  e     and  r  only. 

This  implies  that  a  domain  G  can  be  found  whose  size  L  is  large 
compared  to   Z,  that  is,  that  the  Reynolds  number  is  very  high. 

Dimensional  analysis  permits  then  the  definition  of  the  form  of  the 
moments  of  the  statistical  variable  w.   In  fact  it  is  readily  apparent 

that,  for  B,,(r)   not  to  be  dependent  on   v,   3(— )  must  be  of  the  form 


p(f)  -  cte  *  f f 


whence 


2 
BrtH(r)  =  C(er)3  (29-1) 


C  being  a  universal  constant  (figure  10). 

The  moment  B    must  have  a  similar  form.   The  equation  of  continuity 

indicates  then  that  B   (r)  =  — B^^r). 

r 
On  the  other  hand,  when  —  is  small, 

Z 


Bnn(r)  =  |B"nn(°)r2  =  2Bdd(r)  (29-2) 


Likewise,  one  must  have 

Bddd(r)  =  Cte  X  er 
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But,  if  r  is  small 

Bddd(r)  =  Cte  X  r5  (29-3) 

-,3 


"because  in  that  case,   ^ddd(r)  4   Sd(0) 


Bdd(^) 


p 
,  in  first  approximation. 


Since  I      depends  on  v,  the  nondimensional  quantities  such  as  Sd(  — 

must  be  constant,  for  values  of  r  much  greater  than  I.      On  the  other 
hand,   Sd(0)   is  also  a  universal  constant,  and  the  experiment  seems  to 

indicate  that  these  two  constants  are  identical. 

To  tell  the  truth,  there  is  only  one  series  of  measurements  of  S,(r) 

made  by  Townsend  in  a  region  where  the  turbulence  was  isotropic,  by 
measuring  the  triple  correlation  function 


c(r)  =  &L 
3 


u 


0 


The  relations 


u')3  =  _6u<)3c(r)      (u  -  u')2  =  Uq2  [l  -  f (rji     (29-4) 


are  utilized. 


The  experimental  points  are  placed  near  Sd  =  -O.38  for  values  of 

—  of  the  order  of  kOl .     The  absolute  value  of  S  decreases  slowly 
I 

afterward . 

Instead  of  arguing  about  the  correlation  functions,  one  may  just  as 
well  apply  the  dimensional  analysis  to  the  spectral  functions  F(k)  for 
values  of  k  which  are  sufficiently  large  without  being  too  large,  how- 
ever.  Since  F(k)  must  be  dependent  only  on  k  and  e,  it  is  found  that 

2  _2 
F(k)  must  be  proportional  to  e\.   3  in  the  domain  under  consideration. 
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.1 
This  is  the  spectral  law  in  k  5  spoken  of  earlier  during  the  discussion 
of  Heisenberg's  theory  (section  25). 

Using  the  detailed  hypotheses  of  Heisenberg,  we  have  then  shown  that 

2 
there  is  equivalence  between  the  law  of  correlation  in  r^  and  the  spec- 

_5 
tral  law  in  k  -> .      It  is  now  seen  that  this  equivalence  was  inevitable, 
and  that  it  is  not  the  result  of  the  particular  form  of  the  law  of  energy 
transfer  in  the  spectrum,  but  solely  due  to  the  existence  of  a  general 
similarity  hypothesis  among  the  hypotheses  the  utilization  of  which  was 
unimportant . 

By  using  the  equations  of  motion  the  value  of  Bddd(r)   can  be  defined 

and  a  relation  established  between  the  constants  C  and  S,  which  is  sus- 
ceptible to  experimental  verifications . 

The  equation  corresponding  to  that  of  Karman  (l8-5)  reads  here 

■**   =  B'ddd  +  ^ddd  "  6v(B"dd  +  ^'ddj  (29-5) 


The  time  factor  does  not  appear,  because  the  hypothesis  of  local 

homogeneity  contains  that  of  steadiness  with  time,  so  that  — uq   must 

dt 

be  taken  equal  to  -  ^e . 

3 

The  equation  (29-5)  multiplied  by  r4",  is  integrated  in  the  form 


6vB'dd  -  Bddd  "  ^  (29-6) 


If  r  is  small,  Bddd  can  ^e  neglected,  so  that  one  finds  again 
the  relation 


*  =  ^vB"  „(0) 


ddv 
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For  the  values  of  r  in  the  interval  Z  «  r  «  L,  B, , ,  prevails 


over  B'dd,  because  Bddd  is  proportional  to  r  and  B'dd  to  r 
Hence 


1 
3 


Bddd  =  ~*r 


k 

—t 

5 


But  in  (29-3)  it  was  seen  that  in  first  approximation: 


Bddd(r)  =  Sd(0)[Bdd(r)j2  =  Sd(0)c26r 


Consequently 

C=l--^—f  (29-7) 

Therefore,   C  can  be  calculated,  if  the  value  of  Sd(0)   is  known 
and,  particularly,  if  the  experimental  value  Sd(0)  =  -O.38  is  assumed. 
In  that  case 

2         %  2 


Bdd(r)  =  c(er)3  =  1.64 (er) 5 


Conversely,  the  experimental  study  of  the  function  B,j(r)  enables 

C  to  be  measured  and  the  obtained  value  to  be  compared  with  the  value 
computed  from  Sd(0). 

30.  Validity  of  the  similarity  laws: 

All  second-order  moment  measurements  have  been  made  in  isotropic 
turbulence.   Hence  it  is  useful  to  examine  what  becomes  of  the  theory  of 
local  isotropy  if  applied  to  isotropic  turbulence  proper. 


NACA  TM  1377  119 


The  first  point  is  the  determination  of  the  space -time  range  of 
validity.  The  discussion,  necessarily  rather  vague,  has  been  made  "by 
Batchelor  (ref.  6).   It  is  briefly  summarized  here: 

From  the  spatial  point  of  view,  isotropy  appears  to  be  assured  for 

J  "loo 
f  dr. 
0 

The  difficulty  arises  from  the  estimation  of  the  time  interval.   An 
evaluation  of  the  characteristic  period  of  the  vortices  indicates  that 
the  ratio  between  the  period  of  the  largest  and  that  of  the  smallest 

vortices  is  of  the  order  of  ^,  where  A  is  the  length  of  dissipation. 

A 
This  ratio  must,  therefore,  be  great.   Experience  indicates  an  order  of 

magnitude  of  —  in  what  Batchelor  calls  the  "initial  phase  of  decay  of 

A 
turbulence."   (it  is  to  be  defined  in  chapter  V.)   For  the  present,  it 
is  simply  stated  that  it  is  produced  by  the  turbulence  behind  a  grid  of 

mesh  size  M,  when  the  Reynolds  number  R,  =  —  of  the  grid  (U  =  mean 

velocity)  is  sufficiently  high  (1,000  to  300,000)  and  the  measurements 
are  made  at  a  distance  from  the  grid  of  the  order  of  1,000  M.   One  finds 

T  S."X  Uf)^ 

then  that  —  is  of  the  order  of  — -,  where  R-\  =  -^—  may  be  called  the 
A  10        ^    v 

Reynolds  number  of  turbulence.  During  the  phase  in  question,   R^  is 

/1m 

constant  and  equal  to  v  /  — ,  where  A  is  a  number  dependent  on  the 

V  A 
constitution  of  the  grid. 

For  values  of  R,   not  quite  up  to  200,  which  still  gives  only  values 
—A 

of  ratio  -  of  the  order  of  20.  Townsend  has  measured  the  correlation 

A 


function  g(r).   Since 

1  _  g(r)  =  ^miil 

2UQ2 


1  -  g(r) 

must  be  constant. 

2 
r3 
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The  curves  plotted  for  various  values  of  Fix*  and  —  have  similar 

-*1      M 

shapes.   They  present  a  rather  flat  maximum  between  r  =  0.1L  and 

r  =  0.3L,  from  which  the  value  of  the  constant  C   (section  29)  can  he 

deduced.   The  average  of  the  obtained  values  is  1.53>  with  discrepancies 

of  less  than  3  percent.  Considering  the  accuracy  of  the  measurements  of 

Sd(0),  the  agreement  with  the  value  1.64  deduced  from  S^(0)   is 

sufficient. 

Other  measurements,  made  at  a  lower  Reynolds  number,  yield  C  =  1.33> 
while  those  by  Dryden  produced  C  =  I.50.   In  conclusion,  although  the 
range  of  a  validity  of  the  law  is  too  restricted  for  the  aspect  of  the 
curves  to  be  convincing,  the  agreement  of  the  values  of  C  is  quite 
satisfactory. 

The  lower  limit  of  validity  of  the  law  of  local  isotropy  in  these 
experiments  is  of  the  order  of  magnitude  of  A,  that  is,  well  superior  to 

A    hi — 
2,  since  —  =  <yl5R^  is  here  of  the  order  of  20  to  25- 

L 

It  seems  that  the  Reynolds  numbers  used  are  still  too  low  for  the 
field  of  application  of  Kolmogorof f ' s  second  hypothesis  to  be  extended 
very  far.   The  first  hypothesis,  less  restrictive,  is  much  better  verified. 
Townsend's  correlation  measurements  show,  in  fact,  that  the  curves  repre- 
senting g(r)   are  when  r  «  A  functions  of  —  only,  at  different  dis- 

A 

tances  from  the  grid.   The  similarity  is  extended  to  values  of  —  higher 

A 

than  unity  for  the  highest  values  of  R^s,  where  the  second  hypothesis 

possesses  a  domain  of  validity.   This  validity  is  tied  to  the  fact  that 
IK   is  constant  at  the  beginning  of  the  turbulent  development.   (Compare 

chapter  V.)  This  information  on  the  development  of  turbulence  can, 
obviously,  not  be  deduced  from  the  theory  of  local  isotropy  which  assumes 
e  constant. 

On  the  other  hand,  the  absence  of  similarity  in  the  correlation 

curves  at  high  values  of  —  suggests  that  the  range  of  application  of 

A 

the  similarity  law  is  very  restricted,  and  also  belief  in  the  fact  that 

the  ratio  of  correlation  length  L  to  dissipation  length  A  varies  in 

reasonable  fashion  during  the  evolution  of  turbulence,  even  when  R^ 

T   Ea 
remains  reasonably  constant,  because  the  relation  —  =  — —  is  only  a 

A   10 
first  approximation. 

In  reality,  the  domain  of  statistical  equilibrium  is  extremely 
extended,  but  this  does  not  appear  in  the  structure  of  the  spatial  corre- 
lation tensor.   For  a  clear  definition  of  the  individual  behavior  of  each 
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scale  of  vortex  dimensions  the  introduction  of  the  spectral  tensor  is 
necessary.   Heisenberg's  theory  is  therefore  found  capable  of  exceeding 
the  stage  reached  by  the  theory  of  local  isotropy,  and  it  likewise 
permits  interpretation  of  the  different  stages  of  development  of  turbu- 
lence discussed  in  paragraph  31>  and  in  chapter  V. 


31.  Interpretation  of  the  laws  of  statistical  equilibrium  in  spectral 
terms  -  Weizsacker's  and  Heisenberg's  theories  (refs.  k-2  and  23 , 
respectively) : 

Kolmogorof f ' s  theory  assumes  a  locally  steady  equilibrium  and  its 
spatial  domain  of  application  is  necessarily  restricted.   However,  exami- 
nation of  the  correlation  curves  and  of  the  laws  of  variation  of  the 
correlation  length  L  during  the  decay  of  isotropic  turbulence  behind  a 
grid  leads  one  to  believe  that  the  scope  of  similarity  will  be  extended 
farther  than  anticipated. 


But  this  notion  of  similarity  can  not  appear  clearly  unless  the 
influence  of  the  different  vortices  is  separated.  This  means  changing 
from  the  correlation  functions  to  spectral  functions . 

It  will  be  remembered  (compare  sections  11  and  12)  that  to  the 


spatial  correlation  tensor 
tensor  cpag(A)   defined  by 


RapU) 


uau 


o  there  corresponds  the  spectral 


W6) 


1 
/<Pa0(?v)e 


P=l 


Vi 


dA 


If  there  is  isotropy,  the  spectral  tensor  depends  on  a  single  spectral 
function  F(k),  just  as  the  spatial  correlation  tensor  depends  on  a  single 
function  R(r)   for  example. 

To  the  properties  of  R(r)  for  small  values  of  r  correspond  those 
of  F(k)   for  great  values  of  k,  as  results  from  formula  (14-^)  which 
links  F  to  R,  and  which  is  recalled  here: 


R(r)  =  2 


sin  r 
rk 


^F(k) 


dk 
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It  has  been  shown  (section  29)  that  to  the  law  of  correlation  in 
2 

rJ  for  R(r),  applicable  in  the  domain  I  «  r  «  L,  and  resulting  from 
Kolmogorof f ' s  second  hypothesis,  there  corresponds  a  spectral  law  of  the 
form 

F(k)  =  A£5k  5 


where  A  is  a  universal  constant.   This  law  is  valid  in  the  domain 
kQ  <  k<  k   and  the  numbers   2k_^   Lie-,  are  of  the  order  of  unity. 

In  the  most  general  case  of  Kolmogorof f ' s  first  hypothesis,  one  can 
only  put : 

F(k)  =  F0*(£) 

for  k  >  kg,  Fq  and  kQ  being  functions  of  e  and  v  only  and  <J>  a 
universal  function. 

Pk 
Putting     f(k')dk'  =  T(t)   in  (22-8),  gives 

Jo 

§-  /  F(k')dk'  =  -T(k)  -  v    2k,2F(k')dk*         (31-1) 

at  Jo         "     Jo 


The  function  T(k)  represents  the  transfer  of  energy  by  turbulence 
of  the  large  eddies  to  smaller  eddies. 

In  the  case  in  which  the  dissipation  by  viscosity  is  negligible, 
the  conditions  of  statistical  equilibrium  are  written  simply  T(k)  =  c^e  . 

3 
5 

A  dimensional  analysis  affords  the  spectral  law  in  k   .   In  fact,  it  can 

.  .  OUn 

be  said  that  T(k)  has  the  dimensions  of  averages  so  that  u-iUot-^  is 
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v  3k,  where  v^.  Is  a  certain  characteristic  velocity  of  vortices  having 

k  for  wave  number.   Hence  one  must  have   v^3k  =  c^e,  in  the  domain 

_1 
ks  «  k«  kg,  and  v^  is  proportional  to  k  5. 


For  velocity  v^,  the  quantity  J\     (u1^  -  u^p  which,  according 

2 

to  Kolmogorof f ' s  theory,  is  proportional  to  r*  in  the  chosen  conditions, 
can  be  assumed.  Weizsacker's  law  is  equivalent  to  Kolmogorof f ' s  law 

2  1 

because  F(k),  having  the  dimensions  of  — — ,  is  proportional  to  k  *> . 

k 

(Compare  section  25-) 

The  variation  with  k  of  other  quantities  associated  with  different 
vortices  can  also  be  determined.   For  example,  the  characteristic  period 

of  the  vortices  is  proportional  to  -i— ,  a  quantity  which  has  the  dimen- 


kvk 


2 


sions  of  a  time  interval,  that  is,  to  k   .   It  decreases  when,  k 
increases,  which  is  in  agreement  with  the  hypothesis  according  to  which 
the  characteristic  periods  of  development  of  the  small  vortices  are 
shorter  than  those  of  large  size  vortices. 

It  can  likewise  be  verified  that  the  viscous  dissipation,  which  is 

1 

o  5 

proportional  to  k  F,  varies  like  k  .  Hence,  it  cannot  be  neglected  for 

great  values  of  k. 

If  now  the  viscosity  is  no  longer  neglected,  the  form  of  the  func- 
tion T_(k)  must  be  specified.   This  is  what  Heisenberg  did  who,  as 
shown  in  section  25,  has  assumed  that  T(k)  could  be  put  in  the  same 
form  as  the  losses  of  energy  by  viscosity,  the  form  of  the  "turbulent 
viscosity"  coefficient  resulting  from  similarity  considerations. 


3 

The  comparison  with  experiment  of  the  spectral  law  k^  has  not  as 

yet  been  carried  out  on  flows  in  which  the  Reynolds  number  is  high  enough 
so  that  the  conclusions  are  well  defined.   Nevertheless,  its  range  of 

2 
validity  seems  more  extended  than  that  of  the  law  r^  for  the  correla- 
tions, because  the  influence  of  the  various  vortex  sizes  is  different  due 
to  the  use  of  the  spectral  analysis. 
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Townsend  attempted  to  check  the  validity  of  the  expression  of  the 
function  T(k)   given  by  Heisenberg,  and  according  to  which  F  should  be 

_7 
proportional  to  k  '   at  great  values  of  k.   The  mean  value  of 


is  proportional  to   I  k°F(k)dk.   If  the  law  in  k~'  were  verified,  this 

Jo 

integral  should  be  infinite,  which  signifies  that  the  integral 

Ik,, 

k°F(k)dk  should  increase  considerably  together  with  its  upper 


2 


limit  kg.   However,  the  experiment  seems  to  indicate  that  this  integral 
has  a  limiting  value  independent  of  kp,  and  variable  during  the  decay  of 

turbulence.   Therefore,  the  law  in  k~7  cannot  be  true  for  the  highest 
values  of  k,  which  limits  the  range  of  validity  of  Heisenberg's  theory. 


CHAPTER  V 


DECAY  OF  THE  TURBULENCE  BEHIND  A  GRID 


32.    History: 

The  theory  of  local  isotropy  shows  that  the  laws  of  turbulent  flows 
are  the  same,  regardless  of  the  origin  of  turbulence,  provided  that  the 
conditions  are  "local"  and  the  Reynolds  numbers  high  enough.   This  chapter 
deals  with  the  particular  laws  of  a  turbulent  flow  of  great  importance, 
that  which  is  produced  in  the  test  chamber  of  wind  tunnels,  and  in  which 
turbulence  is  largely  due  to  a  grid,  that  is,  an  obstacle  of  periodic 
structure.   Turbulence  arises  from  the  mixing  of  wakes  of  grid  elements, 
and  is  dissipated  progressively  as  it  moves  farther  away  from  the  grid. 
For  an  observer  carried  along  with  the  velocity  of  the  ensemble,  the 
turbulence  is  therefore  subjected  to  a  "decay"  as  function  of  the  time. 
The  laws  of  this  decay  are,  naturally,  compatible  with  the  local  laws  of 
Kolmogoroff ,  but  today  they  are  known  with  a  high  degree  of  accuracy  in  a 
much  more  extended  domain,  and  it  can  be  stated  that  it  is  in  the  study 
of  this  phenomenon  that  the  theory  of  turbulence,  steadily  checked  by 
experiment,  has  made  the  fastest  and  most  constructive  advance. 
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The  basis  of  this  theory  was  established  by  G.  I.  Taylor,  who,  by 

considerations  of  similarity,  had  indicated  that,  if  the  Reynolds  number 

1  ?  1 

is  high,  the  turbulent  energy  — u^   must  decrease  as  t    as  a  function 

of  the  time . 

The  measurements  made  since  by  numerous  researchers  have  shown  that 
this  law  is  close  to  experimental  results;  nevertheless  the  experimental 
curves  were  rather  of  the  type 

UQ2tn  =  Cte 


n  ranging  between  1  and  2.  As  the  precision  of  the  measurement  is  low 
as  soon  as  the  turbulence  intensity  becomes  low,  the  interval  to  which 
the  comparisons  refer,  is  restricted. 

The  variation  of  the  length  of  dissipation  A  in  terms  of  time  is 
tied  directly  to  that  of  the  intensity  of  turbulence  by  the  relation  (19-1) 

1  duo2    10v 


uq2  dt     A2 


2  -n 

If  u^   varies  as  t   ,  one  deduces  from  it : 


A2  =  ^t 
n 


But,  up  to  the  last  few  years,  the  measurements  of  A  offered  no 
possibility  of  defining  the  value  of  n.   The  value  of  A  is,  in  effect, 
deduced  from  the  curvature  at  the  origin  of  the  curve  representing  the 
correlation  function  f(r).   The  measurements  are  made  with  two  anemom- 
eters, of  which  the  distance  r  cannot  be  reduced  indefinitely  for 
manifold  reasons  (geometrical  and  physical).   As  a  result,  the  curve  f(r) 
is  not  well  known  for  small  values  of  r. 

However,  an  examination  of  the  correlation  curves  seems  to  indicate 
that  their  shape  remains  constant  in  terms  of  the  distance  from  the  grill, 
during  the  decay  of  turbulence. 

All  these  questions  have  taken  a  decisive  forward  step,  as  a  result  of 
A.  A.  Townsend's  new  methods  which  enable  direct  measurement,  by  electrical 
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methods,  of  the  momentary  or  mean  values  of  the  quantities  u,   — ,   

^     ax2 

and  those  of  the  products  of  two  or  three  quantities,  as  well  as  the 
direct  determination  of  the  statistical  frequency  curve  of  a  quantity  in 
terms  of  the  velocity. 

The  use  of  these  methods  enabled  G.  K.  Batchelor  and  A.  A.  Townsend 
to  determine  the  law  of  decay  of  turbulence  produced  by  a  grill  with  a 
good  degree  of  accuracy  over  a  wide  range  of  values  of  the  Reynolds  num- 
bers of  the  grill,  defined  by 


*M  '  f  (32-X) 


where  M  is  the  width  of  the  mesh,   U  the  overall  velocity. 
When  the  turbulence  is  isotropic 


J^2=  _10v^!=  "10v(^r  (32-2) 


dt  U        A2       \c 


du-^2 


X 


The  direct  measurement  of  [ j   affords  the  evaluation  of  A 

much  more  quickly  and  more  accurately  than  by  the  procedure  employed 
previously,  and  consequently  defines  A  in  terms  of  time,  that  is,  of 
the  distance  from  the  turbulence -producing  grid. 

53-  Initial  and  final  phase  of  turbulence  (Batchelor  and  Townsend,  refer- 
ences 8  and  11) : 

Townsend 's  measurements  made  it  possible  to  distinguish  in  the  decay 
of  turbulence  an  initial  phase  and  a  final  phase,  separated  by  a  less 
well  defined  intermediary  phase . 

The  initial  phase,  clearly  proved  by  these  experiments,  corresponds 
to  grill  distances  of  less  than  100  or  I50  mesh  widths  (variable  according 
to  the  Reynolds  number).   It  is  characterized  by  the  fact  that  the  decay  of 

turbulent  energy  can  be  represented  by  an  expression  proportional  to  t"n 

where  the  exponent  n  differs  from  unity  by  less  than  10  percent.   It  is 

2  -1 

said  that  Uq   decreases  as  t 
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In  the  final  phase,  which  corresponds  to  the  limiting  state  of  decay 

_5 
o  p 

of  turbulence,  experiments  prove  that  u_   decreases  as  t 

The  initial  phase  is  also  characterized  by  the  fact  that  the  Reynolds 
number  of  turbulence,  defined  by 


V 

v 


RA  =  -°-  (35-I) 


remains  constant  throughout  the  period  of  decay  to  which  it  corresponds , 

If,  in  fact,   un2  =  At"n,  it  is  seen  that  A2  =  i^t,  and  that 

n 

1-n 
lOAv  +  2 


^K   = 

Thus  Ra   remains  constant  if  n  =  1,  and  in  this  case  only. 

To  define  the  structure  of  the  two  essential  phases  of  turbulence, 
we  now  introduce  the  rotation  (Batchelor  and  Townsend,  reference  10)  of 
which  the  relations  with  the  correlation  functions  and  the  dissipation 
length  A  are  simple,  at  least  when  the  turbulence  is  isotropic.   The 
general  equations  relevant  to  the  rotation  are  simplified  by  the  Karman- 
Howarth  equation. 

The  rotation  of  the  velocity,  defined  by  the  formulas 

Sux   cftio  cSu-i    du*  duo   Sut 

^  =  _2  -  J*     ^  =  ^_^     ^.J*.^:   (35_2) 

0x2   oxv  0X3   ox-l  ox-j_   0x2 


satisfies  two  well  known  equations  which  are  obtained  by  simple  combina- 
tions of  the  equations  of  motion 

Sui   V~  bu±       1  dp 

—  +  2_  ui^—  =  -  -—  +  v  Au±  (33-3) 

ot     1        lSxz    P  oXj_ 
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If   (33-3)    is   derived  with  respect  to     x^,    since  one   interchanges 
i     and     j      and  subtracts   the   two  equations   obtained,    the   result  is 


dm,-        v —       doL-        v —       cW 

^  +  Z_  V"1  =  Z_  %P~  +  v  ^  (33-^) 

ot  k         dxk         k  ox-^ 


The  two  members  of  (33-^ )  ^xe   multiplied  by  2(x>±,    or  added  with 
respect  to  i,  with  due  allowance  for  incompressibility,  the  averages 


taken  and  homogeneous  turbulence  assumed.   Putting  3u>  =  7   «i?    leaves 


2   V-   2 


3^-o)2  =  2  23  ^-"V—  +  2v  XI  ^—^|         (33-5) 

st  ik    x  ^^         ik     axk^ 

If,  in  addition,  the  turbulence  is  isotropic,  then  (section  9): 

2 

1 


or  -■-   a>]_  =  ay>  =  0^3  =  5-^ 

\2 


In  this  equation,  the  first  term  of  the  second  member  represents 
the  increment  of  the  vorticity  due  to  the  elongation  of  the  vortex 
filaments,  and  the  second  the  decrement  due  to  the  influence  of  the 
viscosity.   It  is  verified  that  this  last  term  really  is  negative  because 


5  oij_   d  oij_ 

"**  2  "  v  2~  ' 
c*k    3xk 

and  —  is  zero  if  the  turbulance  is  homogeneous. 

Sxk 

Batchelor  and  Townsend  measured  the  three  terms  of  (33-5)  directly. 
But  it  is  convenient  to  give  equation  (33-5)  a  slightly  different  and 
easier-to-interpret  aspect.   The  Karman-Howarth  equation  (l8-5)  can  be 
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used  to  evaluate  the  first  member.   Introducing  the  limited  developments 
of  f,   c,  up  to  the  terms  in  r^  (cf .  section  9) 


f(r)  =  1  +  iif'(o)  +  i(]°(o)  +  .  •  . 

2        2k 


t   \         3  c'"(0) 
:(r)  =  r^  — - — -  + 

6 


and  identifying  the  terms  in  r2,  gives 


Wfl"{0) 


^vuo2f(^(0)  +  luoV'^O) 


(33-6) 


Now,  the  derivatives  of  the  functions  f  and  c  can  be  expressed 
in  terms  of  the  average  values  bearing  on  the  velocity  derivatives.   The 
calculation,  as  done  before  for  the  rotation,  gives 


w? 


V  f"(°)  =  -kr1 


Sx- 


U02*(10(0) 


^2Ul\2 


u05c'"(0) 


,  SXn 


If  therefore  the  nondimensional  quantities 


(33-7) 


s  =  - 


fe)J 

J-               "1 

3 

[fe)1 

2 

G  =  Uq 


;-|2 


/au,\2 


(33-8) 
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are  introduced  of  which  the  first  has  already  played  a  prominent  part  in 
section  28  (coefficient  of  dissymmetry,  for  r  =  0),  it  is  seen  that 


c»  ,  =  -aV'(o)    g  =  f(j°(°>  =  xM^(o)   (33-9) 

f"2(0) 


s  = 

3  -   ^..2 

2 


-f"(0) 


A  being,  as  will  be  remembered,  a  "length  of  dissipation"  defined  by 

f(o)  =-JL. 

A2 

—    „  2 

Since  olt  =  5-^— >  equation  (33-6)  can  be  written 
A2 


.3 


i^  ,  X(^2s  _  |£  (33-IO) 

dt    3\/5     \   V 


where  R>  =  —^—  is  the  Reynolds  number  of  turbulence  already  defined  by 


UqA 
v 
equation  (33-2). 


So  in  turbulence  language,  equation  (33-10)  is  a  consequence  of  the 
Karman-Howarth  equation  of  turbulence,  equivalent  to  (33-5),  but  superior 
for  comparisons  with  experiments . 

The  measurements  of  the  various  terms  of  (33-IO)  during  the  initial 
period  of  decay  of  turbulence  confirm,  first,  that 


*te       ^2 


R  =  Cze  Ad  =  10 vt  (33-11) 


t  being  counted  from  any  convenient  origin. 
Consequently  one  may  write: 


u2    ,  R2 


of   =  5-2-  =  _  =L.  (33-12) 


A2    20  t2 
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and  (33-10)  takes  the  form 


2  =  |*as  +  y  (33-13) 


It  has  been  shown  (section  28)  that  the  measurements  of  the  coeffi- 
cient of  asymmetry  S  prove  that  the  latter  is  independent  of  R,   and 

close  to  O.38.  The  direct  measurements  of  G  are  in  good  agreement  with 
the  above  formula,  in  which  S  =  O.38,  which  constitutes  a  check  on  the 
qualities  of  the  measurements. 

These  results  make  it  possible  to  evaluate  the  ratio  of  the  two 
terms  of  the  second  member  of  the  equations  (33-5)  or  (33-10).  The  quo- 
tient of  the  term  representing  the  dissipation  through  viscosity  to  that 
which  expresses  the  inertia  effects  is  equal  to 


^i-  =  1  +  -i°_  (33-1M 

SRA      7SRA 
or,  replacing  S,  by  its  experimental  value 


1  +  §£  (33-15) 


h 


This  ratio  is  constant  throughout  the  initial  period  of  decay  of 
turbulence.   If  —  is  considerably  less  than  unity,  the  contributions 


hs 


-n  2 
am 


to  of  the  viscosity  and  of  the  inertia  terms  are  of  comparable 

bt 
importance  and  the  velocity  resulting  from  the  variation  of  the  vorticity 
is  of  an  order  of  magnitude  lower  than  each  one  of  these  two  terms.   This 
is  the  range  of  high  Reynolds  numbers  R,  ,  where  the  theory  of  statistical 

equilibrium  (Kolmogoroff ,  Weizsacker)  applies.  For  large  enough  values  of 

_5 
the  wave  number  k,  the  spectral  law  k  3     can  t)e  observed  there .   The 
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measurements  show,  in  fact,  that,  in  proportion  as  R.   increases,  the 

correlation  curves  approach  the  limiting  shape  with  tangent  vertical  to 
the  origin  as  implied  by  the  formula 


2 

f  =  1  -  cr5 


The  original  parabolical  region  becomes  consistently  smaller  as  R. 
becomes  higher  (Figure  11;  compare  also  the  end  of  section  25)- 

If  —  approaches  unity  (Reynolds  numbers  neither  too  high  nor  too 
^A 
low),  the  results  deteriorate  and  new  partial  or  total  similarity  hypoth- 
eses must  be  made,  resulting  in  solutions  of  the  type  of  those  studied  in 
section  21  from  the  purely  mathematical  point  of  view,  and  which  are  not 
repeated  here. 

op 

If  —  is  great  (low  Reynolds  numbers),  the  effect  of  the  triple 

^A 
correlations  is  negligible  against  the  viscous  dissipation.   The  mechanism 
of  the  decay  of  turbulence  is  mainly  due  to  the  viscosity,  and  the  solu- 
tions of  section  20  can  be  invoked.   That  is  what  happens  in  particular  in 
the  final  phase  of  decay  which  is  discussed  below. 


3^4- .  Concepts  regarding  the  structure  of  the  final  phase  of  turbulence 
(Batchelor  and  Townsend,  reference  12;  Batchelor,  reference  8): 

It  is  necessary  to  predict  by  theory  the  law  of  the  final  phase  of 

5 

2  "5" 

decay  of  turbulence,  according  to  which  Uq   decreases  as  t   ,  and, 

consequently 

A2  =  kvt  (3^-1) 


When  the  turbulence  is  sufficiently  attenuated,  the  fluctuation  is 
slight,  and  the  inertia  effects  negligible  against  the  effects  of  vis- 
cosity. As  a  result,  the  nonlinear  terms  in  the  equations  of  motion  can 
be  discounted,  and  one  may  simply  write 

Sui   1  op 

— -  = £—  +  v  Au-j 

St    P  bx± 
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In  these  conditions  the  equations  of  motion  can  be  integrated,  the 
velocity  components  computed  in  terms  of  time  and  the  initial  conditions, 
and  the  correlations  deduced.  This  is  the  program  set  forth  in  section  16, 
and  it  cannot  be  accomplished  except  in  very  special  cases .   Thus  the 
time  correlations  can  be  calculated  just  as  well  as  the  classical  space 
correlations . 

Only  the  case  of  isotropic  turbulence  is  analyzed.  To  neglect  the 
inertia  terms  is  to  neglect  the  triple  correlations  in  the  "fundamental 
equations."  Thus  one  reverts  to  the  conditions  of  section  20,  or  turns  to 
the  possibility  of  expressing  the  spectral  equation  in  the  simplified  form 


^  +   2vk2F  =  0  (3^-2) 

St 


The  general  solution  of  this  equation  reads 


J2(- 


P(k,t)  =  F(k,t1)e'2vk  ?'**)  (34-3) 

t-^  being  the  instant  that  makes  the  start  of  the  final  phase. 

One  may  also  use  equation  (20-1),  equivalent  in  correlation  terms 
to  the  spectral  equation  (3^-2),  which,  with  the  triple  correlations 
disregarded,  reads 

|£  =  2v  AR 
St 

of  which  the  general  solution  (20-9)  is  a  little  harder  to  write. 

But  these  solutions  must  have  an  asymptotic  character.   It  is  of 
no  interest  here,  except  when  t  -  t-j_  is  sufficiently  great.   F(k,t) 

approaches  zero  when  t  -  t-j_  approaches  infinity,  as  is  natural,  since 

the  turbulence  decays  progressively.   So,  when  t  -  t-,   is  great,  the 

only  regions  of  the  spectrum  continuing  to  exist  are  those  for  which  k 

is  small,  k.     t  -  t-])  remaining  practically  finite.   It  has  been  seen 

(in  section  18)  that,  for  small  values  of  k,   F(k)  was  necessarily  of 

the  form  Ck  ,   C  being  a  constant  independent  of  the  time,  and  tied 
to  the  Loitsiansky  invariant. 
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Consequently,  for  the  large  values  of  t,  the  spectral  function 
takes  the  form 

,    ,     k  -2vk2ft-t1')  ,   ' 

F(k,t)  =  Ck^e     ^   -1/  (3k-k) 

valid  for  small  values  of  k,  that  is,  for  large-size  vortices. 

Example  1  in  section  20  gave  the  corresponding  form  of  the  correla- 
tion functions  R(r,t)   and  f(r,t).   In  particular  it  is  recalled  that 

r2 


f(r,t)  =eMt-tl) 


and  that 


5  JO         52  ^2 


vft  -  tj) 


.1 
2 


This  is  the  law  of  decay  that  had  to  be  explained.  The  proof  rests 
on  the  hypotheses  of  section  20  and  section  18,  particularly  on  the  fact 
that  the  Loitsiansky  invariant  is  finite  and  other  than  zero.  The  agree- 
ment of  theory  and  experiment  confirms  them. 

In  this  calculation,  it  was  supposed  that  the  turbulence  was  isotro- 
pic.  On  the  other  hand,  only  the  large  size  vortices  continued  to  exist 
because  the  energy  contained  in  the  region  of  the  spectrum  corresponding 
to  substantial  values  of  k  becomes  negligible.   As  the  constant  C  is 
independent  of  time,  the  shape  of  the  spectrum  for  small  values  of  k  is 
independent  of  the  state  of  decay  of  turbulence.   The  large-size  vortices 
contain  very  little  energy,  but,  in  the  final  phase,  the  energy  of  the 
small  vortices  is  dissipated  and  only  the  large  vortices  remain  visible 
whose  energy,  as  feeble  as  it  is,  has  become  very  great  proportionally. 

Now,  while  the  turbulence  is  always  locally  isotropic  and  exhibits 
a  tendency  to  isotropy  for  the  small  vortices  (large  values  of  k),  it  is 
no  longer  the  same  for  the  large -size  vortices,  in  which  the  geometrical 
dissymmetries  of  the  wind  tunnel  are  exhibited. 
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A  resumption  of  the  calculations,  but  without  introducing  the  hypoth- 
esis of  isotropy,  produces  similar  conclusions;  the  portion  of  the  spec- 
trum (that  is,  of  the  spectral  tensor)  relative  to  small  values  of  k 
remains  constant  during  the  decay,  while  the  rest  of  the  spectrum  vanishes 
progressively.   It  follows  that  the  anisotropy  of  turbulence  that  results 
from  the  geometrical  and  mechanical  configuration  of  the  flow,  and  which 
is  practically  concealed  at  the  beginning  of  the  decay  of  turbulence,  must 
appear  in  the  final  phase  (Batchelor,  ref.  9)-   Unfortunately,  it  is  rather 
difficult  to  check  it,  since  the  turbulence  is  then  very  weak. 


35. 


The  concept  of  "dynamic  statistical  equilibrium"  (Heisenberg, 
Batchelor,  ref.  8): 


ref.  23, 


In  the  problem  of  the  decay  of  turbulence  behind  a  grill,  the  study 
of  the  shapes  of  spectral  curves  and  correlation  curves  led  Heisenberg  to 
believe  that  the  notion  of  similarity  could  be  extended  to  vortices  con- 
taining the  major  part  of  the  turbulent  energy.  The  overall  spectrum, 
with  exception  of  the  region  of  small  wave  numbers  or  large  vortices 
(spectrum  in  k^"),  thus  would  be  in  a  dynamic  statistical  equilibrium  or 
quasi -equilibrium,  in  which  the  energy  distribution  in  the  spectrum  is 
modified  during  the  decay  of  turbulence,  the  maximum  shifting  toward  the 
small  wave  numbers,  but  the  general  shape  of  the  curve  remaining  the  same, 

According  to  the  static  theory  of  local  isotropy,  the  spectral  func- 
tion F(k)   can  be  written  in  the  form 


F(k) 


Wit 


(35-D 


where  —  is  a  length  that  fixes  the  transition  between  the  spectral 

ks 
region  of  predominant  inertia  forces  and  that  where  the  viscosity  forces 
become  comparable  to  them.   0  is  a  nondimensional  function.  Fq  and  k 

are  dependent  on  e  and  v  only,  and  dimensional  analysis  shows  that 


51 
FQ  - 


vTe 


1  I 
ks  =  v  V> 


(35-2) 


It  was  shown  (in  section  29)  that,  when  the  viscosity  is  negligible, 

_5 
.F  becomes  proportional  to  k  J,   which,  in  this  limiting  case,  fixes 
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the  form  of  the  function  0.   In  the  general  case,  the  form  of  the  func- 
tion $  depends  upon  that  of  the  function  $  (k)   or  T(k)  which  appears 
in  the  second  term  of  the  fundamental  equation  (22-8)  or  (31-1),  and  which 
represents  the  transfer  of  energy  through  turbulence  from  the  large  to  the 

v 
small  size  vortices.   <D  is  in  all  cases  a  universal  function  of  -=-,   and 

the  similarity  law  which  it  defines  is  valid  in  a  domain  k^  <  k  <  kg, 
where  lu  is  an  approximate  lower  limit.   This  domain  is  the  more  extended 

Uq 

the  higher  the  Reynolds  number  R~,  =  -^—-     This  number  characterizes  the 

"°   kov 
importance  of  the  portion  of  the  spectrum  in  absolute  equilibrium.   This 
portion  is  maintained  during  the  decay  of  turbulence,  so  that  Rq  remains 

constant,  as  well  as  the  nondimensional  ratio  -^— .   It  is  found  that 


lev 

this  number  intervenes  when  the  region  of  absolute  equilibrium  is 
left.   The  quasi -equilibrium  therefore  depends  on  the  parameters   e  and 

3~  2 

2 

which  the  dissipation  e  is  the  derivative  with  respect  to  time .   By 


v,  but,  in  addition,  on  a  third  parameter,  the  total  energy  E  =  — Uq  ,  of 


u  ^  u„A 


is  proportional  to  the  Reynolds  number  R,  =  ~^- 
\fev~  -A    v 


of  turbulence.   For  a  quasi -equilibrium  to  be  possible,   R>   must  there- 
fore be  constant  during  the  development  of  turbulence .   This  is  what  the 
experiment  proves  during  the  initial  phase  of  decay  of  turbulence,  and 
this  result  supplies  an  argument  in  favor  of  the  quasi-equilibrium  hypoth- 
esis.  If  R,   varies,  dynamic  similarity  ceases,  but  that  does  not  pre- 
vent the  possible  existence  of  a  region  of  local  isotropy,  variable  during 
the  development . 

By  virtue  of  the  known  formula  (ik- 11) 

2      ur>2 


2  2 

The  parameters   v,   €  and  u~.   can  be  replaced  by  v,   A  and  u^  .   If 

UqA 

R-*  =  is  constant  in  a  defined  regime  of  decay,  the  spectral  func- 

— a    v 

tion  F(k,t)  depends  only  on  a  variable  parameter,   A  for  example,  equi- 
valent to  the  time,  and  in  addition,  proportional  to  Jt. 
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We  write 


F(k,t)  =  F(k1,t)$(r1)  (35-5) 

where  r\  -  —     is  a  nondimensional  parameter  dependent  on  the  fixed  wave 

kl 
number  k]_.   The  nondimensional  function  $  is  completely  defined  by  the 
value  of  R^,  that  is,  by  the  initial  conditions  at  the  grid.   The  length 

of  dissipation  A  can  serve  as  reference  length  — .  The  dimensional 

kl 
analysis  proves  then  that  F(k]_,t)   is  the  product  resulting  from  uq2A, 

or  — — ,  a  quantity  proportional  to  —  multiplied  by  a  quantity  con- 

A  ft 

stant  during  the  decay.   Introducing  a  velocity  vQ  and  a  reference 

length  — ,  one  may  put: 


« - kA  =  k\/?    '*'*>  ■  ^fe f  «">       <55-4> 


The  form  of  the  function  <J>  is  now  defined  by  means  of  the  funda- 
mental spectral  equation  (22-8),  which  we  recall  here: 

^  +  2vk2F  =  I  (35-5) 

St 


We  examine  first  what  happens  in  case  of  small  values  of  parameter  t\. 
We  shall  show  that   <j>(t})   is,  according  to  (35_10  >  infinitely  small  of  the 
first  order  with  respect  to  tj,  and  that  the  position  of  the  tangent  at 
the  origin  to  the  spectral  curves  F(k,t)  is  independent  of  t.   For  the 
sake  of  simplification 


F(k,t)  =  -i-a^Ti!)       ^  =  k\/t 
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By  (35-5) 

1 
r\19\   -  4L  +  Wt^2^  =   2t2f  (35-6) 

Assuming  that  $-j_  allows  a  development  of  the  form 

,  /   x       n       n+1       n+2 

I^Ul)  =  oqTj-l  +  o-xTVl    +  c^tj-l    +  .  .  . 

the  first  member  of  (35-6)  becomes,  except  for  the  factor  t  2, 


(n  -  1)oqTi111  +  na1Tl1n+     +      (n  +  1)02  +  ^va^ 


f+2  +  .  .  .  (35-7) 


It  was  seen  in  section  22,  that  $(k)  was  infinitely  small  as  k  , 


that  is,  as  i\    .      It  follows  that,  either 


n  =  1      a,]_  =  0      (n  +  l)cLg  +   ^voq  =  0,  .  .  . 
where  the  first  term  of  (35-7)  must  be  the  term  in  t\   ;   or  else: 


n  =  6 


In  the  hypothesis  n  =  6,  there  would  result,  for  F(k,t)  a  law  k° 
which  is  incompatible  with  the  exact  law  k^"  and  with  experiment. 

For  n  =  1,  there  is  for  Ffat)  a  law  in  k  which  is  not  compatible 

with  experiment  when  k — >0     but  is  verified  for  small  enough  values  of  k, 
except  at  the  limit.   In  that  case 
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and 


F(k,t)  =  o^k  +  a£tk5  +  .  .  .  (35-8) 

Oq,   cip  being  constants  independent  of  t.   It  follows  that 


f  =  °0  +  5a2tk2  + 


approaches  a  limit  independent  of  t  when  k — >0.      All  spectral  curves 
compatible  with  the  quasi -equilibrium  have  the  same  tangent  at  the  origin. 

Therefore,  if  the  part  of  the  spectrum  in  k   (represented  by- 
dotted  curves  in  figure  12)  is  neglected,  the  spectral  curves  are  approxi- 
mately homothetic  with  respect  to  the  origin,  the  maximum  decreasing  and 
shifting  toward  the  small  values  of  k  when  the  time  increases,  that  is, 
with  progressing  decay  of  turbulence. 

The  similarity  is  obtained  again  on  the  correlation  curves  f  and 
g.   In  fact,  (Ik-k) 


J  sin  -r— Tin 

"loo  poo         !—   '-L 

£HLi^(k,t)dk  =  %         l^1(Tll)dTll   (35-9) 

0  fc  Jo   _£_ni 


^t 


The  functions  Rt,  and  hence  }      f,      g,  depend  therefore  on  the 

uo2 

variable  — ^—  and  not  separately  on  r  and  t.   In  other  words,  the 
shape  of  the  correlation  curves  is  maintained,  provided  that  -^—  or  — 
is  taken  as  variable. 
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Incidentally,  —  can  be  taken  as  variable  also,  because 


dk  = 


2uo2 


P"»l(Tl1)  d^ 
Jo   ^1    s[t 

(35-10) 

Since  in  the  region  where  R-,  =  C   ,   Uq   is  inversely  proportional 

to  t,   L  is  proportional  to  v/t~  as  A,  and  —  is  constant. 

A 

But  the  similarity  with  L  is  more  quickly  destroyed  during  the  decay 
of  turbulence  than  the  similarity  with  A,  because  in  the  expression  of  L, 
the  values  of  $-,   corresponding  to  the  small  values  of  r\-,      are  of  greater 

2  P    te 

weight  than  the  others .   The  law  Uq  A =  c    prevails  therefore  much 

2  2    te 
longer  than  the  law  Uq  L  =  c   .   The  explanation  is  that  similarity  is 

approached  with  respect  to  L,  while  it  still  seems  to  be  verified  for  A. 
It  can  be  seen  that  the  similarity  discrepancies  of  the  correlation  curves 
appear,  first,  for  the  great  values  of  r,  and  that  they  afterward  reach 
the  central  zone .   The  reason  for  this  zone  to  start  changing  is  that  the 

maximum  energy  of  the  spectral  curve  reaches  the  spectral  region  in  k  . 
The  law  of  decrease  of  energy  can  no  longer  be  of  the  form  Uq  t  =  C*e . 

The  turbulence  reaches  the  final  period  of  decay,  in  which  the  energy  UQ^t 

_5 
decreases  as  t   . 

In  this  whole  analysis  it  is  not  necessary  to  assume  that  R\   is 

high.   The  main  point  is  that  the  initial  conditions  are  such  that  the 

region,  in  which  the  spectral  function  F  varies  as  k\  be  relatively 
small.   If  R,   is  very  high,  there  exists  a  region  of  the  spectral  curve 

_5 
in  which  F  varies  as  k  3.  ^h±s   region  does  not  exist  when  Ik   is 
medium. 

If  it  is  desired  to  go  farther  and  define  the  form  of  the  function  F, 
the  form  of  the  energy  transfer  $(k)   (Heisenberg,  reference  2k)   must  be 
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chosen.   If  $(k)   has  the  form  given  to  it  by  Heisenberg  (section  25), 
for  instance 


Pk 


T(k)  = 


Jo 


jCk'Jdk1  =  20     /F(k')  dk'    k,2F(k')dk'   (35-11) 


k  V  k 


,3 


the  reduced  variables  (35-^)  are  employed,  and  it  is  proved  that  (35-5) 
leads  to  the  reduced  equation 


2>(V)dTT  -  — ti* ( tj )  =  2 


A  +  C 


MvidV 


V2»(rr)dV    (35-12) 


where  A  =  v—  is  the  inverse  of  a  Reynolds  number. 
v0 

Equation  (35-12)  is  easily  reduced  to  a  differential  equation  which, 
contrary  to  what  occurs  in  the  case  of  the  steadiness  hypothesis  (25-6), 
seems  unsolvable.   Heisenberg  (reference  2k)   has  pointed  out  approximate 
solutions,  valid  for  great  or  small  values  of  A.   (The  reader  is  referred 
to  Heisenberg 's  report  and  also  to  Batchelor,  reference  8). 

36.  Synthesis  of  the  results  relating  to  the  structure  of  the  spectrum 
of  turbulence: 

Obviously  only  isotropic  turbulence  is  involved.   The  results  col- 
lected regarding  Heisenberg 's  spectral  function  F(k)  which,  incidentally, 
cannot  be  measured  directly,  but  is  deduced,  by  simple  transformations, 
from  the  directly  measurable  spectral  functions  (spectral  functions  of 
Taylor,  insofar  as  they  define  a  spectral  distribution  in  space,  and  not 
in  time )  or  from  the  correlation  functions . 

Whatever  the  structure  of  turbulence  may  be,  the  spectral  laws  tend, 
when  the  Reynolds  number  becomes  high  enough,  toward  an  absolute  limiting 

law,  where  F(k)   is  proportional  to  k  3   (absolute  equilibrium).   The 


corresponding  correlation  functions  are  of  the  form  f(r)  =  1  -  C  r*. 
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In  the  problem  of  the  decay  of  turbulence  behind  a  grid,  an  initial, 
an  intermediate  and  a  final  period  are  distinguished.   For  small  values  of 

k,   F(k)   is  always  =  Ck\   C  being  a  constant  independent  of  time 
(bound  to  the  invariant  of  Loitsiansky) . 

For  the  average  values  of  k,  and  in  the  initial  period,  a  state  of 
"dynamic  equilibrium"  exists  (figure  13)- 

The  corresponding  portions  of  the  spectral  curves  are  homothetic  with 

respect  to  the  origin  and,  when  the  time  of  decay  increases,  their  maximum 

approaches  the  origin.   Their  shape  is  known  up  to  the  origin,  but  the 

straight  part  that  corresponds  to  the  small  values  of  k  and  which  is 

fixed  during  the  decay  of  turbulence,  has  no  physical  reality  and  must  be 

k 
replaced  by  the  part  of  the  spectrum  in  k  . 

If  the  Reynolds  number  is  high  enough,  the  portion  of  the  spectrum 
corresponding  to  great  values  of  k  is  defined  by  the  laws  of  absolute 

5N 


equilibrium  \spectrum  in  k  J.     When  k — >°°  the  spectrum  terminates  in 

_5 
a  curve  k~7   (experimentally  doubtful).   The  region  k  5  disappears 
if  the  Reynolds  number  is  not  very  high. 

If  the  Reynolds  number  is  small  enough,  the  region  in  dynamical 

-2k  vt 
equilibrium  is  numerically  known:   F(k)   is  proportional  to  ke~     , 

o 
provided  that  k  vt  is  not  too  great. 

In  this  initial  period  the  energy  of  the  fluctuations  varies  as   t"  . 

In  the  final  period  the  zone  of  dynamic  similarity  disappears,  and 

the  maximum  of  the  spectral  curve  intrudes  on  the  region  k^"   (fig- 
ure Ik).      Starting  from  an  instant  t-j_  marking  the  beginning  of  this 

phase  we  can  write 

F(k)  =  CkV21^-^ 

_5 
and  the  energy  of  the  fluctuations  decreases  as   (t  -  t-A      . 
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Index  of  Principal  Notations,  Fundamental 

Formulas,  and  Dimensional  Equations 

Certain  letters  have  been  employed  successively,  in  various  paragraphs, 
to  represent  different  quantities  without  having  to  be  afraid  of  confusion. 
Others,  on  the  contrary,  always  represent  the  same  quantity;  at  least  from 
chapter  II  on.  We  shall  enumerate  them  and  recall  the  principal  formulas 
into  which  they  enter  and  to  which  one  has  frequently  to  refer : 

X"l,X2,x^  represent  the  coordinates  of  a  point  M 

|  ,|  ,%  designate  the  differences  of  the  coordinates  of  two 

points  M,  M',  of  which  the  distance  is  r.   (in 
the  first  chapter,  through  a  few  lines,  r  repre- 
sents, by  exception,  a  correlation  coefficient.) 

f(r),g(r)  designate  Karman's  double  correlation  functions. 

(in  chapter  I,   f  is  used  to  represent  certain 
statistical  frequencies  and  probability  densities.) 

These  functions  are  derived  from  the  correlation  ten- 

sor  Rop(^i^2,63,t)   ln  space>  of  which  R  =  /   Rqq 
^         '  a 

is  the  scalar  invariant.   (In  the  first  chapter,   R 

denotes  the  density  of  probability  of  position  and 

velocity  of  a  statistical  point.)  R  should  not  be 

confused  with  R  which  denotes  a  Reynolds  number 

particularized  by  a  subscript:   FL,  R, ,  .  .  . 

a(r) ,b(r) ,c(r)      denote  the  triple  correlation  functions  of  Karman 

(section  9)  which  are  deduced  from  the  triple- 
correlation  tensor  To 

A-j_,A2,A^  are  the  vector  components  of  the  wave  number  (spatial 

frequency)  in  the  space  A  of  the  wave  numbers 

k  is  the  length  of  the  vector  A^ 

F(k)  is  Heisenberg's  spectral  function,  associated  with 

the  spectral  tensor  <p  R(^i>^2'^3) 

E  denotes  the  total  energy  of  the  fluctuation  of  tur- 

bulence per  unit  mass 
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£  denotes  the  energy  dissipated  as  heat  per  unit  mass 

v  denotes  the  coefficient  of  kinematic  viscosity 

The  letter  u  always  designates  a  velocity.   From  chapter  II  on  u-j_,Up,u^ 

represent  the  components  of  the  velocity  of  the  fluctuations  at  the 
point  M,  and,  if  the  turbulence  is  isotropic,   uq  is  the  mean  square 

value  of  Uj_ . 

A  is  the  dissipation  length,   L  the  correlation  length 

h  always  represents  a  statistical  function  with  orthog- 

onal increments 

The  quantities  enumerated  above  have  the  following  dimensional  equations: 

F(k)  =  L5T"2 
E  =  L2T~2 

v  =  L2T_1 

2  -2 

R  =  L  T 


2^-3 


t  =  L^T 


By  way  of  comparison,  the  dimensional  equations  of  classical  mechanical 
quantities  read: 


Velocity  =  LT_1 
Acceleration  =  LT 


Force  =  MLT"2 
Work  =  ML2T~2 


Power  =  ML2T"5 
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There  follows  a  list  of  the  principal  formulas,  referred  to  con- 
stantly in  chapters  IV  and  V.   The  sections  in  which  these  formulas 
appear  for  the  first  time  are  shown  in  parentheses . 


Rap^i,62^5*t)  =  ua(x)up(x') 


=  u. 


(Section  8)  < 


^a^  +  g6aP 


R  -  ZZRaa 
a=l 


TaP7  =  u-a(x)up(x)u7(x  +  |) 


(isotropic  turbulence) 


A  -./- 


(Section  9)   < 


g"(0) 


L  =  f(r)dr 


r  bf 


-  f  +  (incompressibility) 

2  dr 


(Section  13) 


E  "  \£  =    /     ?(*)<& 


e   =  2v         k2F(k)dk 


R(r)    -  2 


sin  r 


'0  rk 


^F(k)dk 


(Section  Ik} 


1     P°° 
F(k)    =  -  rk  sin  rk  R(r)dr 

*  Jo 
A2  =  10v£ 
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(Section  19)  — 
dt 


=  -e 


a  2 

— Uq 

dt 


-lOv 


"0 


(Section  22 )< 


—  +  2vk2F  =  f 
St 


ijr  is  a  function  of  k  and  of  t,  infinitely  small  as  kc 


and  such  that 


$  dk  =  0 


fo 


X* 


u 


The  following  symbols  and  abbreviations  are  used: 

is  the  conjugate  complex  of  X. 

is  the  mean  value  of  u.   Beginning  with  chapter  II  the  mean 
values  refer  to  a  point,  considered  in  the  first  chapter  as 
conditional  averages . 

The  bold-faced  letters  in  the  first  paragraph,  such  as  K,  represent 
operators . 

For  simplifying  the  writing,  notations  such  as  f(u;  x)  are  used  in 
place  of  what  should  be  written  unabbreviatedly  f  ('u-^^jU*;  x-]_,X2,x^\; 

likewise,   If  (u;  x)   is  the  abbreviation  for  I  I  If  (u^u^u^; 
x1,x2,x5\du1  du2  du^. 

The  formulas  are  numbered  in  every  section;  for  example,  form- 
ula (15-4)  is  the  fourth  formula  in  section  15. 


Translated  by  J.  Vanier 
National  Advisory  Committee 
for  Aeronautics 
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APPENDIX 


Some  experimental  results . 

In  the  foregoing  analytical  study  of  the  various  recent  theories  on 
turbulence,  reference  was  made  frequently  to  experimental  results,  mostly 
of  English  origin,  "but  numerical  values  were  given  only  rarely.  The  data 
contained  in  this  appendix  are  intended  to  remedy  this  to  some  extent,  by 
way  of  typical  examples  of  several  correlation  curves  and  spectrum  curves 
obtained  by  direct  observations  in  the  wind  tunnel. 

These  curves  have  been  determined  by  A.  Favre  and  his  collaborators 
at  the  Laboratory  of  Mechanics  of  the  Atmosphere  of  the  Institute  of  Fluid 
Mechanics,  at  Marseille,  for  the  official  French  Aeronautical  Research 
Establishment  (0.  N.  E.  R.  A.).   They  are  copied  from  reports  by 
A.  Favre, ^3  except  the  spectral  curve  .F(k)   (fig.  20)  which  was  calcu- 
lated from  other  data. 

The  turbulent  velocity  was  measured  with  considerably  modified  hot- 
wire instruments  of  the  Datwyler  type  which  had  undergone  important  modi- 
fications at  the  Marseille  laboratory. 

The  time -correlation  curve  was  obtained  by  means  of  a  recording 
device  with  Tolana  magnetic  tape,  remodeled  in  the  laboratory  so  that  it 
could  be  used  for  aerodynamic  purposes.   The  electric  current  from  the 
hot  wire  is  recorded  at  a  point  on  the  magnetic  tape  driven  at  uniform 
forward  speed  V.   Then  the  record  is  read  at  two  points  separated  by  a 
distance  D,  before  which  the  tape  unrolls.  The  currents  obtained  corre- 
spond to  two  turbulence  records  separated  by  a  time  displacement  —  and 

then  conjugated,  as  customary  in  space -correlation  measurement,  for  the 
two  currents  coming  simultaneously  from  two  distinct  hot  wires.   The 
method  also  makes  it  possible  to  measure  the  correlations  with  time  and 
space  displacement.   So,  from  the  experimental  point  of  view,  the  problem 
of  time  correlations  is  solved.   In  this  respect  experiment  is  ahead  of 
theory  which  up  to  now  provides  no  clear  prediction  about  the  time- 
correlation  curves.   It  is  to  be  noted  that  the  correlations  at  a  point 
with  time  displacement  are  more  precise  to  measure  than  the  space  corre- 
lations, because  only  one  hot  wire  is  used,  and  the  time  displacement 
can  be  diminished  as  desired,  with  the  employed  instrumentation.   By  way 
of  contrast,  for  the  space -correlations,  it  is  difficult  to  get  two  per- 
fectly identical  hot  wires,  and  it  is  impossible  to  bring  them  closer 
together  beyond  a  certain  limit. 

13a.  Favre:   Statistical  time -correlation  measurements;  VII. 
International  Congress  of  Applied  Mechanics,  London,  19^8. 

New  time -space  correlation  measurements  downstream  from  a  turbulence 
grid,  0.  N.  E.  R.  A.   Report  No.  8/522-A,  Dec.  31,  19^+9- 
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All  the  curves  reproduced  here  correspond  to  an  overall  speed  of 
12.2  m/sec.   For  the  time -correlation  curve  (fig.  18)  the  intensity  of 

turbulence  was  1.92  x  10  ~2. 

The  grid  producing  the  turbulence  in  figure  l6  was  of  3-25  inch 
(8.3cm)  mesh  width;  this  value  was  chosen  in  order  to  facilitate  the 
comparison  with  H.  L.  Dryden's  measurements.   In  the  other  graphs  a 
grid  of  1  inch  was  involved.   In  all  the  cases  the  spacing  of  the  hot 
wires  from  the  grid  was  kO   meshes. 

Figure  15  represents  a  space-correlation  curve .   The  distance  r 
of  the  hot  wires,  which  were  placed  on  a  single  horizontal  perpendicular 
to  the  tunnel  axis,  serves  as  abscissa,  Karmans's  dimensionless  func- 
tion g(r)   as  ordinate. 

Figure  l6  represents  the  same  function  g(r)   for  a  different  grid. 
The  crosses  correspond  to  a  direct  measurement;  the  round  points  corre- 
spond to  the  same  measurement,  but  after  recording  on  the  time -correlation 
instrument,  and  rectification  with  zero  time  displacement. 

This  produces  a  certain  check  on  the  accuracy  of  the  time -correlation 
measurements.   The  curve  was  extended  up  to  near  200  mm.   It  remained 
reasonably  constant  and  below  the  axis  of  the  r. 

Figure  17  shows  the  representative  curve  of  the  longitudinal  corre- 
lation function  f(r),  computed  from  g(r)  by  Karman's  formula,  and  the 

curve  representative  of  the  function  — (f  +  2g),  proportional  to  the  func- 

3 
tion  R(r),  and  obtained  by  calculation  from  the  curve  f  and  the 
curve  g  of  figure  16 .  The  values  of  R  and  f  are  very  small  and 
negative  starting  from  r  =  120  mm  and  r  =  200  mm,  respectively.   (Com- 
pare formula  (18-12).) 

Figure  18  represents  a  time -correlation  curve.   It  is  seen  to  be 
strongly  negative  for  settings  above  5  milliseconds. 

Figure  19  represents  the  spectral  function  F(n)   of  turbulence 
obtained  by  transformation  of  the  time -correlation  curves.   It  is  there- 
fore the  Fourier  transform  in  cosines  of  the  time -correlation  curve. 

■*-%4ore  precisely,  according  to  Taylor's  formulas,   F(n)   is  defined 

so  that,  if  t   is  the  time  displacement,   /  F(n)  cos  2itnT  dn  is  the 

JO  ~ 

time -correlation  function,  and,  in  particular,  so  that   /  F(n)dn  =  1. 

do  ~ 
In  this  case,  infinity  corresponds  practically  to  1,000  periods  per  second, 
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It  was  limited  to  the  low  frequencies  of  around  16  periods  per  second. 
The  theory,  as  far  as  it  can  be  applied  to  the  spectral  energy  distribu- 
tion curves  in  time  at  a  point,  suggests  that  the  spectral  curve  passes 
through  the  origin.   This  corresponds  to  the  fact  that  the  areas  of  the 
positive  and  negative  parts  of  the  time  correlation  curve  are  equal.   It 
is  verified  approximately  on  figure  18.   This  verification  is  not  alto- 
gether rigorous,  but  it  does  seem  that  for  values  of  more  than  30  milli- 
seconds in  time  displacement  the  curve  becomes  positive  again,  although 
remaining  very  close  to  the  time  axis . 

Figure  20  gives,  computed  from  the  curves  of  figure  17,  the  integral 


,    .   ,  f(r)  +  2g(r)  , 

kr  sin  kr— — - — — L   dp 


equal  to  the  quotient  of  the  spectral  function  F(k)  by  the  constant  3uq2. 

Therefore,  with  suitable  units  as  ordlnates,  it  is  the  representative  curve 
of  F(k)j  k  is  the  inverse  of  a  length,  while  in  figure  19,   n  is  the 
inverse  of  a  time  interval.   A  priori,  Taylor's  spectral  function  F(n), 
which  represents  the  spectral  energy  distribution  in  time,  and  Heisenberg's 
spectral  function  F(k)  which  represents  the  spectral  energy  distribution 
in  space,  are  different. 
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Figure  1.-    Recording  of  the  turbulent  fluctuation  velocity  in  a  wind  tunnel 
of  20  cm  x  30  cm.    Main-flow  velocity:    20  m/sec.    Intensity  of  tur- 

bulence:    5-10     .    The  recording  duration  is  0.03  sec.     (Photograph 
Heubes.) 
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Figure  15.-    Transverse  correlation  function   g(r)    in  space. 
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Figure  16.-    Transverse  correlation  function  g(r)    in  space. 
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Figure  17.-    Longitudinal  correlation  function    f(r).    Total  correlation 

function   M  +  MA. 
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Figure  18.-    Time  correlation  function. 
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Figure  19.-    Taylor's  spectral  function   F(n). 
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Figure  20.-    Spectral  function   F(k).  (Value  of  the  integral 


f-  J     kr  sin  kr  i-^  dr) . 
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